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('I do')%读取音乐片段，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('幅度');grid on;sound(y,fs); 实验二[y,fs,bit]=wavread('I do')y1=y(:,1);n1=length(y1);D=2;%设置抽样间隔 y2=y1(1:D:n1);%减抽样n2=length(y2);%减抽样后信号长度t2=(0:n2-1)/fs;%设置横坐标figure subplot(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('幅度');grid on;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('幅度');grid on;sound(y3,fs/D)实验三[y,f...
