MATLAB实验报告(5-9)(4)

wm=2; wc=wm; Ts=2.5; n=-100:100; nTs=n*Ts;

fs=((1+cos(nTs))/2).*(uCT(nTs+pi)-uCT(nTs-pi)); t=-4:0.1:4;

ft=fs*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t)))); t1=-4:0.1:4;

f1=((1+cos(t1))/2).*(uCT(t1+pi)-uCT(t1-pi)); subplot(311); plot(t1,f1,':'),hold on axis([-4,4,-0.1,1.1]) xlabel('nTs'),ylabel('f(nTs)');

title('抽样间隔Ts=2时的抽样信号f(nTs)') hold off subplot(312) plot(t,ft),grid on axis([-4,4,-0.1,1.9]) xlabel('t'),ylabel('f(t)');

title('由f(nTs)信号重建得到有失真的升余弦脉冲信号') error=abs(ft-f1);

subplot(313) plot(t,error),grid on xlabel('t'),ylabel('error(t)');

title('重建信号与原升余弦脉冲信号的绝对误差')

(4).离散时间信号在MATLAB中的表示：

离散时间信号是指在离散时刻才有定义的信号，简称离散信号或者序列。离散序列通常用x(n)表示，自变量必须是整数。离散时间信号的波形绘制在MATLAB中，一般用stem函数，stem函数的基本用法和plot函数一样，亚绘制的波形的每个样本上都有一个小圆圈，默认是空心的。

1. 单位取样序列：

Function y =impDT(n)

Y=(n==0); //当参数为0时候的冲激为1，否则为0；

调用此函数的时候，n必须为整数或者整数向量。 例题12-1：

n=-3:3; x=impDT(n);

stem(n,x,'fill'),xlabel('n'),grid on title('单位冲击序列') axis([-3,3,-0.5,1])

2. 单位阶跃序列：

Function y = Udt(n)

Y=n>=0; //当参数为非负时输出为1; 调用此函数的时候，n必须为整数或者整数向量。 例12-2：

n=-3:5;

x=uDT(n);

stem(n,x,'fill'),xlabel('n'),grid on title('单位阶跃序列') axis([-3,5,-0.5,1.5])

3. 矩形序列

例12-3：

n=-3:8;

x=uCT(n)-uCT(n-5);

stem(n,x,'fill'),xlabel('n'),grid on title('矩形序列') axis([-3,8,-1,2])

4. 单边指数序列 5. 正弦序列

(5). 离散时间信号的基本运算

对离散时间序列实型基本运算可得到新的序列，这些基本运算包括加减乘除、、移位和反折等。

例12-7：

a=0.8;N=8;n=-12:12; x=a.^n.*(uCT(n)-uCT(n-N)); n1=n;n2=n1-3;n3=n1+2;n4=-n1; subplot(411) stem(n1,x,'fill'),grid on