2、length(x)-1)]-fs/2;%解析法计算信号幅度谱y=4*(sinc(2*f1)).^2-(sinc(f1)).^2;%绘图指令subplot(2,1,1)plot(f1,abs(y));xlabel('f')title('解析法求得的x(t)幅度谱')subplot(2,1,2)plot(f,fftshift(abs(X1)));xlabel('f')title('根据采样定理(数值法)求得的x(t)幅度谱')另：FFT算法程序matlab上没有，需要在网上找到，并保存到matlab函数库中。下面是图形截屏：二：模拟信号幅度调制解
3、调及噪声对其影响1，：仿真思路（1）首先利用载波对信号进行调制，采用常规的AM调制方式。V(t)=（1+amn(t)）coswct220.5（2）利用包络检波方式解调。v(t)的包络v(t)=(uc(t)+us(t))（3）加入噪声。2：程序代码echoont0=0.15;%信号m(t)ts=0.001;%采样间隔fc=250;%载频a=0.85;%调制指数fs=1/ts;%采样频率t=[0:ts:t0];%时间取值df=0.25;m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*
4、ts)+1)];%信号c=cos(2*pi*fc.*t);%载波m_n=m/max(abs(m));%单位化信号[M,m,df1]=fftseq(m,ts,df);%傅里叶变化f=[0:df1:df1*(length(m)-1)]-fs/2;%频率向量u=(1+a*m_n).*c;%已调信号[U,u,df1]=fftseq(u,ts,df);%傅里叶变化env=env_phas(u);%包络检波dem1=2*(env-1)/a;signal_power=spower(u(1:length(t)));%已调信号的功率noise_power=si
5、gnal_power/100;%噪声功率noise_std=sqrt(noise_power);noise=noise_std*randn(1,length(u));r=u+noise;%添加噪声[R,r,df1]=fftseq(r,ts,df);env_r=env_phas(r);%包络检波dem2=2*(env_r-1)/a;%有噪声时的解调信号pausesubplot(2,1,1)plot(t,m(1:length(t)))axis([00.15-2.12.1])xlabel('时间')title('消息信号')pauseclfsubp
6、lot(2,1,1)plot(t,u(1:length(t)))axis([00.15-2.12.1])xlabel('时间')title('已调信号')subplot(2,1,2)plot(t,env(1:length(t)))xlabel('时间')title('已调信号的包络')pauseclfsubplot(2,1,1)plot(t,dem1(1:length(t)))xlabel('时间')title('解调后的信号')subplot(2,1,2)plot(t,dem2(1:length(t)))xlabel('时间')title(
7、'有噪声的解调后的信号')pauseclfsubplot(2,1,1)length(m)plot(f,fftshift(abs(M)));axis([-400,400,0,150]);title('原始信号的幅度谱')subplot(2,1,2)plot(f,fftshift(abs(U)));axis([-400,400,0,100]);title('已调信号的幅度谱');仿真结果截图：三：单边带调幅系统的建模仿真1：程序代码：t0=0.15;%信号持续时间ts=0.001;%采样时间间隔Fc=250;%载波频率Fs=1/ts;%采样频率d
8、f=0.3;%频率分辨率t=[0:ts:t0];%时间矢量m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/
