给予MATLAB的DSB调制解调代码

n=1024;fs=n; %设取样频率fs=1024Hz
m=640*pi; %产生调制信号s(t)
i=0:1:n-1;
t=i/n;
t1=(t-1/8).^2;t3=(t-3/8).^2;t4=(t-4/8).^2;
t6=(t-6/8).^2;t7=(t-7/8).^2;
s=exp(-m*t1)+exp(-m*t3)+exp(-m*t4)+exp(-m*t6)+exp(-m*t7);
c=cos(2*pi*300*t); %产生载波信号载波频率fc=300Hz
x=s.*c; %正弦波幅度调制(DSB)
y=x.*c; %解调
wp=0.1*pi;ws=0.12*pi;Rp=1;As=15;%设计巴特沃思数字低通滤波器
[N,wn]=buttord(wp/pi,ws/pi,Rp,As);
[b,a]=butter(N,wn);
s1=filter(b,a,y); %滤波
s1=2*s1;
S=fft(s,n); %求上述各信号及滤波器的频率特性
C=fft(c,n);
X=fft(x,n);
Y=fft(y,n);
[H,w]=freqz(b,a,n,'whole');
n2=300;
f=(-n/2:1:n/2-1);
subplot(521),plot(s);
title(‘s(t)’) ;
subplot(523),plot(c);
subplot(525),plot(x);
subplot(527),plot(y); subplot(529),plot(f,abs(fftshift(H)));
subplot(522),plot(f,abs(fftshift(S)));
subplot(524),plot(f,abs(fftshift(C)));
subplot(526),plot(f,abs(fftshift(X)));