用matlab实现数字信号处理实现调幅音乐信号的同时解调.docVIP

clear all; close all;clc [x,fs,bit]=wavread(F:\费玉清-一剪梅 00_01_23-00_01_28.wav); x=x(:,1); % sound(x,fs); N=length(x); figure subplot(2,1,1); plot(x) F1=fft(x,N); grid on title(原始音乐信号时域波型) xlabel(Time) ylabel(Magnitude) w1=2/N*[0:N-1]; subplot(2,1,2); plot(w1,abs(F1)) grid on title(原始音乐信号频域波型) xlabel(Frequency/Hz) ylabel(Magnitude) n=0:N-1; x1=cos(n*0.95); figure subplot(2,2,1); plot(x1) title(余弦信号时域波型) xlabel(Time) ylabel(Magnitude) F2=fft(x1,N); subplot(2,2,2); plot(w1,abs(F2)) grid on title(余弦信号频域波型) xlabel(Frequency/Hz) ylabel(Magnitude) X=x.*x1; subplot(2,2,3); plot(X) F3=fft(X,N); grid on