matlab调制程序.docVIP

clear; clear all; %μ÷??D?o? fs=100000; fm=1000; dt=1/fs; t=0:dt:0.01; mt=cos(2*pi*fm*t); subplot(2,1,1); plot(t,mt);grid on; title(μ÷??D?o?mt=cos(pi*t)μ?2¨D?); %??2¨D?o? fc=10000; Ct=cos(2*pi*fc*t); subplot(2,1,2); plot(t,Ct); title(??2¨D?o?2¨D?); figure; %AMμ÷?? A=1; SAM=(A+mt).*Ct; subplot(3,1,1); plot(t,SAM);grid on;hold on; plot(t,A+mt,r); axis([0 0.01 -2 2]); title(AMμ÷??oóμ?2¨D?); %DSBμ÷?? DSB=mt.*Ct; subplot(3,1,2); plot(t,DSB,b);grid on; axis([0 0.01 -2 2]); title(dsbμ÷??oóμ?2¨D?); subplot(3,1,3); plot(t,mt); axis([0 0.01 -2 2]); title(μ÷??D?o?); figure; %SSB?àò?·¨ SSB1=(1/2)*mt.*Ct+(1/2)*sin(2*pi*fm*t).*sin(2*pi*fc*t);%??±?′? SSB2=(1/2)*mt.*Ct-(1/2)*sin(2*pi*fm*t).*sin(2*pi*fc*t);%é?±?′? %???÷D?o??ó?μ?××?fft±??? N=2048; %2é?ù?μ?êoíêy?Yμ?êy n=0:N-1; t=n/fs; %ê±??DòáD y=fft(mt,N); %??D?o???DD?ì?ùFourier±??? mag=abs(y); y1=fft(Ct,N); mag1=(abs(y1)); y2=fft(DSB,N); mag2=(abs(y2)); y3=fft(SAM,N); mag3=(abs(y3)); y4=fft(SSB1,N); mag4=(abs(y4)); y5=fft(SSB2,N); mag5=(abs(y5)); f=n*fs/(1000*N); %?μ?êDòáD subplot(321); plot(f,mag); %??3????μ?ê±??ˉμ???·ù axis([0 20 0 400]); title(μ÷??D?o??μ?×); xlabel(?μ?ê/kHz); ylabel(??·ù);grid on; subplot(322); plot(f,mag1); axis([0 20 0 400]); title(??2¨D?o??μ?×); xlabel(?μ?ê/kHz); ylabel(??·ù);grid on; subplot(323); plot(f,mag2); axis([0 20 0 400]); title(DSBò?μ÷D?o??μ?×); xlabel(?μ?ê/kHz); ylabel(??·ù);grid on; subplot(324); plot(f,mag3); axis([0 20 0 400]); title(AMò?μ÷D?o??μ?×); xlabel(?μ?ê/kHz); ylabel(??·ù);grid on; subplot(325); plot(f,mag4); axis([0 20 0 400]); title(SSBò?μ÷??±?′?D?o??μ?×); xlabel(?μ?ê/kHz); ylabel(??·ù);grid on; subplot(326); plot(f,mag5); axis([0 20 0 400]); title(SSBé?±?′?ò?μ÷D?o??μ?×); xlabel(?μ?ê/kHz); ylabel(??·ù);grid on; figure; %????μíí¨??2¨?÷?μ?ê?ìó| wc=1.5*2*pi*fm/fs; N=1000; B=fir1(16,wc/pi); [h,w]=freqz(B, 1, N); x=w*fs/(2*pi)/1000; y=20*log10(abs(h)); B; plot(x,y);grid on;