matlab上机实验报告

4 AM调制音乐信号的同步解调

clear all;cla;close all

[a,fs,bit]=wavread('I do片段'); y1=a( : ,1);%去单列数据进行分析 f1=fft(y1); n=length(f1); tn=(0:n-1)/fs; w=2/n*[0:n-1]; %sound(y1,fs); figure(1)

subplot(2,2,1);plot(tn,y1);%绘制音乐信号时域波形图 grid on;

title('音频'); xlabel('时间'); ylabel('幅度');

subplot(2,2,2);plot(w,abs(f1));%绘制音乐信号频谱图 grid on;

title('频谱'); xlabel('频率/pi'); ylabel('幅度');

t=[0:n-1];

y2=cos(pi*1/2*t);%载波函数 y3=y1.*y2';%信号调制

ty3=(0:length(y3)-1)/fs;

subplot(2,2,3);plot(ty3,y3);%绘制调制后信号波形图 grid on;

title('AM调制音频信号'); xlabel('时间'); ylabel('幅度'); f3=fft(y3); n2=length(f3);

w2=2/n2*[0:n2-1];

subplot(2,2,4);plot(w2,abs(f3));%绘制调制后信号频谱图 grid on;

title('AM调制频谱'); xlabel('频率/pi'); ylabel('幅度'); %矩形窗和布莱克曼窗

N=33;wc=0.3*pi;%基于经验的指标，其中N为理想低通滤波器阶数，wc为截止频率 hd=ideal(N,wc);%调用理想低通滤波器函数 w1=boxcar(N);%产生各种窗函数 w2=blackman(N);

h1=hd.*w1';%加窗设计各种FIR滤波器