信号实验.docx
- 文档编号:8804132
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:12
- 大小:154.32KB
信号实验.docx
《信号实验.docx》由会员分享,可在线阅读,更多相关《信号实验.docx(12页珍藏版)》请在冰豆网上搜索。
信号实验
1)相加
clearall;
t=0:
0.0001:
3;
b=3;
t0=1;u=stepfun(t,t0);
n=length(t);
fori=1:
n
u(i)=b*u(i)*(t(i)-t0);
end%产生一个斜坡信号
y=sin(2*pi*t);
f=y+u;
plot(t,f);
xlabel('时间(t)');ylabel('幅值f(t)');title('连续信号的相加');
2)相乘
clearall;
t=0:
0.0001:
5;
b=3;
t0=1;u=stepfun(t,t0);
n=length(t);
fori=1:
n
u(i)=b*u(i)*(t(i)-t0);
end
y=sin(2*pi*t);
f=y.*u;
plot(t,f);
xlabel('时间(t)');ylabel('幅值f(t)');title('连续信号的相乘');
3)移向
clearall;
t=0:
0.0001:
2;
y=sin(2*pi*(t));
y1=sin(2*pi*(t-0.2));
plot(t,y,'-',t,y1,'--');
ylabel('f(t)');xlabel('t');title('信号的移位');
4)翻转
clearall;
t=0:
0.02:
1;t1=-1:
0.02:
0;
g1=3*t;
g2=3*(-t1);
gridon;
plot(t,g1,'--',t1,g2);
xlabel('t');ylabel('g(t)');
title('信号的反折');
5)尺度变换
clearall;
t=0:
0.001:
1;
a=2;
y=sin(2*pi*t);
y1=sin(2*a*pi*t);
subplot(2,1,1)
plot(t,y);
ylabel('y(t)');xlabel('t');
title('尺度变换');
subplot(2,1,2)
plot(t,y1);
ylabel('y1(t)');xlabel('t');
6)倒相
clearall;
t=-1:
0.02:
1;
g1=3.*t.*t;
g2=-3.*t.*t;
gridon;
plot(t,g1,'-',t,g2,'--');
xlabel('t');ylabel('g(t)');title('反相');
7)微分
clearall;
t=-1:
0.02:
1;
g=t.*t;
d=diff(g);
subplot(2,1,1);
plot(t,g,'-');
xlabel('t');ylabel('g(t)');title('微分');
subplot(2,1,2);
plot(d,'--');xlabel('t');ylabel('d(t)');
8)积分
clearall;
t=-1:
0.2:
1;symst
g=t*t;
d=int(g);
subplot(2,1,1);
ezplot(g);
xlabel('t');ylabel('g(t)');title('积分')
subplot(2,1,2)
ezplot(d);xlabel('t');ylabel('d(t)');
9)综合
symst
f=sym('(t/2+1)*(heaviside(t+2)-heaviside(t-2))');
subplot(2,3,1);ezplot(f,[-3,3]);
y1=subs(f,t,t+2);subplot(2,3,2);ezplot(y1,[-5,1]);
y2=subs(f,t,t-2);subplot(2,3,3);ezplot(y2,[-1,5]);
y3=subs(f,t,-t);subplot(2,3,4);ezplot(y3,[-3,3]);
y4=subs(f,t,2*t);subplot(2,3,5);ezplot(y4,[-2,2]);
y5=-f;subplot(2,3,6);ezplot(y5,[-3,3]);
10)卷积
(1)函数卷积计算
%函数卷积计算
a=1000;
t1=-5:
1/a:
5;
f1=stepfun(t1,0);
f2=stepfun(t1,-1/a)-stepfun(t1,1/a);
subplot(2,3,1);
plot(t1,f1);axis([-5,5,0,1.2]);%xlabel('时间(t)');ylabel('f1(t)');title('单位阶跃函数');
subplot(2,3,2);plot(t1,f2);ylabel('f2(t)');title('单位阶跃函数');
y=conv(f1,f2);r=2*length(t1)-1;t=-10:
1/a:
10;
subplot(2,3,3);plot(t,y);axis([-5,5,0,1.2]);title('f1与f2的卷积');
ylabel('y(t)');
f11=conv(f1,f1);f22=conv(f2,f2);
subplot(2,3,4);plot(t,f11);title('f1与f1的卷积');
ylabel('f11(t)');axis([-5,5,0,5000]);
subplot(2,3,5);plot(t,f22);title('f2与f2的卷积');ylabel('f22(t)');
(2)连续函数的卷积计算1
%计算连续信号的卷积积分
%f:
函数的样值向量
%k:
对应的时间向量
%:
采样时间间隔
s=0.01;
k1=0:
s:
2;
k2=k1;
f1=3*k1;
f2=3*k2;
f=conv(f1,f2);
f=f*s;
k0=k1
(1)+k2
(1);
k3=length(f1)+length(f2)-2;
k=k0:
s:
k3*s;
subplot(3,1,1);
plot(k1,f1);
title('f1(t)');
subplot(3,1,2);
plot(k2,f2);
title('f2(t)');
subplot(3,1,3);
plot(k,f);
title('f(t)');
(3)连续函数的卷积计算2
clearall;
T=0.1;t=0:
T:
10;f=sin(t);
h=0.5*(exp(-t)+exp(-3*t));Lf=length(f);Lh=length(h);
fork=1:
Lf+Lh-1y(k)=0;
fori=max(1,k-(Lh-1)):
min(k,Lf)
y(k)=y(k)+f(i)*h(k-i+1)
end
yzsappr(k)=T*y(k)
end
subplot(3,1,1);
plot(t,f);title('f(t)');
subplot(3,1,2);
plot(t,h);title('h(t)');subplot(3,1,3);
plot(t,yzsappr(1:
length(t)));title('卷积近似计算结果');xlabel('时间');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 实验