x5=cos(pi*n/12)+cos(3*pi*n/8); subplot(3,1,3) stem(n,x5,'fill')

xlabel('n'),ylabel('x5'),title('x5[n]的时域波形图'); n=0:30; x6=0.2*(0.8).^n;

x7=0.2*exp((-1/12+1j*pi/6)*n); x61=0.2*(1.2).^n; x62=0.2*(-0.8).^n; figure(1) subplot(2,4,1) stem(n,x6,'fill');

title('x6[n] 底数为0.8'); xlabel('n'); ylabel('x6[n]');

subplot(2,4,2) stem (n,x61,'fill');

title('x61[n] 底数为1.2'); xlabel('n'); ylabel('x61[n]');

subplot(2,4,3) stem (n,x62,'fill');

title('x62[n] 底数为-0.8'); xlabel('n'); ylabel('x62[n]');

subplot(2,4,5) stem (n,real(x7),'fill');

title('x7[n]实部');xlabel('n'); ylabel('实部');

subplot(2,4,6) stem (n,imag(x7),'fill');

title('x7[n]虚部'); xlabel('n'); ylabel('虚部')

subplot(2,4,7) stem (n,abs(x7),'fill');

title('x7[n]幅度');xlabel('n'); ylabel('幅度')

subplot(2,4,8) stem(n,angle(x7),'fill');

title('x7[n]相位'); xlabel('n'); ylabel('相位');

（3）产生2ASK、2FSK、2PSK调制信号的程序（拓展要求） （4）产生DTMF信号的程序（拓展要求） f1=941; f2=1336; t=0:0.01:5;

m=sin(2*pi*f1*t)+sin(2*pi*f2*t); figure(1); plot(t,m);

title('数字0的时域波形图');xlabel('时间t');

（5）高斯白噪声序列的产生程序（扩展要求） n=0:100; N=length(n);

x8=0.15+sqrt(0.1)*rand(1,N); figure(1); subplot(3,1,1); stem(n,x8)

title('白噪声n，均值为0.15，方差为0.1') xlabel('n'); ylabel('x8[n]'); N2=41;

n2=-20:20;

x2=[zeros(1,20) ones(1,21)]; subplot(3,1,2); stem(n2,x2,'fill')

xlabel('n'),ylabel('x2[n]'),title('x2[n]的时域波形图'); x2_noise(1:20)=x2(1:20);

x2_noise(21:41)=x2(21:41)+x8(1:21); subplot(3,1,3); stem(n2,x2_noise,'fill')

xlabel('n'),ylabel('x2[n]_noise'),title('x2[n]加噪后的时域波形图');

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

（1）信号处理Toolbox中的要求掌握及了解的函数列表。 函数 用法 x=sawtooth(t,width) 功能 width是0到1之间的标量。在0到2π×width区间内，x的值从-1线性变化到1；在2π×width~2π区间内，x的值又从1线性变化到-1。sawtooth(t,1)和sawtooth(t)是等价的。 sinc square y = sinc(x) y=square(t,duty) 产生sinc函数 产生一个周期为2 π,幅值为±1的周期性方波, duty表示占空比(duty cycle) fft Y = fft(X,n) 返回 n 点 DFT。如果未指定任何值，则 Y 的大小与 X 相同。 ifft X = ifft(Y) 使用快速傅里叶变换算法计算 Y 的逆离散傅里叶变换。X 与 Y 的大小相同。 residuez [r,p,k] = residuez(b,a) [b,a] = residuez(r,p,k) 找出两个多项式b（z）和a（z）的比值的部分分数展开式的残基，极点和直接项。 sawtoooth x = sawtooth(t) 矢量b和a是以z的递减幂指定离散时间系统b（z）/ a（z）的多项式系数。 tf2zp zp2tf grpdelay [Z,P,K]=TF2ZP(NUM,DEN) 传递函数转换为零极点增益 [num,den]=zp2df(z,p,k) [gd,w] = grpdelay(b,a) 零极点增益转换为传递函数 返回由输入向量b和a指定的离散时间滤波器的群时延响应gd。 输入向量是分子b和分母a的系数,H(z)中的多项式。 zplane zplane(b,a) 其中b和a是行向量，首先使用根找到零和由分子系数b表示的传递函数的极点和分母系数a。 unwrap Q = unwrap(P) 当 P 的连续元素之间的绝对跃变大于或等于 π 弧度的默认跃变容差时，Q = unwrap(P) 将通过在向量 P 中增加 ±2π 的倍数来更正弧度相位角。如果 P 是矩阵，unwrap 将按列运算。如果 P 是一个多维数组，则 unwrap 对第一个非单一维度进行运算。

（2）x1[n]、x2[n]、y1[n]、y2[n]、x3[n]、x4[n]、x5[n]、x6[n]、x7[n]的时域波形