Matlab画瀑布图,福利叶变换,频谱图代码.docxVIP

瀑布图t=0:0.01:1; f=1:5;[t,f]=meshgrid(t,f);y=sin(2*pi*f.*t);waterfall(t,f,y);%就这样，你可以照样画葫芦xlabel(t);ylabel(f);zlabel(y);x = 1:24;y = 1:3;[x_grid,y_grid] = meshgrid(x,y);z = sin(2*pi/10*x_grid+pi*y_grid/3);waterfall(x,y,z);colorbar用MATLAB画瀑布图时怎样消隐啊，就是前面的频谱截面要挡住后面的。感觉不消隐的话，看起来很不舒服，没有效果。下面是源程序，大家帮忙看看怎么改，能消隐。org_data=distilldata(data.txt); %读取txt数据文件for j=1:1:8??%取第第j列数据? ?data=org_data(:,j);? ?length=size(data)/4096; %%%数据分段?? ?x=0:1:2047; %频率分辨率? ?WFIntval=1; %y轴间距为 1s? ?for i=1:1:2048? ?? ?y(i)=0;? ?end? ?for i=1:1:length? ?? ?z=data((i*4096-4095):(i*4096));? ?? ?z=z-mean(z); %去直流? ?? ?z=abs(fft(z,4096))/2048;? ?? ?z=z;? ?? ?figure(j);? ?? ?plot3(x,y+WFIntval*i,z(1:2048));? ?? ?xlabel(频率/Hz);? ?? ?ylabel(时间/s);? ?? ?zlabel(幅值/mm/s^2);? ?? ?title(瀑布图);? ?? ?view(30,30);? ?? ?axis([0,2000,-inf,inf,-inf,inf]);? ?? ?hold on? ?end? ?gridend连续的傅里叶变换t=-1:0.01:1;x=5*sin(2*pi*10*t);N=length(x);fx=fft(x);df=100/N;n=0:N/2;f=n*df;subplot(211);plot(t,x);grid; subplot(212);plot(f,abs(fx(n+1))*2/N); grid;离散傅里叶变换loadABC3.txtA=ABC3;X=A(1:3:30000);Fs=2000;Y=fft(X,10000);Pyy = Y.* conj(Y) / 10000;f=2000*(0:5000)/10000;subplot(211);plot(X);grid; subplot(212);plot(f,Pyy(1:5001)/5000); grid;x=1:0.1:10;y=f(x);时域：plot(x,y);y0=fft(y);plot(x,abs(y0));plot(x,angle(y0))追问频谱呢？y0=fft(y);plot(x,abs(y0));plot(x,angle(y0))怎样用matlab画出语音信号的时域波形和频谱图和画出加噪声-3db的白噪声后的波形、频谱，滤波图片符号编号排版地图您提交的参考资料超过50字，请删除参考资料：匿名提交回答预计用到的函数有：wavread();读入wav格式的语音信号fft()快速傅里叶变换plot()绘制二维图形randn() 高斯白噪声示例：x=wavread(file.wav);%读取波形文件获得数据x=x(1:1024);%取前1024点作为处理使用数据fx=fft(x);figure(1);subplot(211);plot(x);subplot(212)plot(abs(fx));snr=.3;x1=x+snr*randn(1,1024);%添加高斯白噪声fx1=fft(x1);figure(2);subplot(211);plot(x1)subplot(212);plot(abs(fx1))这里是我写的程序：t3=-20:0.001:20;x3=sinc(t3/pi);x31=fft(x3);figure(1);subplot(2,2,1),plot(t3,x3),title(Sampling signal);grid on;subplot(2,2,2),plot(t3,abs(x31)),title(Sampling signal abs);grid on;subplot(2,2,3),plot(t3,angle(x31)),title(Sampling signal angle);grid on;我来帮他解答 图片符号编号排版地图满意回答问题在这一句：subplot(2,2,2