matlab通信原理实验.pdf

一： 低通信号抽样定理 1： 仿真思路： t+2 -2≤t≤-1 信号 x(t)= 1 -1t≤1 -t+2 1t≤2 0 其他 2 2 信号 x(t)经傅里叶变换可得 X(f)=4sin c (2f)-sinc (f) 2:程序代码 %参数设置 ts=0.2; fs=1/ts; df=0.01; f1=[-2.5:0.001:2.5]; %根据采样定理进行抽样 x=[zeros(1,10),[0:0.2:1],ones(1,9),[1:-0.2:0],zeros(1,10)]; %对采样信号进行 FFT [X,x,df1]=fftseq(x,ts,df); X1=X/fs; f=[0:df1:df1*(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 函数库 中。 下面是图形截屏： 二：模拟信号幅度调制解调及噪声对其影响 1，：仿真思路 （1） 首先利用载波对信号进行调制，采用常规的AM 调制方式。V (t) （1+am (t)） n cosw t c 2 2 0.5 （2） 利用包络检波方式解调。v(t)的包络v(t)=(u (t)+u (t)) c s （3） 加入噪声。 2：程序代码 echo on t0=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*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=signal_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; %有噪声时的解调信号 pause sub