信号与系统课设.docx
- 文档编号:26797495
- 上传时间:2023-06-22
- 格式:DOCX
- 页数:12
- 大小:127.64KB
信号与系统课设.docx
《信号与系统课设.docx》由会员分享,可在线阅读,更多相关《信号与系统课设.docx(12页珍藏版)》请在冰豆网上搜索。
信号与系统课设
一、
1.正弦信号
A=input('inputA=');%给正弦信号的幅度A赋值
w=input('inputw=');%给正弦信号的频率w赋值
theta=input('inputtheta=');%给正弦信号的初始相位theta赋值
disp(['这个信号是周期信号'])
T=2*pi/w
t=0:
0.01:
3*T;%定义时间点
ft=A*sin(w*t+theta);%th计算函数值
plot(t,ft);%画图
title('正弦信号');%为图像加标题注释
gridon;%在图上画方格
2.复指数信号
j00=sqrt(-1);%定义复数j
a=input('inputa=');%复指数信号赋值
w=input('inputw=');
K=input('inputK=');
ifa==0
disp('这是一个周期信号')
T=2*pi/w
elseifa>0
disp('这不是一个周期信号')
else
disp('这不是一个周期信号')
end
end
t=-1.5*abs(a):
0.01:
1.5*abs(a)%定义时间点
ft=K*exp((a+j00*w)*t);
subplot(2,2,1);plot(t,real(ft));title('实部');%画图
subplot(2,2,2);plot(t,imag(ft));title('虚部');
subplot(2,2,3);plot(t,abs(ft));title('模');
subplot(2,2,4);plot(t,angle(ft));title('相角');
3.实指数信号
A=input('inputA=')%给指数信号中的A赋值
a=input('inputa=')%给指数信号中的a赋值
t=0:
0.01:
10;%定义时间点
ft=A*exp(a*t);%计算函数值
plot(t,ft);%以时间为横轴,函数值为纵轴画图
title('实指数信号');%为图像加标题注释
gridon;%在图上画方格
三、信号的频谱图
w=0:
0.025:
5;
b=input('请依次输入bn,bn-1...b0:
');
a=input('请依次输入am,am-1...a0:
');
H=freqs(b,a,w);
subplot(2,1,1);
plot(w,abs(H));
gridon;
xlabel('\omega(rad/s)');
ylabel('|H(j\omega)|');
title('H(j\omega)的幅频特性');
subplot(2,1,2);
plot(w,angle(H));
gridon;
xlabel('\omega(rad/s)');
ylabel('\phi(\omega)');
title('H(j\omega)的相频特性');
四、连续时间信号卷积
dt=0.01;t=-1:
0.01:
4;
f1=input('inputf1=');
f2=input('inputf2=');
f=conv(f1,f2);%卷积
n=length(f);
tt=(0:
n-1)*0.01-2;
subplot(3,1,1);
plot(t,f1);
title('f1(t)');
subplot(3,1,2);
plot(t,f2);
title('f2(t)');
subplot(3,1,3);
plot(tt,f);
title('f1(t)*f2(t)');
离散时间信号卷积
nx=-1:
5;
nh=-2:
10;
x=input('inputf1=');
h=input('inputf2=');
y=conv(x,h);
ny1=nx
(1)+nh
(1);
ny=ny1+(0:
(length(nx)+length(nh)-2));
subplot(311)
stem(nx,x,'fill'),gridon
axis([-41603])
xlabel('n'),title('x(n)')
subplot(312)
stem(nh,h,'fill'),gridon
xlabel('n'),title('h(n)')
axis([-41603])
subplot(313)
stem(ny,y,'fill'),gridon
xlabel('n'),title('y(n)=x(n)*h(n)')
axis([-41603])
六、采样信号
Ts=1;
dt=0.1;
t1=-4:
dt:
4;
ft=input('inputft=');
subplot(221)
plot(t1,ft),gridon
axis([-44-0.11.1])
xlabel('Time(sec)'),ylabel('f(t)')
title('抽样前信号:
')
N=500;
k=-N:
N;
w=pi*k/(N*dt);
Fw=dt*ft*exp(-j*t1'*w);
subplot(222)
plot(w,abs(Fw)),gridon
axis([-1010-0.21.1*pi])
xlabel('\omega'),ylabel('F(w)')
title('采样前信号频谱')
t2=-4:
Ts:
4;
fst=input('inputfst=');
subplot(223)
plot(t1,ft,':
'),holdon
stem(t2,fst),gridon
axis([-44-0.11.1])
xlabel('Time(sec)'),ylabel('fs(t)')
title('采样后的信号:
'),holdoff
Fsw=Ts*fst*exp(-j*t2'*w);
subplot(224)
plot(w,abs(Fsw)),gridon
axis([-1010-0.21.1*pi])
xlabel('\omega'),ylabel('Fs(w)')
title('采样信号频谱:
')
七、拉普拉斯变换
f=input('inputf=')
F=laplace(f)%Laplace变换
b=input('以z的幂的降序输入分子系数b:
');%以z的幂的降序输入分子系数
a=input('以z的幂的降序输入分母系数a:
');%以z的幂的降序输入分母系数
[z,p,k]=tf2zp(b,a);%计算零极点和增益
plot(real(z),imag(z),'o');%画零点
holdon
plot(real(p),imag(p),'x');%画极点
xlabel('Re-axis');
ylabel('Im-axis');
title('zero-poleplot');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 系统