信号与系统MATLAB实验课后习题问题详解Word文档下载推荐.docx
《信号与系统MATLAB实验课后习题问题详解Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《信号与系统MATLAB实验课后习题问题详解Word文档下载推荐.docx(12页珍藏版)》请在冰豆网上搜索。
subplot(2,2,3);
impulse(sys,t);
title('
ImpulseResponse'
);
xlabel('
Time(sec)'
ylabel('
Amplitude'
subplot(2,2,4);
x_delta=zeros(size(t));
x_delta(t==0)=100;
[y1,t]=lsim(sys,x_delta,t);
y2=y1;
plot(t,y2);
2、
函数int1如下:
function[F,tF]=int1(f,tf,a)
T=tf
(2)-tf
(1);
F=zeros(size(tf));
tF=zeros(size(tf));
tF=tf;
forn=1:
length(tf)-1;
F(n+1)=F(n)+T*f(n);
end
验证如下:
t=[-1:
4];
e=zeros(size(t));
e=(t>
-1/2&
t<
1);
[z,zz]=intl(e,t,-1);
plot(zz,z);
第四章练习题
T1=1;
N1=10000;
t1=linspace(0,T1-T1/N1,N1)'
;
f1=1-2*t1;
OMG=32*pi;
K1=100;
omg=linspace(-OMG/2,OMG/2-OMG/K1,K1)'
X1=T1/N1*exp(-j*kron(omg,t1.'
))*f1;
fs1=OMG/2/pi/K1*exp(j*kron(t1,omg.'
))*X1;
T2=5;
N2=10000;
t2=linspace(0,T2-T2/N2,N2)'
fs2=0*t2;
f2=sawtooth(t2*2*pi,0);
X2=T2/N2*exp(-j*kron(omg,t2.'
))*f2;
fs2=fs2+OMG/2/pi/K1*exp(j*kron(t2,omg.'
))*X2;
plot(omg,abs(X1),'
r'
Frequency'
),ylabel('
)
单个锯齿周期幅频特性曲线'
plot(t1,fs1,'
Time'
Functionafterrecovered'
plot(omg,abs(X2),'
五个锯齿周期幅频特性曲线'
plot(t2,fs2,'
fsana函数如下:
functionF=fsana(t,f,N);
omg1=2*pi/(max(t)-min(t));
k=[-N:
N]'
F=1/length(t)*exp(-j*kron(k*omg1,t.'
))*f;
fssyn函数如下:
functionf=fssyn(F,t)
N=floor(length(F)/2);
N];
f=exp(j*kron(t,k*omg1))*F;
验证如下:
clc
clear
closeall
N1=256;
t=linspace(0,T1-T1/N1,N1)'
f=1-2*t;
subplot(3,1,1);
plot(t,f);
验证原函数'
N=25;
F1=fsana(t,f,N);
subplot(3,1,2);
stem(abs(F1),'
s'
前N项傅立叶级数系数幅度曲线'
f2=fssyn(F1,t);
subplot(3,1,3);
plot(t,f2);
time[s]'
傅立叶逆变换后时域函数'
第五章练习题
(a)Residue计算
a1=[1,5,6];
b1=[4,5];
[r1,p1,k1]=residue(b1,a1);
e1=zeros(size(t));
size(r1);
e1=e1+r1(n)*exp(p1(n)*t);
end;
subplot(1,2,1);
plot(t,e1);
residue计算'
t/s'
u/v'
lism仿真
sys1=tf([4,5],[1,5,6]);
delta=zeros(size(t));
delta(t==0)=100;
h1=lsim(sys1,delta,t);
subplot(1,2,2);
plot(t,h1);
lism仿真'
Residue计算和lism仿真结果一样
〔b〕Residue计算
t=[0,0.01,10];
e2=sin(t);
subplot(2,1,1);
plot(t,e2);
sys1=tf([1,0,2],[1,0,1]);
h2=lism(sys1,delta,t);
subplot(2,1,2);
plot(h,t2);
axis([0,10,-1,1]);
S=isstable(sys)函数:
Functions=isstable(sys);
X=ploe(sys);
S=1;
Forn=1:
Size(x)
Ifx(n)>
S=0;
break;
End;
稳定系统:
Sys=tf(1,[1,2]);
S=isstable(sys);
S=
1
不稳定系统:
Sys=tf(1,[1,-2]);
S=
第七章练习题
a=[1,0.5,-0.2,-0.1];
b=[1,-0.3];
n=[0:
10]'
[hi,t]=impz(b,a,n);
stem(n,hi);
u=(n>
=0);
hn=filter(b,a,u);
stem(n,hn);
n1=[0:
9]'
n2=[10:
19]'
x1=(n1>
x2=-(n2>
=10);
a1=[1,-0.2,-0.1];
a2=[1,-0.2,0.5];
b=[1,0.01];
[y1,wf1]=filter(b,a1,x1,[0,1]);
[y2,wf2]=filter(b,a2,x2,wf1);
stem(n1,y1);
holdon;
stem(n2,y2);