通信原理 matlab pdf,matlab通信原理实验.pdf

matlab通信原理实验matlab通信原理实验.pdf

一： 低通信号抽样定理

1： 仿真思路：

t+2 -2≤t≤-1

信号 x(t)= 1 -1

-t+2 1

0 其他

信号 x(t)经傅里叶变换可得 X(f)=4sin c2(2f)-sinc2(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+amn(t))

coswct

2 2 0.5

(2) 利用包络检波方式解调。v(t)的包络 v(t)=(uc (t)+us (t))

(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

subplot(2,1,1)

plot(t,m(1:length(t)))

axis([0 0.15 -2.1 2.1])

xlabel('时间')

title('消息信号')

pause

clf

subplot(2,1,1)

plot(t,u(1:length(t)))

axis([0 0.15 -2.1 2.1])

xlabel('时间')

title('已调信号')

subplot(2,1,2)

plot