matlab上机实验4(2)

matlab上机实验报告

h2=hd.*w2';

th1=(0:32)/fs; th2=(0:32)/fs; M=21184;

fh1=fft(h1,M);%矩形窗频谱函数 w=2/M*[0:M-1];

fh2=fft(h2,M);%布莱克曼窗频谱函数 figure(2)

subplot(2,2,1);plot(th1,h1) title('矩形窗时域');

subplot(2,2,2);plot(w,abs(fh1)); title('矩形窗频域');

subplot(2,2,3);plot(th2,h2); title('布莱克曼窗时域');

subplot(2,2,4);plot(w,abs(fh2)); title('布莱克曼窗频域') %解调后信号

n3=length(y3) t2=0:n3-1;

y4=cos(pi*1/2*t2)

y5=y3.*y4';%调制后的信号函数

%滤波处理

y6=conv(h1,y5);%用矩形窗对调制后信号进行滤波 f6=fft(y6); n4=length(f6); ty6=(0:n4-1)/fs; w3=2/n4*[0:n4-1]; %sound(y6,fs); figure(3)

subplot(2,2,1);plot(ty6,y6); title('矩形窗滤波后音频')

subplot(2,2,2);plot(w3,abs(f6)); title('矩形窗滤波后频谱')

y7=conv(h2,y5);%用布莱克曼窗对调制后的信号进行滤波 f7=fft(y7); n5=length(f7); ty7=(0:n5-1)/fs; w4=2/n5*[0:n5-1];

%sound(y7,fs);

subplot(2,2,3);plot(ty7,y7); title('布莱克曼窗滤波后音频')

subplot(2,2,4);plot(w4,abs(f7)); title('布莱克曼窗滤波后频谱')