西安邮电大学
专业课程设计报告书
院系名称：电子工程学院学生姓名：李群学号05113096
专业名称：光信息科学与技术班级：光信1103
实习时间：2014年4月8日至2014年4月
18日
一、课程设计题目：
用matlab 仿真光束的传输特性。
二、任务和要求
1、用matlab 仿真光束通过光学元件的变换。
① 设透镜材料为k9玻璃，对1064nm 波长的折射率为1.5062，镜片中心厚度为3mm ，凸面曲
率半径，设为100mm ，初始光线距离透镜平面20mm 。用matlab 仿真近轴光线（至少10条）经过平凸透镜的焦距，与理论焦距值进行对比，得出误差大小。
② 已知透镜的结构参数为101=r ，0.11=n ，51=d ，5163.121=='
n n （K9玻璃），
502-=r ，0.12='
n ，物点A 距第一面顶点的距离为100，由A 点计算三条沿光轴夹角分别为10、20、
30的光线的成像。试用Matlab 对以上三条光线光路和近轴光线光路进行仿真，并得出实际光线的球差大小。
③ 设半径为1mm 的平面波经凸面曲率半径为25mm ，中心厚度3mm 的平凸透镜。用matlab 仿
真平面波在透镜几何焦平面上的聚焦光斑强度分布，计算光斑半径。并与理论光斑半径值进行对比，得出误差大小。（方法：采用波动理论，利用基尔霍夫—菲涅尔衍射积分公式。） 2、用MATLAB 仿真平行光束的衍射强度分布图样。（夫朗和费矩形孔衍射、夫朗和费圆孔衍射、夫朗和费单缝和多缝衍射。） 3、用MATLAB 仿真厄米—高斯光束在真空中的传输过程。（包括三维强度分布和平面的灰度图。） 4、（补充题）查找文献，掌握各类空心光束的表达式，采用费更斯-菲涅尔原理推导各类空心光束在真空中传输的光强表达式。用matlab 对不同传输距离处的光强进行仿真。
三、理论推导部分
第一大题
（1）十条近轴光线透过透镜时，理想情况下光线汇聚透镜的焦点上，焦点到像方主平面的距离为途径的焦距F ，但由于透镜的折射率和厚度会影响光在传输过程中所走的路径（即光程差Δ）。在用MATLAB 仿真以前先计算平行光线的传输路径。，R 为透镜凸面的曲率半径，h 为入射光线的高度，θ1为入射光线与出射面法线的夹角，θ2为出射光线与法线的夹角，n 为透镜材料的折射率。设透镜的中心厚度为d ，则入射光线经过透镜的实际厚度为：L=(R-d) 光线的入射角为：sinq1=h/R 折射角度满足：sinq2=nsinq1 而实际的光束偏折角度为：θ2-θ1。
由此可以看出，当平行光线照射透镜时，在凸面之前光线平行于光轴，在凸面之后发生了偏折，于光轴交汇一点，这一点成为焦点f ，折线的斜率为（-tan(θ2-θ1)）。
（2）根据题意可得，本题所讨论的是与光轴夹角不同的三条光线，经过透镜的两次反射后的成像问题。利用转面公式计算。
入射光线与光轴的交点A到球面定点O的距离L;
入射光线与光轴的夹角U;
像方对应的用L’,U’表示；
根据折射定律可得 ,可由入射角求得折射角I’
sin I’=n/n’sin I
L’=r+rsin I’/rsin U’
计算完第一面以后，其折射光线就是第二面的入射光线。
转面公式 U2=U1’ L2=L1’-d1
（3）为了计算焦面上光强分布和光斑的大小，必须采用波动理论，利用基尔霍夫-菲涅耳衍射积分公式进行计算。
；
其中，、分别是、与之间的夹角。
推论从点光源Q0发射的单色光波，其波扰的数值大小与传播距离成反比，在位置以方程
表达为。又在其发射出的球面波的波前任意位置，与同向，
夹角。设定比例常数，，则可得到菲涅耳衍射积分公式。
第二大题
夫朗和费衍射装置
P 的光场，可以看做时开孔处入社波面各点次波波源发出的球面次波在P 点产生光场的叠加。由夫朗和费近似下的基尔霍夫公式
（1）矩形衍射：衍射孔为矩形孔 a= kax/2f b=kby/2f
P （x ，y ）点光强I=Io （sin a/a ）^2*(sin b/b)^2 沿x 轴的光强分布I=Io （sin a/a ）^2 沿y 轴的光强分布I=Io(sin b/b)^2 极大值
(2)单缝衍射：如果矩形孔一个方向的尺寸比另一个方向大得多，如b >>a ，则该矩形孔的衍射就变成一个单缝衍射
相应 P 点的光强为 I=Io （sin a/a ）^2
①当a ＝0，对应于thea ＝0 的衍射位置是光强中央主极大值(亮条纹)；
a=m*pi ，对应于满足
的衍射角方向为光强极小值(暗条纹)。
（3）多缝衍射：多缝是指在一块不透光的屏上，刻有 N 条等间距、等宽度的通光狭缝。当每个单缝等宽时，各套衍射条纹在透镜焦平面上完全重叠，其总光强分布为它们的干涉叠加 相应于 P 点的光强度为
2
20
02m ax sin lim
I I A I u p
===→α
α
β
β
ααsin sin ~),(~
022
2
11/)(11?
==?
?
--+-E dy dx e C y x E b b a a f yy xx ik sin 1, 2, (51)
m
m a
λ
θ==±±2
2
0sin sin 2() (56)
sin 2
N I P I ?α?α?
?
????=?
? ?
??
????
多缝衍射现象包含有衍射和干涉双重效应， 个狭缠的衍射光强关系式中包含有两个因子：一个是单缝衍射因子(sina/a)2,另外一个因子是[sin(N /2)/sin(/2)]2
多缝衍射主极大强度为 它们是单缝衍射在各级主极大位
置上所产生强度的N2 倍。在两个主极大之间，有(N －1)个极小。
（4）圆孔衍射
设圆孔半径为a, 圆孔中心O1位于光轴上，则圆孔上任一点Q 的位置坐标为ρ11，与相
应的直角坐标x1, y1的关系为 P(ρ,)点的光强为
式中Φ=ka J1(Φ)- -一阶贝塞尔
(bessel)函数 P 点的相对光强分布
极大、极小位置，令
()0/d =Φ
d I I o 求出I/O I 的极值点 四、Matlab 仿真部分
第一题（1） clear all
r=100;%透镜的曲率半径为100mm ， n1=1.5;%透镜的折射率n1=1.5 n2=1;%空气的折射率n2=1 d=3; %中心厚度
c=20;%距离透镜20mm 出入射
2
2
M 0sin (59)I N I αα??= ???1
11111sin ,cos ?ρ?ρ==y x 2
202
10212
22)(,])(2[])(2[
)(),(f A S I J I ka ka J C a I λθθπ?ρ=Φ
Φ==()2
10
]
2[
Φ
Φ=J I I p
%x=77:0.1:320;
figure(1)
for n=-5:5
y1=0.5*n;
%hold on;
%plot(x1,y1);
a1=asin(y1/r);%入射角
a2=asin(n1/n2*(y1/r));%折射角
a=a2-a1;
k=tan(a);%出射光线的斜率
x1=sqrt(r^2-y1^2);
%圆X^2+Y^2=R^2与入射光线y=y1的相交点，即为出射光线经过的一点 x=(r-c-d):0.01:x1;%距离透镜20mm出入射
hold on;
plot(x,y1); %平行光束
x2=x1:0.01:320;
y=-k*(x2-x1)+y1;
%出射光线
hold on;
plot(x2,y);
end
第一题（2）
clear all;
n=1.5163;
d=5;
r1=10;
r2=-50;
L=-100;
x0=-100;
u=[-pi/180 -pi/90 -pi/60];%三个角度
i=asin((L-r1)*sin(u)/r1)%第一入射角
i1=asin(sin(i)/n)%第一折射角
u1=u+i-i1%第一折射光与光轴夹角
s=u+i%两个夹角之和
g=cos(s)*r1%计算误差时的距离
h=sin(s)*r1%第一高度
t=r1-g%小误差
u2=u1%砖面公式
L1=r1+r1*sin(i1)/sin(u1)%第一相距
L2=L1-d%第二物距
i2=asin((L2-r2)*sin(u2)/r2)%第二入射角
i3=asin(n*sin(i2))%第二折射角
u3=u2+i2-i3%最终像与光轴夹角
h2=h-d*u2;%第二高度
x=linspace(-100,100,1000);
figure;
for a=1:3;
i(a)=asin((L-r1)*sin(u(a))/r1)
i1(a)=asin(sin(i(a))/n)
u1(a)=u(a)+i(a)-i1(a)
s(a)=u(a)+i(a)%he jiao du
g(a)=cos(s(a))*r1%zong chang
t(a)=r1-g(a)%xiao chang
h(a)=sin(s(a))*r1%gao du
u2(a)=u1(a)
L1(a)=r1+r1*sin(i1(a))/sin(u1(a))
L2(a)=L1(a)-d
i2(a)=asin ((L2(a)-r2)*sin(u2(a))/r2) i3(a)=asin (sin(i2(a))*n)
u3(a)=u2(a)+i2(a)-i3(a)
h2(a)=h(a)-d*u2(a)
for b=1:1000;
if x(b)<=t(a)
y(b)=tan(-u(a))*(100+x(b));
else if x(b)>=5
y(b)=tan(-u3(a))*(x(b)-5)+h2(a);
else
y(b)=tan(-u2(a))*(x(b)-t(a))+h(a); end
end
end
plot(x,y)
hold on;
end
% axis([-20 20 -1 2])
title('llll');
xlabel('\Delta');
第一题（3）
clear all
n=1.5062;%K9玻璃的折射率
d=3;%透镜的中心厚度
R=25;%透镜凸面曲率半径
f=R/(n-1);%透镜焦距
R0=1;%入射光束半径
lambda=1.064e-3;%波长
k=2*pi/lambda;
phy=lambda*0.61/R0;%角半径
z=f;
rmax=3*f*phy;%艾利斑半径
r=linspace(0,rmax,100);%产生从0到rmax之间的100点行矢量将衍射半径100等分
eta=linspace(0,2*pi,100);%将0到2*pi100等分
[rho,theta]=meshgrid(r,eta);%生成绘制3D图形所需的网格数据
[x,y]=pol2cart(theta,rho);%衍射斑某点的坐标转换极坐标到直角坐标
r0=linspace(0,R0,100);%将入射光束半径100等分
datr0=r(2)-r(1);%dat r
eta0=linspace(0,2*pi,100);
[rho0,theta0]=meshgrid(r0,eta0);
[x0,y0]=pol2cart(theta0,rho0);
deta=datr0*2*pi/100;%dat theta入射在透镜以前的光程差
for dx=1:100%都是为了建立网格
for dy=1:100
Rrho=sqrt((x-x0(dx,dy)).^2+(y-y0(dx,dy)).^2+z^2);%r 入射平面上网点与衍射平面上的距离
Rtheta=z./Rrho;%cos theta
opd=exp(j*k*(n-1)*(sqrt(R^2-rho.^2)-(R-d)));%引入透镜后的引入的光程差
Ep=-j./(lambda*Rrho)/2*exp(Rrho*j*k).*(1+Rtheta).*deta.*opd;%菲涅耳任意一点光场幅振幅表达式，取得平面光波振幅为1
E2(dx,dy)=sum(Ep(:));%积分公式的求和表达
end
end
Ie=conj(E2).*E2;%光强表达式
figure(1);
surf(x,y,Ie);
figure(2)
plot(x(50,:),Ie(50,:));
第二题（1）矩形衍射
clear all;
% 矩形孔
lamda=500e-9;
a=1e-3;
b=1e-3;
f=1;
m=500;
ym=8000*lamda*f;
ys=linspace(-ym,ym,m);
xs=ys;
n=255;
for i=1:m
sinth1=xs(i)/sqrt(xs(i)^2+f^2);
sinth2=ys./sqrt(ys.^2+f^2);
angleA=pi*a*sinth1/lamda;
angleB=pi*b*sinth2./lamda;
B(:,i)=(sin(angleA).^2.*sin(angleB).^2.*5000./(angleA.^2.*angleB.^2)); end
subplot(1,2,1)
image(xs,ys,B)
colormap(gray(n))
subplot(1,2,2)
plot(B(m/2,:),ys,'k')
第二题（2）单缝衍射
clear;
f=0.5;%透镜焦距
d=0.04e-3; %缝长
lumda=500e-9; %透镜波长
I0=8;
b=0.08e-3;%缝宽
x=-0.2:0.0001:0.2; %自变量向量
n=-2:0.0001:2;% 自变量向量
arf=n*pi/180;% 衍射角
a=(pi*b*sin(arf))/lumda;% 衍射因子角度
I1=I0*(sin(a)./a).^2;% 单缝衍射光强度
a(a==0)=eps;%当arf为0.取值为无穷小
figure;% 创建图形窗口
c=256;%颜色大小
br=(I1/max(I1))*c; % 提高颜色亮度
subplot(2,1,1)
image(a,I1,br)% 画图
colormap(gray(c));% 形成图形
subplot(2,1,2)
plot(arf,I1, 'b');
多缝衍射
clear all;
lamda=500e-9; %波长
N=2; %缝数，可以随意更改变换
a=2e-4;D=5;d=5*a;
ym=2*lamda*D/a;
xs=ym;
n=1001;
ys=linspace(-ym,ym,n);
for i=1:n
sina=ys(i)/D;%衍射角
alpha=pi*a*sina/lamda;%衍射因子的角度
beta=pi*d*sina/lamda;
B(i,:)=(sin(alpha)./alpha).^2.*(sin(N*beta)./sin(beta)).^2; B1=B/max(B);
end
NC=256; %确定灰度的等级
Br=(B/max(B))*NC;
subplot(1,2,1)
image(xs,ys,Br );
colormap(gray(NC)); %色调处理subplot(1,2,2)
plot(B1,ys,'k');
圆孔衍射
clear ;
lamda=600e-9;
a=0.0005;
f=3;
m=300;
ym=4000*lamda*f;
ys=linspace(-ym,ym,m);
xs=ys;
n=100;
for i=1:m
r=xs(i)^2+ys.^2;
sinth=sqrt(r./(r+f^2));
x=2*pi*a*sinth./lamda;
hh=(2*BESSELJ(1,x)).^2./x.^2; b(:,i)=(hh)'.*5000;
B=b/max(b);
end
subplot(1,2,1)
image(xs,ys,b)
colormap(gray(n))
%figure
subplot(1,2,2)
plot(xs,B,'R-')
title('夫琅禾费圆孔衍射光的强度分布','Fontsize',20,'Fontname','黑体')% 标题xlabel('xs','Fontsize',20,'Fontname','黑体')% 标记横轴
ylabel('强度B','Fontsize',20,'Fontname','黑体')% 标纵轴
五、画出仿真图形
第一题（1）
第一题（2）
第一题（3）
第二题（2）单缝衍射
第二题（4）圆孔衍射
六、仿真过程中发生的问题及解决的方法
（1）注意在书写程序时，中，英文字符输入，尤其是分号以及冒号。
（2）在写程序之前，要先将题目的含义搞懂，了解在计算过程中所需要的公式以及定义，逐个分析。
（3）特别要注意此软件编写语言虽然与C语言相似但有区别在编写时要注意。例：函数的调用以及for语句。
（4）注意一旦出现错误先查书写。
七、分析和总结
通过本次课程设计实验，使我对MATLAB软件有了深入的了解，MATLAB是一个功能庞大的商业数学软件，除了矩阵运算、绘制函数、数据图像等常用功能外还有很多的功能,利用它可以高效快捷的处理很多数学等方面问题。在以后的学习中我们要善于使用类似于MATLAB的其他工具软件来帮助我们解决学习上的一些问题，同时本次课程设计不仅是我们巩固了过去的知识也同时接触到了新的实物。
这次研究的课题是对过去所学的光学知识的应用，但我在初次接触MATLAB软件时完全不知从何下手，很茫然，又加之对过去所学知识有些淡忘，同时也掌握的不是很到位，使得我们很困惑，迷茫，一致于实验开始前几天可以说几乎毫无进展，后来花费了很长时间去看旧的知识及MATLAB软件应用相关的文档，情况才得以好转。在以后的学习中，要尽量避免这样的情况，在课题下来之后，首先应但是“看”，而不是“想”。
其实在这次实习中，我感触最深的是：“书到用时方恨少”，学习知识不应该囫囵吞枣，马马虎虎，应当细嚼慢咽，融会贯通，否则直至用时才发现知识体系支离破碎，面对问题毫无头绪。同时应该多跟身边的人讨论讨论，集思广益，在得到新的思路的同时也可以认识到自身的不足。
八、参考文献
【1】秦襄培 MATLAB 图像处理及界面编程宝典电子工业出版社
【2】袁东详解MATLAB快速入门及应用电子工业出版社
【3】安连生应用光学北京理工大学出版社
【4】石顺祥物理光学与应用光学西安电子科技大学出版社
九、过程考核表和成绩鉴定表
西安邮电大学电子工程学院专业课程设计过程考核表
西安邮电大学电子工程学院专业课程设计成绩鉴定表
华东交通大学MATLAB程序设计报告书 课题名称：基于MATLAB的粒子群优化算法的实现 姓名： 学号：20160280800014 专业：控制科学与工程 2016年 11月 20日
基于MATLAB的粒子群优化算法的实现 一、课程选题目的 本次课程设计的课题为《基于MATLAB的粒子群优化算法的实现》，主要为学会运用MATLAB对实际算法编程，加深对粒子群优化算法的理解，并为今后熟练使用MA TLAB进行系统的分析仿真和设计奠定基础。数值计算分析可以帮助更深入地理解理论知识，并为将来使用MA TLAB进行各领域数值分析分析和实际应用打下基础。 此次课程主要是为了进一步熟悉对MATLAB软件的使用，以及学会利用MA TLAB对数值运算这种实际问题进行处理，将理论应用于实际，加深对它的理解。 二、粒子群优化算法原理 优化是科学研究、工程技术和经济管理等领域的重要研究工具。它所研究的问题是讨论在众多的方案中寻找最优方案。例如，工程设计中怎样选择设计参数，使设计方案既满足设计要求又能降低成本；资源分配中，怎样分配有限资源，使分配方案既能满足各方面的基本要求，又能获得好的经济效益。在人类活动的各个领域中，诸如此类，不胜枚举。优化这一技术，正是为这些问题的解决，提供理论基础和求解方法，它是一门应用广泛、实用性很强的科学。近十余年来，粒子群优化算法作为群体智能算法的一个重要分支得到了广泛深入的研究，在路径规划等许多领域都有应用。 2.1 粒子群优化算法的起源 粒子群优化（PSO）算法是由Kennedy和Eberhart于1995年用计算机模拟鸟群觅食这一简单的社会行为时，受到启发，简化之后而提出的。 设想这样一个场景：一群鸟随机的分布在一个区域中，在这个区域里只有一块食物。所有的鸟都不知道食物在哪里。但是他们知道当前的位置离食物还有多远。那么找到食物的最优策略是什么呢。最简单有效的方法就是追寻自己视野中目前离食物最近的鸟。如果把食物当作最优点，而把鸟离食物的距离当作函数的适应度，那么鸟寻觅食物的过程就可以当作一个函数寻优的过程。鱼群和鸟群的社会行为一直引起科学家的兴趣。他们以特殊的方式移动、同步，不会相互碰撞，整体行为看上去非常优美。生物学家CargiReynolds提出了一个非常有影响的鸟群聚集模型。在他的模拟模型boids中，每一个个体遵循：避免与邻域个体相冲撞、匹配邻域个体的速度、试图飞向感知到的鸟群中心这三条规则形成简单的非集中控制算法驱动鸟群的聚集，在一系列模拟实验中突现出了非常接近现实鸟群聚集行为的现象。该结果显示了在空中回旋的鸟组成轮廓清晰的群体，以及遇到障碍物时鸟群的分裂和再度汇合过程。由此受到启发，经过简化提出了粒子群优化算法。 2.2粒子群优化算法的原理 在粒子群优化算法中，每个优化问题的潜在解都是搜索空间中的一只鸟，称之为“粒子”。所有的粒子都有一个由被优化的函数决定的适应值，每个粒子还有一个速度决定他们飞翔的方向和距离。然后粒子们就追随当前的最优粒子在解空间中搜索。优化开始时先初始化为一群随机粒子（随机解）。然后通过迭代找到最优解。在每一次迭代中，粒子通过跟踪两个极值来更新自己。第一个极值就是整个种群目前找到的最优解。这个极值是全局极值。另外也可以不用整个种群而只是用其中一部分作为粒子的邻居，那么在所有邻居中的极值就是局部极值。第二个极值是粒子本身所找到的最优解，称为个体极值。这是因为粒子仅仅通过跟踪全局极值或者局部极值来更新位置，不可能总是获得较好的解。这样在优化过程中，粒子在追随全局极值或局部极值的同时追随个体极值则圆满的解决了这个问题。这就是粒子群优化
控制系统仿真课程设计 （2011级） 题目控制系统仿真课程设计学院自动化 专业自动化 班级 学号 学生姓名 指导教师王永忠/刘伟峰 完成日期2014年6月
控制系统仿真课程设计一 ———交流异步电机动态仿真 一 设计目的 1．了解交流异步电机的原理，组成及各主要单元部件的原理。 2. 设计交流异步电机动态结构系统； 3．掌握交流异步电机调速系统的调试步骤，方法及参数的整定。 二 设计及Matlab 仿真过程 异步电机工作在额定电压和额定频率下，仿真异步电机在空载启动和加载过程中的转速和电流变化过程。仿真电动机参数如下： 1.85, 2.658,0.2941,0.2898,0.2838s r s r m R R L H L H L H =Ω=Ω===， 20.1284Nm s ,2,380,50Hz p N N J n U V f =?===，此外，中间需要计算的参数如下： 21m s r L L L σ=-，r r r L T R =，22 2 s r r m t r R L R L R L +=，10N m TL =?。αβ坐标系状态方程： 其中，状态变量： 输入变量： 电磁转矩： 2p m p s r s L r d ()d n L n i i T t JL J βααωψψβ=--r m r r s r r d 1d L i t T T ααβαψψωψ=--+r m r r s r r d 1d L i t T T ββαβψψωψ=-++22s s r r m m m s r r s s 2r r r r d d i R L R L L L L i u t L T L L ααβαα σψωψ+=+-+22 s s r r m m m s r r s s 2 r r r r d d i R L R L L L L i u t L T L L ββαββ σψωψ+=--+[ ] T r r s s X i i αβαβωψψ=[ ] T s s L U u u T αβ=()p m e s s s s r n L T i i L βααβ ψψ=-
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);
2DPSK调制与解调系统的仿真 设计原理 (1) 2DPSK信号原理 1.1 2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。 图1.1 2DPSK信号 在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。定义为本码元初相与前一码元初相之差，假设： →数字信息“0”； →数字信息“1”。 则数字信息序列与2DPSK信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK信号相位：0
或 ： 1.2 2DPSK 信号的调制原理 一般来说，2DPSK 信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。 图1.2.1 模拟调制法 2DPSK 信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0” 时接相位0，当输入数字信息为“1”时接pi 。 图1.2.2 键控法调制原理图 1.3 2DPSK 信号的解调原理 2DPSK 信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。 码变换 相乘 载波 s(t) e o (t)
、 北京理工大学珠海学院 课程设计任务书 2010 ～2011 学年第 2学期 学生姓名：林泽佳专业班级：08自动化1班指导教师：钟秋海工作部门：信息学院一、课程设计题目 ： 《控制系统建模、分析、设计和仿真》 本课程设计共列出10个同等难度的设计题目，编号为：[0号题]、[1号题]、[2号题]、[3号题]、[4号题]、[5号题]、[6号题]、[7号题]、[8号题]、[9号题]。 学生必须选择与学号尾数相同的题目完成课程设计。例如，学号为8xxxxxxxxx2的学生必须选做[2号题]。 二、课程设计内容 （一）《控制系统建模、分析、设计和仿真》课题设计内容|
！ " [2 有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见（二）。 （二）《控制系统建模、分析、设计和仿真》课题设计要求及评分标准【共100分】 , 1、求被控对象传递函数G(s)的MATLAB描述。（2分） 2、求被控对象脉冲传递函数G(z)。（4分） 3、转换G(z)为零极点增益模型并按z-1形式排列。（2分） 4、确定误差脉冲传递函数Ge(z)形式,满足单位加速度信号输入时闭环稳态误差为零和实际 闭环系统稳定的要求。（6分） 5、确定闭环脉冲传递函数Gc(z)形式，满足控制器Dy(z)可实现、最少拍和实际闭环系统稳 定的要求。（8分）
6、根据4、5、列写方程组，求解Gc(z)和Ge(z)中的待定系数并最终求解Gc(z)和Ge(z) 。 （12分） 7、求针对单位加速度信号输入的最少拍有波纹控制器Dy(z)并说明Dy(z)的可实现性。 （3分） ！ 8、用程序仿真方法分析加速度信号输入时闭环系统动态性能和稳态性能。（7分） 9、用图形仿真方法(Simulink)分析单位加速度信号输入时闭环系统动态性能和稳态性能。 （8分） 10、确定误差脉冲传递函数Ge(z)形式,满足单位速度信号输入时闭环稳态误差为零和实际 闭环系统稳定的要求。（6分） 11、确定闭环脉冲传递函数Gc(z)形式，满足控制器Dw(z)可实现、无波纹、最少拍和实际 闭环系统稳定的要求。（8分） 12、根据10、11、列写方程组，求解Gc(z)和Ge(z)中的待定系数并最终求解Gc(z)和Ge(z) 。 （12分） 13、求针对单位速度信号输入的最少拍无波纹控制器Dw(z)并说明Dw(z)的可实现性。（3分） 14、用程序仿真方法分析单位速度信号输入时闭环系统动态性能和稳态性能。（7分） 15、用图形仿真方法(Simulink)分析单位速度信号输入时闭环系统动态性能和稳态性能。 & （8分） 16、根据8、9、14、15、的分析，说明有波纹和无波纹的差别和物理意义。（4分） 三、进度安排 6月13至6月14：下达课程设计任务书；复习控制理论和计算机仿真知识，收集资料、熟悉仿真工具；确定设计方案和步骤。 6月14至6月16：编程练习，程序设计；仿真调试，图形仿真参数整定；总结整理设计、 仿真结果，撰写课程设计说明书。 6月16至6月17：完成程序仿真调试和图形仿真调试；完成课程设计说明书；课程设计答 辩总结。 [ 四、基本要求
MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件
MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等； 3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析， 针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： ①初步了解Matlab、熟悉Matlab界面、进行简单操作； ②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； ③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等； ④使用文本编辑器编辑m文件，函数调用； ⑤能进行简单的信号处理Matlab编程； ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计内容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献（不少于5篇）； ⑦其它必要内容等。 时间安排：1.5周（分散进行） 参考文献： [1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
西安邮电大学 专业课程设计报告书 院系名称：电子工程学院学生姓名：李群学号05113096 专业名称：光信息科学与技术班级：光信1103 实习时间：2014年4月8日至2014年4月 18日
一、课程设计题目： 用matlab 仿真光束的传输特性。 二、任务和要求 1、用matlab 仿真光束通过光学元件的变换。 ① 设透镜材料为k9玻璃，对1064nm 波长的折射率为1.5062，镜片中心厚度为3mm ，凸面曲 率半径，设为100mm ，初始光线距离透镜平面20mm 。用matlab 仿真近轴光线（至少10条）经过平凸透镜的焦距，与理论焦距值进行对比，得出误差大小。 ② 已知透镜的结构参数为101=r ，0.11=n ，51=d ，5163.121==' n n （K9玻璃）， 502-=r ，0.12=' n ，物点A 距第一面顶点的距离为100，由A 点计算三条沿光轴夹角分别为10、20、 30的光线的成像。试用Matlab 对以上三条光线光路和近轴光线光路进行仿真，并得出实际光线的球差大小。 ③ 设半径为1mm 的平面波经凸面曲率半径为25mm ，中心厚度3mm 的平凸透镜。用matlab 仿 真平面波在透镜几何焦平面上的聚焦光斑强度分布，计算光斑半径。并与理论光斑半径值进行对比，得出误差大小。（方法：采用波动理论，利用基尔霍夫—菲涅尔衍射积分公式。） 2、用MATLAB 仿真平行光束的衍射强度分布图样。（夫朗和费矩形孔衍射、夫朗和费圆孔衍射、夫朗和费单缝和多缝衍射。） 3、用MATLAB 仿真厄米—高斯光束在真空中的传输过程。（包括三维强度分布和平面的灰度图。） 4、（补充题）查找文献，掌握各类空心光束的表达式，采用费更斯-菲涅尔原理推导各类空心光束在真空中传输的光强表达式。用matlab 对不同传输距离处的光强进行仿真。 三、理论推导部分 第一大题 （1）十条近轴光线透过透镜时，理想情况下光线汇聚透镜的焦点上，焦点到像方主平面的距离为途径的焦距F ，但由于透镜的折射率和厚度会影响光在传输过程中所走的路径（即光程差Δ）。在用MATLAB 仿真以前先计算平行光线的传输路径。，R 为透镜凸面的曲率半径，h 为入射光线的高度，θ1为入射光线与出射面法线的夹角，θ2为出射光线与法线的夹角，n 为透镜材料的折射率。设透镜的中心厚度为d ，则入射光线经过透镜的实际厚度为：L=(R-d) 光线的入射角为：sinq1=h/R 折射角度满足：sinq2=nsinq1 而实际的光束偏折角度为：θ2-θ1。 由此可以看出，当平行光线照射透镜时，在凸面之前光线平行于光轴，在凸面之后发生了偏折，于光轴交汇一点，这一点成为焦点f ，折线的斜率为（-tan(θ2-θ1)）。 （2）根据题意可得，本题所讨论的是与光轴夹角不同的三条光线，经过透镜的两次反射后的成像问题。利用转面公式计算。
实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间范围0---Tn 。 3、),(T sys step ；表示时间范围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：) ()()()(1)(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.
《计算机仿真及应用》课程设计报告书 学号：08057102，08057127 班级：自动化081 姓名陈婷，万嘉
目录 一、设计思想 二、设计步骤 三、调试过程 四、结果分析 五、心得体会 六、参考文献
选题一、 考虑如下图所示的电机拖动控制系统模型，该系统有双输入，给定输入)(t R 和负载输入)(t M 。 1、 编制MATLAB 程序推导出该系统的传递函数矩阵。 2、 若常系数增益为：C 1＝Ka ＝Km ＝1，Kr ＝3，C2＝0.8，Kb ＝1.5，时间常数T 1＝5， T 2＝0.5，绘制该系统的根轨迹、求出闭环零极点，分析系统的稳定性。若)(t R 和)(t M 分别为单位阶跃输入，绘制出该系统的阶跃响应图。（要求C 1，Ka ，Km ，Kr ，C2，Kb ， T 1，T 2所有参数都是可调的） 一．设计思想 题目分析： 系统为双输入单输出系统，采用分开计算，再叠加。 要求参数均为可调，而matlb 中不能计算未赋值的函数，那么我们可以把参数设置为可输入变量，运行期间根据要求赋值。 设计思路： 使用append 命令连接系统框图。 选择‘参数=input('inputanumber:')’实现参数可调。 采用的方案： 将结构框图每条支路稍作简化，建立各条支路连接关系构造函数，运行得出相应的传递函数。 在得出传递函数的基础上，使用相应的指令求出系统闭环零极点、画出其根轨迹。 通过判断极点是否在左半平面来编程判断其系统是否稳定。 二．设计步骤 （1）将各模块的通路排序编号
（2）使用append命令实现各模块未连接的系统矩阵 （3）指定连接关系 （4）使用connect命令构造整个系统的模型 三．调试过程 出现问题分析及解决办法： 在调试过程出现很多平时不注意且不易寻找的问题，例如输入的逗号和分号在系统运行时不支持中文格式，这时需要将其全部换成英文格式，此类的程序错误需要细心。 在实现参数可调时初始是将其设为常量，再将其赋值进行系统运行，这样参数可调性差，后用‘参数=input('inputanumber:')’实现。 最后是在建立通路连接关系时需要细心。 四．结果分析 源代码： Syms C1 C2 Ka Kr Km Kb T1 T2 C1=input('inputanumber:') C2=input('inputanumber:') Ka=input('inputanumber:') Kr=input('inputanumber:') Km=input('inputanumber:') Kb=input('inputanumber:') T1=input('inputanumber:') T2=input('inputanumber:') G1=tf(C1,[0 1]); G2=tf(Ka*Kr,[0 1]); G3=tf(Km,[T1 1]); G4=tf(1,[T2 1]); G5=tf(1,[1 0]); G6=tf(-C2,1); G7=tf(-Kb,1); G8=tf(-1,1); Sys=append(G1,G2,G3,G4,G5,G6,G7,G8) Q=[1 0 0;2 1 6;3 2 7;4 3 8;5 4 0;6 5 0;7 4 0;8 0 0;]; INPUTS1=1; OUTPUTS=5; Ga=connect(Sys,Q,INPUTS1,OUTPUTS) INPUTS2=8; OUTPUTS=5; Gb=connect(Sys,Q,INPUTS2,OUTPUTS) rlocus(Ga)
课程设计 双闭环直流调速系统设计及仿真验证 学院年级：工程学院08级 组长：陈春明学号200830460102 08自动化1班成员一：陈木生学号 200830460103 08自动化1班 指导老师： 日期： 2012-2-28 华南农业大学工程学院
摘要 转速、电流双闭环调速系统是应用最广的直流调速系统，由于其静态性能良好，动态响应快，抗干扰能力强，因而在工程设计中被广泛地采用。现在直流调速理论发展得比较成熟，但要真正设计好一个双闭环调速系统并应用于工程设计却有一定的难度。 Matlab是一高性能的技术计算语言，具有强大的科学数据可视化能力，其中Simulink具有模块组态简单、性能分析直观的优点，方便了系统的动态模型分析。应用Simulink来研究双闭环调速系统，可以清楚地观察每个时刻的响应曲线，所以可以通过调整系统的参数来得出较为满意的波形，即良好的性能指标，这给分析双闭环调速系统的动态模型带来很大的方便。 本研究采用工程设计方法，并利用Matlab协助分析双闭环调速系统，依据自动控制系统快、准、稳的设计要求，重点分析系统的起动过程。 关键词：双闭环直流调速 Simulink 自动控制
目录 1、直流电机双闭环调速系统的结构分析．．．．．．．．．．．．．．．．．．．．．．． 1.1 双闭环调速系统的组成．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 1.2 双闭环调速系统的结构．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 2 、建立直流电机双闭环调速系统的模型．．．．．．．．．．．．．．．．．．．．．．．．．．．． 2.1 小型直流调速系统的指标及参数．．．．．．．．．．．．．．．．．．．．．．．．． 2.2 电流环设计．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 2.3 转速环设计．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 3、直流电动机双闭环调速系统的MATLAB仿真．．．．．．．．．．．．．．．．．．．． 3.1 系统框图的搭建．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 3.2 PI控制器参数的设置．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 3.3 仿真结果．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 4、结论与总结．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 5、参考资料．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．
MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************
实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：
1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：
matlab音频降噪课程设计报告
燕山大学 医学软件课程设计说明书 题目：基于MATLAB巴特沃斯滤波器的音频去噪的GUI设计 学院（系）：电气工程学院 年级专业： 13级生物医学工程 2 班 学号： 130103040041 学生姓名：魏鑫 指导教师：许全盛 1
院（系）：电气工程学院基层教学单位：生物医学工程系 学号130103040041 学生 姓名 魏鑫 专业（班 级） 13级生 物医学 工程2 班 设计 题目 基于MATLAB音频去噪的GUI设计设 计 技术参数通带截止频率fp=2700;阻带截止频率fs=3000;采样频率FS=48000; 通带衰减不大于1dB;阻带衰减不小于10dB; 设计要求1.实现用MATLAB导入音频； 2.对音频进行频谱分析； 3.设计滤波器去噪并对含噪信号进行滤 2
波并进行功率谱分析； 4.设计能实现上述功能的GUI； 工作量1.完成音频录入及频谱分析相关程序的编写与调试； 2.设计滤波器去噪； 3.用MATLAB软件做GUI界面的设计； 工作计划11.21-11.24 MATLAB软件中GUIDE 工具箱的使用 11.25-11.29 各处理算法模块的编程实现 11.30-12.1 整体程序联调 12.2 撰写课程设计说明书,答辩 参考资料 1. 陈怀琛吴大正 MATLAB及在电子信息课程中的应用[M] 北京电子工业出版社 2006. 章节2.4； 2. 陈亚勇 MATLAB信号处理详解[M] 北京：人民邮电出版社 2000. 第十 3
章； 3.张康刘雅基于Matlab的巴特沃斯 数字低通滤波器的设计[J] 计算机与现代化 2007年 12期 98-100页 指导 教师签字许全盛 基层教学单 位主任签字 彭勇 目录 一、设计目的意义 (1) 1.1绪论 (1) 1.2设计目的 (1) 1.3意义 (1) 二、设计内容 (2) 2.1 设计原理 (2) 2.2 设计内容 (2) 三、设计过程及结果分析 (3) 3.1 设计步骤 (3) 4
北华大学 《MATLAB仿真》课程设计 姓名： 班级学号： 实习日期： 辅导教师：
前言 科学技术的发展使的各种系统的建模与仿真变得日益复杂起来。如何快速有效的构建系统并进行系统仿真，已经成为各领域学者急需解决的核心问题。特别是近几十年来随着计算机技术的迅猛发展，数字仿真技术在各个领域都得到了广泛的应用与发展。而MATLAB作为当前国际控制界最流行的面向工程和科学计算的高级语言，能够设计出功能强大、界面优美、稳定可靠的高质量程序，而且编程效率和计算效率极高。MATLAB环境下的Simulink是当前众多仿真软件中功能最强大、最优秀、最容易使用的一个系统建模、仿真和分析的动态仿真环境集成工具箱，并且在各个领域都得到了广泛的应用。 本次课程设计主要是对磁盘驱动读取系统校正部分的设计，运用自动控制理论中的分析方法，利用MATLAB对未校正的系统进行时域和频域的分析，分析各项指标是否符合设计目标,若有不符合的，根据自动控制理论中的校正方法，对系统进行校正，直到校正后系统满足设计目标为止。我组课程设计题目磁盘驱动读取系统的开环传递函数为是设计一个校正装置，使校正后系统的动态过程超调量δ%≤7%，调节时间ts≤1s。 电锅炉的温度控制系统由于存在非线性、滞后性以及时变性等特点,常规的PID控制器很难达到较好的控制效果。考虑到模糊控制能对复杂的非线性、时变系统进行很好的控制, 但无法消除静态误差的特点, 本设计将模糊控制和常规的PI D控制相结合, 提出一种模糊自适应PID控制器的新方法。并对电锅炉温度控制系统进行了抗扰动的仿真试验, 结果表明, 和常规的PI D控制器及模糊PI D复合控制器相比,模糊自适应PI D控制改善了系统的动态性能和鲁棒性, 达到了较好的控制效果。
南京工程学院 课程设计说明书（论文）题目模拟信号的数字化 课程名称Matlab通信仿真设计 院(系、部、中心)通信工程学院 专业电子信息工程（传感网） 班级 学生姓名X X X 学号 2 0 8 1 1 0 7 3 2 设计地点信息楼C 216 指导教师潘子宇
设计起止时间：2014年1月10日至2014年 1 月14日
目录 一、内容摘要 (1) 二、课程设计目的和要求 (2) 三、课程设计任务 (2) 四、课程设计软件介绍 (3) 五、课程设计原理 (4) 六、PCM编码及仿真参数设置 (9) 七、PCM解码及仿真参数设置 (11) 八、PCM串行传输模型及仿真参数设置 (13) 九、课程设计成品图 (14) 十、SCOPE端的最终波形图 (14) 十一、主要参考文献 (15)
十二、总结与体会 (15) 一、内容摘要 MATLAB软件是矩阵实验室的简称，是美国M a t h W or k s公司出品的商业数学软件, 可用于算法开发、数据分析以及数值计算的高级技术计算语言和交互式环境, 广泛用于数字信号分析，系统识别，时序分析与建模， 神经网络、动态仿真等方面有着广泛的应用。主要包括MATLAB和Simulink两大部分。Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink。 Simulink是MATLAB中的一种可视化仿真工具，是一种基于MATLAB的框图设计环境，是实现动态系统建模、仿真和分析的一个软件包，被广泛应用于线性系统、非线性系统、数字控制及数字信号处理的建模和仿真中。Simulink可以用连续采样时间、离散采样时间或两种混合的采样时间进行建模，它也支持多速率系统，也就是系统中的不同部分具有不同的采样速率。为了创建动态系统模型，Simulink提供了一个建立模型方块图的图形用户接口(GUI) ，这个创建过程只需单击和拖动鼠标操作就能完成，它提供了一种更快捷、直接明了的方式，而且用户可以立即看到系统的仿真结果。 Simulink是用于动态系统和嵌入式系统的多领域仿真和基于模型的设计工具。对各种时变系统，包括通讯、控制、信号处理、视频处理和图像处理系统，Simulink提供了交互式图形化环境和可定制模块库来对其进行设计、仿真、执行和测试。.
实验报告 课程名称：控制理论（乙） 指导老师：林峰 成绩：__________________ 实验名称：MATLAB 仿真实验 实验类型：________________同组学生姓名：__________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 实验九 控制系统的时域分析 一、 实验目的： 1．用计算机辅助分析的办法，掌握系统的时域分析方法。 2．熟悉Simulink 仿真环境。 二、实验原理及方法： 系统仿真实质上就是对系统模型的求解，对控制系统来说，一般模型可转化成某个微分方程或差分方程表示，因此在仿真过程中，一般以某种数值算法从初态出发，逐步计算系统的响应，最后绘制出系统的响应曲线，进而可分析系统的性能。控制系统最常用的时域分析方法是，当输入信号为单位阶跃和单位冲激函数时，求出系统的输出响应，分别称为单位阶跃响应和单位冲激响应。在MATLAB 中，提供了求取连续系统的单位阶跃响应函数step ，单位冲激响应函数impulse ，零输入响应函数initial 等等。 二、实验内容： 二阶系统，其状态方程模型为 ? 1x -0.5572 -0.7814 1x 1 = + u ? 2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 四、实验要求： 1．编制MATLAB 程序，画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应； （1）画出系统的单位阶跃响应曲线； A=[-0.5572 -0.7814;0.7814 0 ]; B=[1;0];
Matlab程序设计任务书 分院（系）信息科学与工程专业 学生姓名学号 设计题目车牌识别系统设计 内容及要求： 车牌定位系统的目的在于正确获取整个图像中车牌的区域，并识别出车牌号。通过设计实现车牌识别系统，能够提高学生 分析问题和解决问题的能力，还能培养一定的科研能力。 1.牌照识别系统应包括车辆检测、图像采集、牌照识别等几 部分。 2.当车辆检测部分检测到车辆到达时,触发图像采集单元，采 集当前的视频图像。 3.牌照识别单元对图像进行处理，定位出牌照位置，再将牌 照中的字符分割出来进行识别，然后组成牌照号码输出。 进度安排： 19周：Matlab环境熟悉与基础知识学习 19周：课程设计选题与题目分析 20周：程序设计编程实现 20周：课程设计验收与答辩 指导教师（签字）： 年月日学院院长（签字）： 年月日 目录
一．课程设计目的 (3) 二．设计原理 (3) 三．详细设计步骤 (3) 四. 设计结果及分析 (18) 五. 总结 (19) 六. 设计体会 (20) 七. 参考文献 (21) 一、课程设计目的 车牌定位系统的目的在于正确获取整个图像中车牌的区域，并识别出车牌号。通过
设计实现车牌识别系统，能够提高学生分析问题和解决问题的能力，还能培养一定的科研能力。 二、设计原理: 牌照自动识别是一项利用车辆的动态视频或静态图像进行牌照号码、牌照颜色自动识别的模式识别技术。其硬件基础一般包括触发设备、摄像设备、照明设备、图像采集设备、识别车牌号码的处理机等，其软件核心包括车牌定位算法、车牌字符分割算法和光学字符识别算法等。某些牌照识别系统还具有通过视频图像判断车辆驶入视野的功能称之为视频车辆检测。一个完整的牌照识别系统应包括车辆检测、图像采集、牌照识别等几部分。当车辆检测部分检测到车辆到达时触发图像采集单元，采集当前的视频图像。牌照识别单元对图像进行处理，定位出牌照位置，再将牌照中的字符分割出来进行识别，然后组成牌照号码输出。 三、详细设计步骤: 1. 提出总体设计方案: 牌照号码、颜色识别 为了进行牌照识别，需要以下几个基本的步骤： a.牌照定位，定位图片中的牌照位置；