【现代信号处理第一次作业】

目录

1.1-1

1.4-1

1.4-2

1.4-3

4. 总结

1.1-1

用 matlab画出【例 1.1-1】所示的连续信号，并画出取样频率分别为

1000Hz和 5000Hz的离散信号，并分析这两个取样频率对应的信号是否有失真。

解：取样时间间隔分别为

连续信号采用plot函数，离散信号采用stem函数。

代码

%%1.1-1

clc;clear;

t1=0:0.001:0.1;

t2=0:0.0002:0.01;

x1=cos(100*pi*t1);

x2=cos(2100*pi*t1);

x3=cos(100*pi*t2);

x4=cos(2100*pi*t2);

subplot(3,2,1);

plot(t1,x1); ylabel ('x1'); xlabel ('t1');

subplot(3,2,2);

plot(t1,x2); ylabel ('x2'); xlabel ('t1');

subplot(3,2,3);

stem(t1,x1); ylabel ('x1'); xlabel ('t1');

subplot(3,2,4);

stem(t1,x2); ylabel ('x2'); xlabel ('t1');

subplot(3,2,5);

stem(t2,x3); ylabel ('x3'); xlabel ('t2');

subplot(3,2,6);

stem(t2,x4); ylabel ('x4'); xlabel ('t2');       

图像显示如下

当取样频率为1000Hz的离散信号，x1的频率为50Hz ，x2的频率为1050Hz，x2失真，对比两组信号波形图像一致的。

当取样频率为5000Hz的离散信号，x1的频率为50Hz ，x2的频率为1050Hz，两组信号都没有失真，对比两组信号波形是一致的。

1.4-1

解;

1.4-2

解：

1. ,采用matlab的conv卷积函数

Matlab代码：

%% 1.4-2

clc;clear;

t=0:5;

x1=[1 1 0 0 0 0];

x2=[1 1 1 0 0 0];

x22=[1 2 4 0 0 0];

y1=conv(x1,x2);

y2=conv(x1,x22);

Y1=y1(1:6);

Y2=y2(1:6);

subplot(2,3,1);

stem(t,x1); ylabel ('x1'); xlabel ('t'); grid on;

subplot(2,3,2);

stem(t,x2); ylabel ('x2'); xlabel ('t'); grid on;

subplot(2,3,3);

stem(t,Y1); ylabel ('Y1'); xlabel ('t'); grid on;

subplot(2,3,4);

stem(t,x1); ylabel ('x1'); xlabel ('t'); grid on;

subplot(2,3,5);

stem(t,x22); ylabel ('x22'); xlabel ('t'); grid on;

subplot(2,3,6);

stem(t,Y2); ylabel ('Y2'); xlabel ('t'); grid on;

图像如下：

（2）采用卷积公式，根据两者的特点分别列出不为零的点

只有这四个点不为零，其他都为零。

只有这四个点不为零，其他都为零。

1.4-3

序列xn=cos0.1πn+cos0.7πn，用 matlab画出xn、cos0.1πn、cos0.7πn的图像。并根据卷积的定义用 matlab计算出系统(a)、（b）的响应yn，并与输入信号作比较。

解：

1. 响应

响应yn波形与x2相似，有滤除x1信号的效果

1. 响应

响应yn波形与x2相似，有滤除x1信号的效果。说明系统响应与系统的参数有关。

Matlab代码如下，切换对应的h：

clc;clear;

t=0:96;

x1=cos(0.1*pi*t);

x2=cos(0.7*pi*t);

x=x1+x2;

h=[1 -3 3 -1];

%h=[1 3 3 1];

N=size(t);

L=size(h);

x3=zeros(1, L(2)-1);

X=[x3,x,x3];

REV_H=h(end:-1:1);

for i=1:N(2)+L(2)-1

    Y(i)=REV_H*(X(i:i+L(2)-1)');

end

subplot(2,2,1);%2*2图形中的第一个

stem(t,x1); ylabel ('x1'); xlabel ('t');

subplot(2,2,2);

stem(t,x2); ylabel ('x2'); xlabel ('t');

subplot(2,2,3);

stem(t,x); ylabel ('x'); xlabel ('t');

subplot(2,2,4);

stem(Y); ylabel ('Y'); xlabel ('t');

4. 总结

        本文完成了现代信号处理第一次作业，后期会分享更多有趣物联网的操作从而实现对外部世界进行感知，充分认识这个有机与无机的环境，科学地合理地进行创作和发挥效益，然后为人类社会发展贡献一点微薄之力。