0123456-1-0.500.51幅值s时间（）信号波形00.20.40.60.811.21.41.61.820200040006000信号频谱数字角频率幅度0123-1-0.500.51幅度s时间（）2:1减抽样信号波形00.511.5201000200030002:1减抽样信号频谱w数字角频率幅度00.511.5-1-0.500.51幅度s时间（）5:1减抽样信号波形00.511.52050010005:1减抽样信号频谱w数字角频率幅度实验一[y,fs,bit]=wavread('Ido')%读取音乐片段，fs是采样率size(y)%求矩阵的行数和列数y1=y(:,1);%对信号进行分列处理n1=length(y1);%取y的长度t1=(0:n1-1)/fs;%设置波形图横坐标figuresubplot(2,1,1);plot(t1,y1);%画出时域波形图ylabel('幅值');xlabel('时间（s）');title('信号波形');subplot(2,1,2);Y1=fft(y1);w1=2/n1*(0:n1-1);%设置角频率plot(w1,abs(Y1));%画频谱图title('信号频谱');xlabel('数字角频率');ylabel('幅度');gridon;sound(y,fs);实验二[y,fs,bit]=wavread('Ido')y1=y(:,1);n1=length(y1);D=2;%设置抽样间隔y2=y1(1:D:n1);%减抽样n2=length(y2);%减抽样后信号长度t2=(0:n2-1)/fs;%设置横坐标figuresubplot(2,2,1);plot(t2,y2);%绘制减抽样信号波形图ylabel('幅度');xlabel('时间（s）');title('2:1减抽样信号波形');Y2=fft(y2);%对y2进行n2点fft谱分析w2=2/n2*[0:n2-1];subplot(2,2,3);plot(w2,abs(Y2));%绘制减抽样信号频谱图title('2:1减抽样信号频谱');xlabel('数字角频率w');ylabel('幅度');gridon;0246-1-0.500.51(s)时间幅度调制后信号0246-1-0.500.51(s)时间幅度调制后信号00.511.520100020003000w数字角频率幅度()调制后信号的频谱高频率调制00.511.520100020003000w数字角频率幅度()调制后信号的频谱低频率调制sound(y2,fs/D);D=5;%设置抽样间隔y3=y1(1:D:n1);%减抽样n3=length(y3);%减抽样后信号长度t3=(0:n3-1)/fs;%设置横坐标subplot(2,2,2);plot(t3,y3);%绘制减抽样信号波形图ylabel('幅度');xlabel('时间（s）');title('5:1减抽样信号波形');Y3=fft(y3);%对y2进行n2点fft谱分析w3=2/n3*[0:n3-1];subplot(2,2,4);plot(w3,abs(Y3));%绘制减抽样信号频谱图title('5:1减抽样信号频谱');xlabel('数字角频率w');ylabel('幅度');gridon;sound(y3,fs/D)实验三[y,fs,bit]=wavread('Ido')y1=y(:,1);n1=length(y1);n3=0:n1-1;b1=cos(0.75*pi*n3);%设置调制信号b2=cos(0.5*pi*n3);c1=b1'.*y1;%对原信号进行调制c2=b2'.*y1;lc1=length(c1);t=(0:lc1-1)/fs;figure%用载波对信号进行调制，并对其做fft变换subplot(2,2,1)%获取频谱，从图中可以观察到，调制后的plot(t,c1);%信号频谱发生搬移xlabel('时间(s)');ylabel('幅度');title('调制后信号');subplot(2,2,2)%获取频谱，从图中可以观察到，调制后的plot(t,c2);%信号频谱发生搬移xlabel('时间(s)');ylabel('幅度');title('调制后信号');w1=2/lc1*[0:lc1-1];%设置角频率WC1=fft(c1);C2=fft(c2);0510-1-0.500.51音频时间幅度0120200040006000频谱/pi频率幅度0510-1-0.500.51AM调制音频信号时间幅度0120100020003000AM调制频谱/pi频率幅度024x105-1-0.500.51时间幅度AM解调音频信号0120100020003000AM解调频谱/pi频率幅度subplot(2,2,3)plot(w1,abs(C1));xlabel('数字角频率w');ylabel('幅度');title('调制后信号的频谱(高频率调制)');gridon;subplot(2,2,4)plot(w1,abs(C2));xlabel('数字角频率w');ylabel('幅度');title('调制后信号的频谱(低频率调制)');gridon;sound(c1,fs);实验四clearall;cla;closeall[a,fs,bit]=wavread('Ido');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,3,1);plot(tn,y1);gridon;title('音频');xlabel('时间');ylabel('幅度');subplot(2,3,4);plot(w,abs(f1));gridon;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,3,2);plot(ty3,y3);%绘制调制后信号波形图gridon;title('AM调制音频信号');xlabel('时间');ylabel('幅度');f3=fft(y3);n2=length(f3);00.511.522.533.500.10.20.30.40.50...