内容发布更新时间 : 2024/11/15 3:00:15星期一 下面是文章的全部内容请认真阅读。

编程题实验一：

假设基带信号为m(t)?sin(2000?t)?2cos(1000?t)，载波频率为20kHz，仿真出AM、DSB-SC、SSB信号，观察已调信号的波形及频谱。（编程）

源文件：

clear all exec t2f.sci ; exec f2t.sci ;

N=2^16; //采样点数 fs=64; //采样频率 Bs=fs/2; //系统带宽 T=N/fs; //截短时间

t=-T/2+[0:N-1]/fs; //时域采样点 f=-Bs+[0:N-1]/T; //频域采样点

f0=1; phi=%pi/3; fc=10;//待观测正弦波的频率、幅度和初相 //以上是初始化参数设置

s0=sin((2*%pi)*f0*t)+2*cos((%pi)*f0*t); //原始基带信号 S0=t2f(s0,fs) ; //基带信号傅里叶变换 c=cos((2*%pi)*fc*t); //载波 //dsb-am调制

s1=s0.*c; //调制dsd-sc信号 S1=t2f(s1,fs) ; //dsb-sc傅里叶变换 //am调制

s2=0.8*s0.*c+c; //调制am信号 S2=t2f(s2,fs) ; //am调制信号傅里叶变换 //ssb调制

c1=sin((2*%pi)*fc*t); M=t2f(s0,fs);

MH=-%i*sign(f).*M; //在频域进行希尔伯特变换 sh=real(f2t(MH,fs)); //希尔伯特变换后的信号 s3=s0.*c-sh.*c1; //调制后ssb信号 S3=t2f(s3,fs); //ssb傅里叶变换

//以上是仿真计算部分，以下是绘图部分

xset ( \ // 原始信号波形 plot ( t ,s0) //

title( \原始信号波形\ xlabel (\

ylabel ( \

mtlb_axis( [0,10, -5,+5] )

xset ( \ //原信号幅度频谱 plot ( f ,abs(S0))

title( \调制信号的频谱图\

xlabel (\ ylabel (\

xset ( \ // dsb-sd调制信号波形及频谱 plot ( t ,s1)

title( \调制信号波形\ xlabel (\

ylabel ( \ mtlb_axis( [0,4, -5,+5] ) xset ( \plot ( f ,abs(S1) )

title( \调制信号的频谱图\ xlabel (\ ylabel (\

xset ( \ // am调制信号波形及频谱 plot ( t ,s2)

title( \调制信号波形\ xlabel (\

ylabel ( \ mtlb_axis( [0,4, -5,+5] ) xset ( \plot ( f ,abs(S2) )

title( \调制信号的频谱图\ xlabel (\ ylabel (\

xset ( \ // ssb调制信号波形及频谱 plot ( t ,s3)

title( \调制信号波形\ xlabel (\ ylabel ( \ mtlb_axis( [0,4, -8,+8] ) xset ( \plot ( f ,abs(S3) )

title( \调制信号的频谱图\ xlabel (\ ylabel (\

实验结果： 原始信号波形：、

原始信号频谱：

AM调制信号波形：