06附录三 参考程序文档格式.docx
《06附录三 参考程序文档格式.docx》由会员分享，可在线阅读，更多相关《06附录三 参考程序文档格式.docx（36页珍藏版）》请在冰点文库上搜索。
plot（x,y,'
r-'
x,z,'
b:
'
x,i,'
）
gridon;
title（'
幅度增长的信号和信号的振幅y=exp（0.3*x）*cos（2*x+0.2）'
）;
结果为：
画出一般复指数信号（A）y=exp（-0.3*x）*cos（2*x+0.2）的模拟信号
z=exp（-0.3*x）;
i=-exp（-0.3*x）;
幅度减小的信号和信号的振幅y=exp（-0.3*x）*cos（2*x+0.2）'
结果为
第19页图1.26
画出一般复指数信号（A）y=exp（0.3*x）*cos（2*x+0.2）的离散信号
f1=50;
f2=120;
fs=5;
t=-10:
1/fs:
10;
n=t*fs;
y=（exp（-0.2*t））.*（cos（1.6*t+0.2））;
stem（n,y）;
幅度减小的离散信号和信号的振幅y=exp（-0.2*x）*cos（1.6*x+0.2）'
画出一般复指数信号（A）y=exp（-0.3*x）*cos（2*x+0.2）的离散信号
y=（exp（0.2*t））.*（cos（1.6*t+0.2））;
幅度增长的离散信号和信号的振幅y=exp（0.2*x）*cos（1.6*x+0.2）'
第44页图1.21
画出信号图形
程序为
ymst
（t+1）*（heaviside（t+2）-heaviside（t+1））+（heaviside（t+1）+heaviside（t）-2*heaviside（t-1））+（2-t）*（heaviside（t-1）-heaviside（t-2））'
ezplot（f,[-3,3]）%书上第1.12题要求画出信号若用分段函数较之用这种方法复杂而且这种方法可简单求出信号的各种变换
第一章完
第63页，图2.8
求俩信号的卷积：
程序：
n=-3:
8;
n1=-6:
16;
x=[0,0,0,1,1,1,1,1,0,0,0,0];
y=[0,0,0,1,2,4,8,16,32,64,0,0];
c=conv（x,y）%对x，y进行卷积
subplot（3,1,1）;
stem（n,x）;
subplot（3,1,2）;
subplot（3,1,3）;
stem（n1,c）;
结果：
c=
Columns1through12
000000137153162
Columns13through23
1241201129664000000
第74页，图2.22
求俩连续信号的卷积
p=0.01;
k1=0:
p:
1;
k2=-3:
-1;
f1=1+0*k1;
f2=-k2-1;
f=conv（f1,f2）;
%对两函数进行卷积
f=f*p;
k3=length（f1）+length（f2）-2;
%设定卷积后的函数的长度
k=0:
k3*p;
plot（k1,f1）;
plot（k2,f2）;
plot（k,f）;
%划出卷积后的函数
第137页，图3.7和第142页，图3.9
求周期方波的付里叶级数及方波付里叶级数表示的收敛吉伯斯现象
k=input（'
inputk'
）%输入k值以确定占空比
N=input（'
inputn'
）%输入n值以确定付利叶极数
0.01:
y=heaviside（x+4+2/k）-heaviside（x+4-2/k）+heaviside（x+2+2/k）-heaviside（x+2-2/k）+heaviside（x+2/k）-heaviside（x-2/k）+heaviside（x-2+2/k）-heaviside（x-2-2/k）+heaviside（x-4+2/k）-heaviside（x-4-2/k）;
subplot（3,1,1）,plot（x,y）;
%画出方波图像
n=-（N-1）/2:
1:
（N-1）/2;
a=sin（n*2*pi/k）./（n*pi）;
a（（N+1）/2）=2/k;
%计算出付利叶极数
subplot（3,1,2）,stem（n,a）;
方波付里叶级数'
%画出付利叶级数
A=1;
symst;
forA=1:
N
ifA==1
f=a（A）*cos（（A-（N-1）/2）*pi*t）;
else
f=f+a（A）*cos（（A-（N-1）/2）*pi*t）;
end
A=A+1;
%循环以计算付利叶级数表示的收敛
end
ezplot（f,[-1,1]）;
方波付里叶级数表示的收敛吉伯斯现象'
%图像有点误差主要误差出现在分步运算要误差小必须K值较大n值较大
执行结果
1．inputk8
k=
8
inputn55
N=
55
2．inputk12
12
inputn91
91
以上图所示可知要视图效果好必须k&
n值取得较得当。
第175页，图3.34
一阶递归离散时间滤波器的频率响应。
a=input（'
输入参数a'
%输入参数为LTI系统y[n]-ay[n-1]=x[n]中的参数
symsw;
f=abs（（1-a*cos（w）+j*a*sin（w））/（1-2*a*cos（w）+a^2））;
%LIT系统的频率响应的模
subplot（2,1,1）;
ezplot（f,[-2*pi,2*pi]）;
W=-2*pi:
2*pi;
x=1-1.2.*cos（W）+a^2;
t=angle（0.6.*sin（W）./x,（1-0.6.*cos（W））./x）;
%LIT系统的频率响应的弧度
subplot（2,1,2）;
plot（W,t）;
1．A=0.6
2．A=-0.6
第177页，图3.35
3点移动平均低通虑波器频率相应的模特性
w=-3*pi:
0.05:
3*pi;
f=abs（（1+2*cos（w））/3）;
%三点移动平均低通虑波器频率响应的模特性
plot（w,f）;
第203页，图4.2
连续时间方波信号
）%输入k以确定方波的占空比
subplot（2,1,1）,plot（x,y）;
y=-4*pi:
4*pi;
ify==0
a=2/k;
%y为0时分母为0必须用特殊方法计算出付利叶级数
else
a=sin（y*2*pi/k）./（y*pi）;
subplot（2,1,2）,plot（y,a）;
连续方波付里叶级数'
get（gca,'
XGrid'
'
on'
结果
第209页，图4.9
付利叶变换
W=input（'
输入参数'
x=-2*W:
2*W;
y=heaviside（x+W）-heaviside（x-W）;
t=-10*pi/W:
10*pi/W;
ift==0
Fjs=W/pi;
Fjs=（sin（W.*t））./（pi.*t）;
subplot（2,1,2）,plot（t,Fjs）;
第269页图5.14
时域和频域之间的相反的关系
x1=[0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0];
x2=[0,0,0,0,1,0,1,0,1,0,1,0,1,0,0,0,0];
x3=[0,0,1,0,0,1,0,0,1,0,0,1,0,0,1,0,0];
n=-8:
y1=sym（'
1+2*cos（2*w）+2*cos（w）'
%信号x1的频率响应
y2=sym（'
1+2*cos（4*w）+2*cos（2*w）'
%信号x2的频率响应
y3=sym（'
1+2*cos（6*w）+2*cos（3*w）'
%信号x3的频率响应
subplot（3,2,1）,stem（n,x1）;
subplot（3,2,2）,ezplot（y1,[-2*pi,2*pi]）;
subplot（3,2,3）,stem（n,x2）;
subplot（3,2,4）,ezplot（y2,[-2*pi,2*pi]）;
subplot（3,2,5）,stem（n,x3）;
subplot（3,2,6）,ezplot（y3,[-2*pi,2*pi]）;
第303页图6.1
相位变换的作用
1+cos（2*pi*t）*1/2+cos（4*pi*t）+cos（6*pi*t）*2/3'
subplot（4,1,1）,ezplot（y,[-pi,pi]）;
1+cos（2*pi*t+4）*1/2+cos（4*pi*t+8）+cos（6*pi*t+12）*2/3'
subplot（4,1,2）,ezplot（y2,[-pi,pi]）;
1+cos（2*pi*t+6）*1/2+cos（4*pi*t-2.7）+cos（6*pi*t+0.93）*2/3'
subplot（4,1,3）,ezplot（y3,[-pi,pi]）;
y4=sym（'
1+cos（2*pi*t+1.2）*1/2+cos（4*pi*t+4.1）+cos（6*pi*t-7.02）*2/3'
subplot（4,1,4）,ezplot（y4,[-pi,pi]）;
第322页图6.19图6.20
一阶系统单位冲击向应
e=input（'
输入一系数'
t=-1:
y1=exp（-t/e）.*heaviside（t）/e;
%一阶系统的单位冲击响应
y2=（1-exp（-t/e））.*heaviside（t）;
%一阶系统的阶跃响应
a=-2:
0.1:
3;
H=1./（j.*e.*e.*10.^a+1）;
bo=20.*log10（abs（H））;
%连续一阶系统的波特图
jiao=angle（H）;
subplot（4,1,3）,plot（a,bo）;
subplot（4,1,4）,plot（a,jiao）;
subplot（4,1,1）,plot（t,y1）;
subplot（4,1,2）,plot（t,y2）;
3266.226.23
不同阻尼下的响应
请输入阻尼系数'
%请输入阻尼系数
p=input（'
请输入自然频率'
%请输入自然频率
t=0:
u=heaviside（t）;
m=p/（2*sqrt（k^2-1））;
ifk>
0&
k<
1
h=p.*exp（-p*k.*t）.*sin（p.*sqrt（1-k^2）.*t）.*heaviside（t）/sqrt（1-k^2）;
%求出二阶系统的单位冲击响应
c1=-k*p+p*sqrt（k^2-1）;
c2=-k*p-p*sqrt（k^2-1）;
s=real（（1+m.*（exp（c1.*t）/c1-exp（c2.*t）/c2））.*heaviside（t））;
%求二阶系统的阶跃响应
elseifk>
h=p.*（exp（sqrt（k^2-1）.*p.*t）-exp（-sqrt（k^2-1）.*p.*t））.*heaviside（t）/2*sqrt（k^2-1）.*exp（-k.*p.*t）;
s=（1+m.*（exp（c1.*t）/c1-exp（c2.*t）/c2））.*heaviside（t）;
elseifk==1
h=（p^2）.*t.*exp（-p.*t）.*heaviside（t）;
s=（1-exp（-p.*t）-p.*t.*exp（-p.*t））.*heaviside（t）;
a=-1:
ifk==1
H=p^2./（j.*p.*（10.^a）+p）.^2
H=m./（j.*p.*（10.^a）-c1）-m./（j.*p.*（10.^a）-c2）;
bo=20*log10（abs（H））;
%求出二阶系统的波特图
jiao=angle（H）;
subplot（4,1,4）,plot（a,jiao）;
subplot（4,1,3）,plot（a,bo）;
subplot（4,1,2）,plot（t,s）
subplot（4,1,1）,plot（t,h/p）;
结果10<
k=0.3
2k>
k=6
k=1
画出绘制连续系统的零极点的函数在后有用
function[p,q]=sjdt（A,B）
%绘制连续系统零极点图程序
p=roots（A）;
%求出H（S）分母的根
q=roots（B）;
%求出H（S）分子的根
q=q'
;
p=p'
x=max（abs（[p,q]））;
%确定纵座标范围
x=x+0.1;
y=x;
%确定|横坐标范围
clf
holdon
axis（[-xx-yy]）;
axis（'
square'
plot（[-xx],[00]）;
%画横坐标
plot（[00],[-yy]）;
%画纵坐标
plot（real（p）,imag（q）,'
x'
plot（real（q）,imag（q）,'
o'
连续系统零极点图'
text（0.2,x-0.2,'
虚轴'
text（y-0.2,0.2,'
实轴'
第382页图7.16
对应于不同的W0值画出原始信号（兰线）和样本（小圆）和重建信号（绿线）
w=input（'
输入一参数'
x=linspace（0,4*pi,13）;
h=0:
a=cos（w*6*x）;
p=w*6;
symsu
b=cos（w*6*u）;
subplot（1,1,1）,ezplot（b,[0,4*pi]）;
holdon%画出原信号
t=interp1（x,a,h,'
spline'
%对采样点进行插值
plot_handles=plot（x,a,'
+'
h,t）;
%画出插值后的信号
stem_handles=stem（x,a）;
holdoff
subplot（1,1,1）,legend_handles=[plot_handles;
stem_handles
（1）];
%将三个图合在一张图上
legend（legend_handles,'
a'
a=cos（p*x）'
运行结果为
输入一参数W=1/6即W0=WS/6
w=
1.1667
可见没有发生混叠
W0=WS5/6时
发生了混叠
第397页图7.35
采样与抽取之间的关系（有问题）
第427页图8.13
幅度调制系统的输出
x=input（'
输入一函数必须包含tx（t）='
s'
g=eval（x,t）;
%将输入的字符转化为函数
y=（a+g）.*cos（50.*t）;
plot（t,y）;
a=0.5时
输入一函数必须包含tx（t）=cos（t-1）
输入一参数1
a=1时
输入一参数2
第475页图9。
2
已知信号求零极点
a=[1,3,2];
b=[1,-1];
sjdt（a,b）;
%利用函数SJDT画出零极点内容请看文件sjdt.m
结果为x为极点0为零点
第503页图10.4
%根据给出的两个极点和ROC可求出H（S）=1/t^2-2*x*t+x^2+y^2令x=2,y=3
a=[1,-4,13]
b=[5]
impulse（b,a）%划出冲击向应
第489页9.20
%已知二阶系统求极点与零点和系统的频率响应
q=input（'
请输入参数q〉0'
请输入参数w'
W=-3*pi:
ifq>
a1=-q*w+w*sqrt（q^2-1）;
a2=-q*w-w*sqrt（q^2-1）;
%求出系统的极点
subplot（2,2,1）,plot（a1,0,'
a2,0,'
elseifq>
q<
x=-q*w;
y1=w*sqrt（-q^2+1）;
y2=-w*sqrt（-q^2+1）;
subplot（2,2,1）,plot（x,y1,'
x,y2,'
else
fprintf（'
参数q输入错误'
H=w^2./（（j.*W）.^2+2*q*w.*W*j+w^2）;
h=abs（H）;
%系统的频率响应的模
Hj=angle（H）;
%系统的频率响应的弧度
subplot（2,2,2）,plot（W,h）;
subplot（2,2,3）,plot（W,Hj）;
q=0.1
w=1
q=0.2
q=0.4
第552页图10.13
w=-2*pi:
H=1./（1-a.*exp（-w.*j））;
%求出频率响应
h=abs（1./（1-a.*exp（-w.*j）））;
%求出频率响应的模
hdj=angle（H）;
%求出频率响应的复角
subplot（2,1,1）,plot（w,h）;
subplot（2,1,2）,plot（w,hdj）；
a=0.5
a=0.95
升级会员 