本帖最后由 多多会迈特来布 于 2016-2-11 05:59 编辑

我现在是需要用水力学的一个叫做Mansard&Funike的方法来做波浪反射分析，就是通过得到的三个位置随时间变化的水面高度数据，用MATLAB程序，过程中用到傅里叶变换之后，求得随频率分布的密度图像。

我使用的是MATLAB R2010b  操作系统 64位 Win7

​

我得到的最后密度分布图像是这样的

然而另外一个软件得到的图像是

pastedImage.png (59.04 KB, 下载次数: 3)

2016-2-11 05:26 上传

这两张图是在相同的数据条件下得到的，而第二章的图应该是对的。就是中间一个频率的图很高，然后两边低。但是由于那个软件在学术上来说已经过期了， 所以我导师要求我在matlab里面吧这个程序弄通。

以下是我有疑问的程序部分

N=2^14; %Number of Records

L=2^14; %Transform Length

SmoothingK=2^0;

%Define all the required Mansard & Fuke variables

fd=zeros(L,1); %Frequency division matrix

%Sampling Frequency

Fs=(N-1)/(Times(N,1)-Times(1,1));

%Find the resolution from the sampling frequency

f=Fs/L;

%Auto spectal density for each probe

SP1=pwelch(P1,(L/SmoothingK),[],L)/(Fs/12.566);

SP2=pwelch(P2,(L/SmoothingK),[],L)/(Fs/12.566);

SP3=pwelch(P3,(L/SmoothingK),[],L)/(Fs/12.566);

%Find amplitude spectra

AP1=zeros(L,1);

AP2=zeros(L,1);

AP3=zeros(L,1);

for i = 1:(L/2)

AP1(i,1)=sqrt((SP1(i,1)*f));

AP2(i,1)=sqrt((SP2(i,1)*f));

AP3(i,1)=sqrt((SP3(i,1)*f));

end

%Generate cross spectral densities

CSD12=cpsd(P1,P2,(L/SmoothingK),[],L);

CSD13=cpsd(P1,P3,(L/SmoothingK),[],L);

PH11=zeros(L,1);

PH12=zeros(L,1);

PH13=zeros(L,1);复制代码

可不可以请大神帮我看看有没有什么问题 为什么会出现和第二章图片差距这么大的第一张图片。由于我是MATLAB初学者， 对于傅里叶变换也只是看过网络上的扫盲，不知道错误在哪儿啊？  另外由于程序不是我写的，我是拿到这个程序然后将它修改然后让他能运行处正确结果。其中我对于pwelch后面Fs/12.566 这个部分不是很理解。  Fs是采样频率

程序中的AP1，AP2，AP3就是这张图

QQ截图20160211060145.png (43.12 KB, 下载次数: 6)

2016-2-11 05:59 上传

中公式中的Ap,k，是频率k/T对应的傅里叶参数。

以下是我的MATLAB软件和一篇讲述这个软件用的算法的论文，比较重要的内容我在文中用高亮标出来了。论文和完整程序我都附件在了帖子里，由于数据文件太大了，所以就没有上传，如果需要我也可以提供

谢谢各位大神了

我的QQ292031713  如果有思路也可以加我细聊 谢谢啦！ 最近好急这个

我根据讲述这个程序用的Mansard&Funke的方法的文章里面写的计算方法检查过，后面的步骤基本都是正确的，比较可能出现问题的部分就是我写在上方的程序部分。