źֵһֱֵ(n+1)Tʱ̵ǰһ˲ʱѵ(n+1)Tʱ̵IJֵһֱֵ(n+2)Tʱ̣ƣӶһбһĽźšΪÿһĽźŵֵΪֵ༴һΪ㣬ʳΪױ
(8)
Spectrum ScopeƵǣҪƵźŵƵźŵĹʡƵʡʧȡ
(9)
źŷ벻ͬź
3.2 ȵƽ
3.2.1 AMƽ漰
ϤSimulinkģ鹦ܵĻϣAMģͼSimulinkӷʵͼ3.1ʾ
ͼ3.1 AMƽʵͼ
źm(t)Ϊ1VƵΪ2Hzزc(t)Ϊ1VƵΪ200Hzزc'(t)Ϊ1VƵΪ200Hzͨ˲ֹƵΪ10Hzزźŵ
ʱΪ1/2000ױΪ1/1000
õźƵͼͼ3.2ѵźƵͼͼ3.3ʾ
ͼ3.2 AMźƵͼ
ͼ3.3 AMѵźƵͼ
Աͼ3.2ͼ3.3ֵ֪ƷʽռƵϿѵźŵƵǻźŵƵ
ͼͼ3.4ʾ
ͼ3.4 AMƽͼ
ͼ3.4еһǵźŲΣмѵźŲΣѽźŲΡ벻ͬĻźŶԵƽ
ͼ3.4пԷ֣ѽźŵķΪźŵ1/2ڱźŵİǵźŵһֱʵֵַAMƺزвֹʺزϣûϢĴͣЧʽϵͣƺͽ̼
.ܽ
4.1AMƵĴ
t=-1:0.00001:1;%ʱ A0=10;%زź A1=5;%Ʋź A2=3;%ѵź f=6000;%زźƵ w0=f*pi;%Ƶ m=0.15;%ƶ Ft=2000;%Ƶ
fpts=[100 120];%߽ͨƵfp=100Hz,ֹƵfs=120Hz mag=[1 0];
dev=[0.01,0.05];%ͨ1%5%
[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,Ft); %kaiscrordƲÿƵFIR˲
b21=fir1(n21,wn21,Kaiser(n21+1,beta));%FIRL˲
[h,w]=freqz(b21,1);%HɢϵͳƵ 0~piΧNƵʵȷֵֵw ΧNƵʵȷֵ㡣
%**********************زź********************** Uc=A0.*cos(w0*t);%زź subplot(5,2,1);
plot(t,Uc);%زź title('زź');
axis([0,0.01,-15,15]);% T1=10*fft(Uc);%Ҷ任 subplot(5,2,2);
plot(abs(T1));%زźƵ title('زźŵƵ');
axis([5800,6200,0,10000000]);%
%**************ź************************** mes=A1*cos(0.001*w0*t);%ź subplot(5,2,3);
plot(t,mes);%ź title('ź');