内容发布更新时间 : 2024/12/25 13:48:54星期一 下面是文章的全部内容请认真阅读。

c2 = sin(2*pi*F2*n/Fs); FSK = a.*c1 - (a-1).*c2; subplot(4,1,3) stem(FSK) %% 2PSK Fc2 = 50;

d1 = sin(2*pi*Fc2*n/Fs); d2 = sin(2*pi*Fc2*n/Fs+pi); PSK = a.*d1 - (a-1).*d2; subplot(4,1,4) stem(PSK)

（4）产生DTMF信号的程序（拓展要求）

N = 200;

n = 0:500*N; t = 0:0.01:5*N; Fs = 44100;

noise = 4*sin(2*pi*10000*n/Fs);

DTMF = sin(2*pi*941*n/Fs) + sin(2*pi*1336*n/Fs);

（5）高斯白噪声序列的产生程序（扩展要求）

n = 0:100;

x8 = 0.1^0.5*randn(1,101)+0.15; subplot(2,1,1) stem(x8) mean(x8) var(x8)

x2 = [zeros(1,50) ones(1,51)]; subplot(2,1,2) stem(x2+x8)

（6）正弦信号频谱分析的程序（扩展要求）

n = 0:1000; Fs = 100;

x3 = sin(pi*n/16); x4 = sin(n/2);

x5 = cos(pi*n/12)+cos(3*pi*n/8);

fs=100;N=128; n=0:N-1;t=n/fs;

x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); % y=fft(x,N); % mag=abs(y); % f=n*fs/N; %

subplot(1,2,1),plot(f,mag/(128/2)); xlabel('Frequency/Hz');

ylabel('Magnitude');title('N=128');grid on;

subplot(1,2,2),plot(f(1:N/2),mag(1:N/2)/(128/2)); xlabel('Frequency/Hz');

ylabel('Magnitude');title('N=128');grid on;

% figure(2)

% plot(f(1:N/256),mag(1:N/256)/(128/2));

% xlabel('Frequency/Hz');

% ylabel('Magnitude');title('N=128');grid on;结果：

（1）信号处理Toolbox中的要求掌握及了解的函数列表。 列表详见前文，对要求掌握的函数列表进行了整理。

（2）x1[n]、x2[n]、y1[n]、y2[n]、x3[n]、x4[n]、x5[n]、x6[n]、x7[n]的时域波形 在?20?n?20内，画出单位下列信号：

(a)单位采样序列x1[n]??[n]和单位阶跃序列x2[n]?u[n]的时域波形图。(b)

y1[n]?x1[n?5]、y2[n]?x2[n?8]的波形。

画出下列信号在0?n?100内的波形。

??n?x3[n]?sin???16??n?x4[n]?sin??

2????n??3?n?x5[n]?cos???cos???12??8?

x3[n]、x4[n]、x5[n]是周期信号。当采样频率为1Hz时，信号的基波周期

是1/32、1/4、1/24和3/16的最小公倍数，即18/48=3/8，单位是秒。

（3）信号的时移：x1[n]与y1[n]、x2[n]与y2[n]之间的关系。 y信号是x信号的时移或反折的结果。 （4）正弦序列x3[n]、x4[n]、x5[n]周期的判断 分析请见（2）。

（5）指数序列底数与序列变化的关系总结 图象如下：

（6）2ASK、2FSK、2PSK调制信号时域波形（拓展要求） 图象如下：

（7）数字“0” DTMF时域波形（拓展要求）

该程序执行时可以播放一段0的音频，老师助教可以尝试一下。

（8）高斯白噪声序列的时域波形；正弦序列加上高斯白噪声后的时域波形

（扩展要求）

平均数和方差指标可以查看m代码执行的输出：