南通大学通信原理课程设计4题全_南通大学题库

其他范文 时间:2020-02-28 04:31:45 收藏本文下载本文
【www.daodoc.com - 其他范文】

南通大学通信原理课程设计4题全由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“南通大学题库”。

第一题 clear clc clf t0=0.15;ts=0.001;fc=250;snr=20;fs=1/ts;df=0.3;t=(0:ts:t0);snr_lin=10^(snr/10);

%信噪比

m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];c=cos(2.*pi.*fc.*t);

%载波 u=[2+0.85*m].*c;[M,m,df1]=fftseq(m,ts,df);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;dt=0.01;signal_power=sum(u.*u)*dt/(length(t).*dt);noise_power=signal_power/snr_lin;noise_std=sqrt(noise_power);n=noise_std*randn(1,length(u));r=u+n;

[N,n,df1]=fftseq(n,ts,df);

%对噪声进行fft N = N/fs;

%频率采样 [R,r,df1]=fftseq(r,ts,df);

%对最后的信号进行fft R = R/fs;signal_power noise_power figure(1)plot(t,m(1:length(t)))xlabel('Time')title('基带信号')figure(2)subplot(2,1,1)plot(t,c(1:length(t)))xlabel('Time')title('载波信号')subplot(2,1,2)plot(t,u(1:length(t)))xlabel('Time')title('调制信号')figure(3)subplot(2,1,1)plot(f,abs(fftshift(M)))xlabel('frequency')title('基带信号频谱')subplot(2,1,2)plot(f,abs(fftshift(U)))title('调制信号频谱')xlabel('frequency')figure(4)

subplot(2,1,1)

plot(t,n(1:length(t)))title('噪声信号')subplot(2,1,2)

plot(f,abs(fftshift(N)))axis([-200 200 0 0.15])title('噪声信号频谱')figure(5)subplot(2,1,1)

plot(t,r(1:length(t)))title('信号加噪声信号')axis([0 0.15-3 3])subplot(2,1,2)

plot(f,abs(fftshift(R)))title('信号加噪声信号频谱')axis([-500 500 0 0.15])

fftseq文件:

function [M,m,df]=fftseq(m,ts,df)fs=1/ts;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;signal_power =

2.0484 noise_power =

0.0205

%噪声时域显示

%噪声频域显示

%加噪调制信号时域显示

%加噪调制信号频域显示

示范课题:

第二题 clear clc clf 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)xlabel('s/n0')title('高斯信道容量')

clear clc clf n1=2.5;sn1=10.^n1;b0=300000;bs=10;B=(3000:bs:b0);b=1./B;ct=B.*(log2(1+sn1.*b));plot(B,ct)xlabel('B')title('信道容量')ct=(1/log(2)*sn1)

第三题 clear all;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)];% aume the transmit signal is “1010”, ie.u2 u1 and u2 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)];% a is the received signal b=a.*s1;c=a.*s2;

for j=1:4

for i=1:n;

d((j-1)*n+i)=sum(b((j-1)*n+1:(j-1)*n+i));

end;end;%求积分 for j=1:4

for i=1:n;

e((j-1)*n+i)=sum(c((j-1)*n+1:(j-1)*n+i));

end;end;%求积分

for j=1:4;

if d(j*n)>e(j*n);

f((j-1)*n+1:j*n)=zeros(1,n);%输出0

else

f((j-1)*n+1:j*n)=ones(1,n);%输出1

end;end;%比较器 figure(1)plot(T,a);xlabel('时间');ylabel('幅度');title('figure a');hold on figure(2)subplot(2,1,1);plot(T,b);ylabel('幅度');title('figure b');subplot(2,1,2)plot(T,c);ylabel('幅度');xlabel('时间');ylabel('幅度');title('figure c');hold on figure(3)subplot(3,1,1);plot(T,d);ylabel('幅度');title('figure d');subplot(3,1,2);plot(T,e);ylabel('幅度');title('figure e');subplot(3,1,3);plot(T,f);xlabel('时间');ylabel('幅度');title('figure f');

第四题(1)、n=15;k=5;g=[1 0 1 0 0 1 1 0 1 1 1];m=mMatrix(k);for i=0:2^k-1

nkm=[m(i+1,:)zeros(1,n-k)];

[q r]=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)end;

for i=1:2^k

weight(i)=sum(c(i,:));end;minW=weight(2)for i=3:2^k;

if weight(i)

minW=weight(i);

else

end;end;

全部系统码字:

c =

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0 1

00

0 10

01

0

0 01

0 01

0 0

01

0

0 0

0 0

0

0 0

000

0 11

0 1

0

0 1

01

0 11

01

0

0 0

00

0

0 1

001

0 0

0 10

01

0

0 1

00

0

0 1

0 010

0 11

0 010

0 01

0 1

0 00

0 10

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

最小重量:minW =7(2)、clear all;n=15;k=5;g=[1 0 1 0 0 1 1 0 1 1 1];m=[1 0 0 0 0

0 1 0 0 0

0 0 1 0 0

0 0 0 1 0

0 0 0 0 1];for i=1:5;

nkm=[m(i,:)zeros(1,n-k)];

[q r]=deconv(nkm,g);

r=mod(r,2);

len_r=length(r);

r=[zeros(1,n-len_r)r];

c(i,:)=mod(nkm+r,2);end;G=c H=[(G(:,k+1:n))' eye(n-k,n-k)] B=[1 1 0 0 1 0 1 0 0 1 1 0 1 1 0];S=mod(B*H',2)

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0 001

0 10

0 1

0 00

0 1

生成矩阵: G =

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

校验矩阵:H =

0

0

0

0

校验子: S =

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0 0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0 0

0 0

0 0

0 0

0 0

0 0

0 0

0 1

0 0

mMatrix文件: function m=mMatrix(k)% m(1,:)=[0 0 0 0 0];m(2,:)=[0 0 0 0 1];for i=2:2^k-1;

t=i;

e=1;

d=2^e;

a=[0 0 0 0 0];

while 1;

if t>d;

e=e+1;

d=2^e;

else

break;

end;

end;

if t

e=e-1;

end;

while 1;

if t-2^e==0;

a(k-e)=1;

break;

end;

if t-2^e>0;

a(k-e)=1;

t=t-2^e;

e=e-1;

else

a(k-e)=0;

e=e-1;

end;

end;

m(i+1,:)=a;

a=[0 0 0 0 0];end;

通信原理 课程设计 全

第一题 clear echo off t0=0.15; ts=0.001; fc=250; snr=20; fs=1/ts; t=(0:ts:t0); dt=0.01; df=0.3; snr_lin=10^(snr/10); m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),......

微机原理课程设计(南通大学)

微机原理 课 程 设 计指 导 书 南通大学计算机科学与技术学院计算机系统教研室二零零九年五月 目 录1 课程设计的目的和任务..................................................

成都大学通信原理课程设计

成都学院(成都大学)课程设计报告数字通信系统的设计与实现摘要:数字通信系统可分为数字频带传输通信系统、数字基带传输通信系统、模拟信号数字化传输通信系统。本设计主要介绍......

通信原理课程设计

数字幅度调制的抗噪声性能摘要:多进制数字振幅调制又称多电平调制。这种方式在原理上是通断键控(OOK)方式的推广。近几年它成了十分引人注目的一种高效率的传输方式。所谓高效......

通信原理课程设计

课设一一、设计题目信号特性分析(如正弦波信号的波形与频谱)二、设计目的通信原理课程设计是《通信原理》理论课的辅助环节。着重体现通信原理教学知识的运用,培养学生主动研究......

下载南通大学通信原理课程设计4题全word格式文档
下载南通大学通信原理课程设计4题全.doc
将本文档下载到自己电脑,方便修改和收藏。
点此处下载文档

文档为doc格式

热门文章
点击下载本文