幅度调制AM信号输出和功率谱程序清单

1、Ml调制信号及其包络0051152253354456AM信号功率谱f%幅度调制AM信号输出和功率谱程序清单dt=0.001;%时间采样频谱fmax=1;%信源最高频谱fc=10;%载波中心频率T=5;%信号时长N=T/dt;t=0:N-1*dt;mt=sqrt(2)*cos(2*pi*fmax*t);%信源A=2;s_am=(A+mt).*cos(2*pi*fc*t);f,Xf=FFT_SHIFT(t,s_am);%调制信号频谱PSD=(abs(Xf).A2)/T;%调制信号功率谱密度figure(1)subplot(211);plot(t,s_am);holdon;%画出AM信号波形plot

2、(t,A+mt,'r-');%表示AM包络title('AM调制信号及其包络,);xlabel('t');subplot(212);%画出功率谱图形plot(f,PSD);axis(-2*fc2*fc01.5*max(PSD);title('AM信号功率谱,);xlabel('f);csm调制信号及其包络2LJL1II|IJ_00.511.522.633.5446比日信号功率借%抑制载波双边带调制DSBdt=0.001;%时间采样频谱fmax=1;%信源最高频谱fc=10;%载波中心频率T=5;%信号时长t=0:dt:T;mt=sqrt(

3、2)*cos(2*pi*fmax*t);%信源s_dsb=mt.*cos(2*pi*fc*t);f,sf=FFT_SHIFT(t,s_dsb);%调制信号频谱PSD=(abs(sf).A2)/T;%调制信号功率谱密度figure(1)subplot(211)plot(t,s_dsb);holdon;%画出DSB信号波形plot(t,mt,'r-');%标示mt波形title('DSB调制信号及其包络,);xlabel('t');subplot(212)plot(f,PSD);axis(-2*fc2*fc0max(PSD);title('DSB信号

4、功率谱,);%单边带调制SSBdt=0.001;%时间采样频谱fmax=1;%信源最高频谱fc=10;%载波中心频率T=5;t=0:dt:T;mt=sqrt(2)*cos(2*pi*fmax*t);%信源s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t);f,sf=FFT_SHIFT(t,s_ssb);%单边带信号频谱PSD=(abs(sf).A2)/T;%单便带信号功率谱figure(1)subplot(211)plot(t,s_ssb);holdon;%画出SSB信号波形plot(t,mt,'r-');%标示mt的包络title('SS

5、B调制信号');xlabel('t');subplot(212)plot(f,PSD);axis(-2*fc2*fc0max(PSD);title('SSB信号功率谱,);xlabel('f);%调用函数functionf,sf=FFT_SHIFT(t,st)%ThisfunctionisFFTtocalculateasignaltransfosrFourier%Input:t:samplingtime,st:signaldata.Timelengthmustgreaterthean2%output:f:samplingfrequency,sf:freq

6、uen%outputisthefrequencyandthesignalspectrumdt=t(2)-t(1);T=t(end);df=1/T;N=length(t);f=-N/2:N/2-1*df;sf=fft(st);sf=T/N*fftshift(sf);%I%显示模拟调制的波形及其解调方法VSB,文件名：VSB.m%Signaldt=0.001;fmax=5;fc=20;T=5;N=T/dt;t=0:N-1*dt;mt=sqrt(2)*(cos(2*pi*fmax*t)+sin(2*pi*0.5*fmax*t);%VSBmodulations_vsb=mt.*cos(2*pi*fc*

7、t);B1=0.2*fmax;B2=1.2*fmax;f,sf=FFT_SHIFT(t,s_vsb);t,s_vsb=vsbmd(f,sf,B1,B2,fc);%PowerSpectrumDensityf,sf=FFT_SHIFT(t,s_vsb);PSD=(abs(sf).A2)/T;%PlotVSBandPSDfigure(1)subplot(211)plot(t,s_vsb);holdon;plot(t,mt,'r-');title('VSB调制信号,);xlabel('t');subplot(212)plot(f,PSD);axis(-2*fc2

8、*fc0max(PSD);title('VSB信号功率谱,);xlabel('f);%functiont,st=vsbmd(f,sf,B1,B2,fc)%Thisfunctionisaresidualbandpassfilter%Inputsf:samplefrequency,sf:frequencyspectrumdata%B1:residualbandwidth,B2:highestfreqofthebasebandsignal%Outputst:sampletime,st:signaldatadf=f(2)-f(1);T=1/df;hf=zeros(1,length(f)

9、;bf1=floor(fc-B1)/df):floor(fc+B1)/df);bf2=floor(fc-B1)/df)+1:floor(fc+B2)/df);f1=bf1+floor(length(f)/2);f2=bf2+floor(length(f)/2);stepf=1/length(f1);hf(f1)=0:stepf:1-stepf;hf(f2)=1;f3=-bf1+floor(length(f)/2);f4=-bf2+floor(length(f)/2);hf(f3)=0:stepf:(1-stepf);hf(f4)=1;yf=hf.*sf;t,st=IFFT_SHIFT(f,yf);st=real(st);%