通信原理编程实验完整代码告各位学弟学妹书：此实验对于大多数人来说甚是烦人，但是我还是要劝你：尽力去学吧，学的时候觉得没什么用，那是因为不了解其博大精深，如果你一不小心顶住压力认真学了（譬如me），日后涉及相关领域时，你就能够触类旁通，瞬间迸发出灵感，撞击出智慧的火花，这种感觉就好像搂着一群漂亮的姑娘行走在春风中一样，原来那些朦胧的理论都变成了清晰的线条，串在一起：你会有一种类似于愚公般的给子子孙孙签订卖身契投身于学术的激进想法，不过这种极左的思想毕竟要不得——除非你想女朋友和你分手，于激励学习便是极好的。实验一：//基带信号m(t)=sin(2000*pi*t)+2cos(1000*pi*t),fc=20khz,求AM,DSB-SC,SSBclearallexect2f.sci;execf2t.sci;fs=800;//采样速率T=200;//截短时间N=T*fs;//采样点数dt=1/fs;//时域采样间隔t=[-T/2:dt:T/2-dt];//时域采样点df=1/T;//频域采样间隔f=[-fs/2:df:fs/2-df];//频域采样点数fm1=1;//待观测正弦波频率，单位KHz，下同fm2=0.5;//待观测余弦波频率fc=20;//载波频率//以上为初始化参数设置m1=sin((2*%pi)*fm1*t);//待观测正弦波部分M1=t2f(m1,fs);//傅里叶变换MH1=-%i*sign(f).*M1;//希尔伯特变换mh1=real(f2t(MH1,fs));//希尔伯特反变换m2=2*cos((2*%pi)*fm2*t);//待观测余弦波部分M2=t2f(m2,fs);//傅里叶变换MH2=-%i*sign(f).*M2;//希尔伯特变换mh2=real(f2t(MH2,fs));//希尔伯特反变换s1=(1+(m1+m2)/abs(max(m1+m2))).*cos((2*%pi)*fc*t);//AM信号时域表达式S1=t2f(s1,fs);//AM信号频域表达式s2=(m1+m2).*cos((2*%pi)*fc*t);//DSB-SC信号时域表达式S2=t2f(s2,fs);//DSB-SC信号频域表达式s3=(m1+m2).*cos((2*%pi)*fc*t)-(mh1+mh2).*sin((2*%pi)*fc*t);//SSB信号时域表达式，以上边带为例S3=t2f(s3,fs);//SSB信号上边带频域表达式//以上是仿真计算部分//以下为绘图部分//AM信号xset('window',1)plot(f,abs(S1))title('AM信号频谱')xlabel('f')ylabel('S(f)')mtlb_axis([-25,25,0,max(abs(S1))]);xset('window',2)plot(t,s1)title('AM信号波形')xlabel('t')ylabel('s(t)')mtlb_axis([-3,3,-3,3]);//DSB-SC信号xset('window',3)plot(f,abs(S2))title('DSB-SC信号频谱')xlabel('f')ylabel('S(f)')mtlb_axis([-25,25,0,max(abs(S2))]);xset('window',4)plot(t,s2)title('DSB-SC信号波形')xlabel('t')ylabel('s(t)')mtlb_axis([-1,4,-3,3]);//SSB信号（以上边带为例）xset('window',5)plot(f,abs(S3))title('SSB信号频谱')xlabel('f')ylabel('S(f)')mtlb_axis([-25,25,0,max(abs(S3))])xset('window',6)plot(t,s3)title('SSB信号波形')xlabel('t')ylabel('s(t)')mtlb_axis([0,6,-3,3])实验二：//m(t)=sin(2000*pi*t)+2cos(1000*pi*t)+4sin(500*pi*t+pi/3),fc=40KHz,Kf=5kHz/V，仿真产生FM，观察频谱和波形//初始化参数设置clearall;exect2f.sci;execf2t.sci;xdel(winsid());//关闭所有图形窗口fs=4000;//采样速率T=16;//截短时间N=T*fs;//采样点数dt=1/fs;//时域采样间隔t=[-T/2:dt:T/2-dt];//时域采样点df=1/T;//频域采样间隔f=[-fs/2:df:fs/2-df];//频域采样点数fm1=1;fm2=0.5;fm3=0.25;fc=40;Kf=5;//频偏常数//仿真计算m1=sin((2*%pi)*fm1*t);phi1=(2*%pi)*Kf*cumsum(m1)*dt;//cumsum函数为求累积和相当于积分m2=2*cos((2*%pi)*fm2*t);phi2=(2*%pi)*Kf*cumsum(m2)*dt;m3=4*sin((2*%pi)*fm3*t+%pi/3);phi3=(2*%pi)*Kf*cumsum(m3)*dt;s=cos((2*%pi)*fc*t+phi1+phi2+phi3);//FM信号时域表达式S=t2f(s,fs);//FM信号频域表达式m=sin(2*%pi*fm1*t)+2*cos(2*%pi*fm2*t)+4*sin(2*%pi*fm3*t+%pi/3);//调制信号表达式M=t2f(m,fs);//绘图//FM已调信号xset('window',1)plot(f,abs(S))title('FM信号频谱')xlabel('f')ylabel('S(f)')mtlb_axis([-80,80,0,max(abs(S))]);xset('window',2)plot(t,s)title('FM信号波形')xlabel('t')ylabel('s(t)')mtlb_axis([0,3,-2,2]);//调制信号xset('window',3)plot(f,abs(M))title('m(t)信号频谱')xlabel('f')ylabel('S(f)')mtlb_axis([-5,5,0,max(abs(M))]);xset('window',4)plot(t,m)title('m(t)信号波...
