通信原理 课程设计 全Word格式.docx
- 文档编号:8012752
- 上传时间:2023-05-09
- 格式:DOCX
- 页数:18
- 大小:93.26KB
通信原理 课程设计 全Word格式.docx
《通信原理 课程设计 全Word格式.docx》由会员分享,可在线阅读,更多相关《通信原理 课程设计 全Word格式.docx(18页珍藏版)》请在冰点文库上搜索。
plot(t,u(1:
调制信号u(t)'
[M,m,df1]=fftseq(m,ts,df);
%
(2)m(t)和u(t)的频谱图
M=M/fs;
[U,u,df1]=fftseq(u,ts,df);
U=U/fs;
[C,c,df1]=fftseq(c,ts,df);
f=[0:
df1:
df1*(length(m)-1)]-fs/2;
figure
(2)
plot(f,abs(fftshift(M)))
frequency'
基带信号频谱M(f)'
plot(f,abs(fftshift(U)))
调制信号频谱U(f)'
signal_power=sum(u.*u)*dt/(length(t).*dt);
%(3)周期T0=t0,求已调信号的功率
signal_power
noise_power=signal_power/snr_lin;
%(4)噪声功率
noise_power
noise_std=sqrt(noise_power);
n=noise_std*randn(1,length(u));
r=u+n;
figure(5)%(5)噪声n(t)和已调信号加噪声(u(t)+n(t))的波形;
subplot(2,1,1)%噪声时域显示
plot(t,n(1:
噪声信号'
subplot(2,1,2)%加噪调制信号时域显示
plot(t,r(1:
信号加噪声信号'
[N,n,df1]=fftseq(n,ts,df);
%对噪声进行fft
N=N/fs;
%频率采样
[R,r,df1]=fftseq(r,ts,df);
%对最后的信号进行fft
R=R/fs;
figure(6)%(6)噪声n(t)和已调信号加噪声(u(t)+n(t))的频谱图;
subplot(2,1,1)%噪声频域显示
plot(f,abs(fftshift(N)))
axis([-20020000.15])
噪声信号频谱'
axis([00.15-33])
subplot(2,1,2)%加噪调制信号频域显示
plot(f,abs(fftshift(R)))
信号加噪声信号频谱'
axis([-50050000.15])
M文件
function[M,m,df]=fftseq(m,ts,df)
n1=fs/df;
n2=length(m);
n=2^(max(nextpow2(n1),nextpow2(n2)));
M=fft(m,n);
m=[m,zeros(1,n-n2)];
df=fs/n;
第二题
b=3000;
sn0=-20;
sn1=30;
ns=0.01;
N=(sn0:
ns:
sn1);
snrlin=10.^(N/10);
ct=b.*log2(1+snrlin.*(1/b));
plot(snrlin,ct)
s/n0'
高斯信道容量'
n1=25;
sn1=10.^(n1/10);
bs=10;
B=(3000:
bs:
300000);
b=1./B;
ct=B.*(log2(1+sn1.*b));
plot(B,ct)
B'
B无限增加时信道容量'
ct=(1/log
(2)*sn1)
ct=456.2202
第三题
clearall;
clf;
Tb=1;
f1=1000/Tb;
f2=f1+1/Tb;
phi=pi/4;
N=12000;
n=N/4;
t=0:
Tb/(n-1):
Tb;
%t=0:
1/(n-1):
1;
T=0:
4*Tb/(N-1):
4*Tb;
%T=0:
4*1/(N-1):
4*1;
s1=[cos(2*pi*f1*t+phi)cos(2*pi*f1*t+phi)cos(2*pi*f1*t+phi)cos(2*pi*f1*t+phi)];
s2=[cos(2*pi*f2*t+phi)cos(2*pi*f2*t+phi)cos(2*pi*f2*t+phi)cos(2*pi*f2*t+phi)];
%assumethetransmitsignalis"
1010"
ie.u2u1andu2
a=[cos(2*pi*f2*t+phi)cos(2*pi*f1*t+phi)cos(2*pi*f2*t+phi)cos(2*pi*f1*t+phi)];
%aisthereceivedsignal
b=a.*s1;
c=a.*s2;
forj=1:
4
fori=1:
n;
d((j-1)*n+i)=sum(b((j-1)*n+1:
(j-1)*n+i));
end;
end;
%求积分
e((j-1)*n+i)=sum(c((j-1)*n+1:
4;
ifd(j*n)>
e(j*n);
f((j-1)*n+1:
j*n)=zeros(1,n);
%输出0
else
j*n)=ones(1,n);
%输出1
%比较器
figure
(1)
plot(T,a);
时间'
);
ylabel('
幅度'
figurea'
holdon
figure
(2)
subplot(2,1,1);
plot(T,b);
figureb'
plot(T,c);
figurec'
figure(3)
subplot(3,1,1);
plot(T,d);
figured'
subplot(3,1,2);
plot(T,e);
figuree'
subplot(3,1,3);
plot(T,f);
figuref'
第四题
n=15;
k=5;
g=[10100110111];
m=mMatrix(k);
fori=0:
2^k-1
nkm=[m(i+1,:
)zeros(1,n-k)];
[qr]=deconv(nkm,g);
r=mod(r,2);
len_r=length(r);
r=[zeros(1,n-len_r)r];
c(i+1,:
)=mod(nkm+r,2)
fori=1:
2^k
weight(i)=sum(c(i,:
));
minW=weight
(2)
fori=3:
2^k;
ifweight(i)<
minW;
minW=weight(i);
c=
000000000000000
000010100110111
000101001101110
000111101011001
001000111101011
001010011011100
001101110000101
001111010110010
010001111010110
010011011100001
010100110111000
010110010001111
011001000111101
011011100001010
011100001010011
011110101100100
100001010011011
100011110101100
100100011110101
100110111000010
101001101110000
101011001000111
101100100011110
101110000101001
110000101001101
110010001111010
110101100100011
110111000010100
111000010100110
111010110010001
111101011001000
111111111111111
minW=
7
m=[10000
01000
00100
00010
00001];
5;
nkm=[m(i,:
c(i,:
)=mod(nkm+r,2);
G=c
H=[(G(:
k+1:
n))'
eye(n-k,n-k)]
B=[110010100110110];
S=mod(B*H'
2)
G=
H=
110101000000000
011010100000000
111000010000000
011100001000000
001110000100000
110010000010000
101100000001000
010110000000100
111110000000010
101010000000001
S=
0101001100
functionm=mMatrix(k)
%
m(1,:
)=[00000];
m(2,:
)=[00001];
fori=2:
2^k-1;
t=i;
e=1;
d=2^e;
a=[00000];
while1;
ift>
d;
e=e+1;
break;
ift<
2^e;
e=e-1;
ift-2^e==0;
a(k-e)=1;
ift-2^e>
0;
t=t-2^e;
a(k-e)=0;
m(i+1,:
)=a;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 通信原理 课程设计 通信 原理