matlab模拟
MATLAB模拟太阳系实验报告
小组成员：邹智康U201313766
王思懋U201313475
陈桃斌U201313531
杨飞U201313530
陈铁U201313522 班级：信卓1301
指导老师：文灏
一、实验目的
1、熟悉MATLAB的基本功能，了解MATLAB的编译环
境。
2、掌握MATLAB的基本绘图函数，熟悉绘图函数的用
法、简单图形标注、简单颜色鉴定。
3、通过组内分工、协商合作，增强组员自主学习的能
力，加强组员间团结合作的能力。
二、预期目标
1、绘制太阳系三维动态图形。
2、模拟八大行星和彗星绕太阳运行动态轨道。
3、模拟月球绕地球运动.小星星带及土星环。
4、模拟出大致的行星大小.轨道大小及运行周期的比例
三、实验前的准备工作
1、对MATLAB的初步了解：MATLAB（矩阵实验室）
是MATrix LABoratory的缩写，是一款由美国The
MathWorks公司出品的商业数学软件。MATLAB是
一种用于算法开发、数据可视化、数据分析以及
数值计算的高级技术计算语言和交互式环境。除
了矩阵运算、绘制函数/数据图像等常用功能外，
MATLAB还可以用来创建用户界面及与调用其它
语言（包括C，C++和FORTRAN）编写的程序。
2、组内进行多次讨论，并明确组员分工：陈铁、王
思懋两人负责前期资料(包括太阳系的各种所需
数据、MATLAB的编译环境等)的收集，陈桃斌负
责ppt的制作，杨飞负责主体代码，邹智康负责
实验报告的书写。由杨飞将主题代码写好后，再
把其余代码分块分给其余组员各自书写，以保证
每个组员的参与度。
四、实验成果展示
1、全部代码
2、%3?ê?éè??
3、figure('name','D??μ?Yê?');%éè??±êìa??×?
4、pausetime=.01;%éè???Yí￡ê±??
5、set(gca,'xlim',[-50 50],'ylim',[-50 30],'zlim',[-50 50]);
6、set(gcf,'doublebuffer','on') %??3y???ˉ
7、xlabel('x?á'),ylabel('y?á'),zlabel('z?á');
8、axis equal;
9、grid on;
10、view([3 5 2]);
11、hold on
12、a=[8.5 12.5 20 30 50 60 80 100 90 1.5];%3¤?á
13、b=[8 12 18 26 45 55 70 90 30 1.5];%?ì?á
14、omga=[4 1.25 1 0.5 0.1 0.05 0.1 0.09 0.1 5];%???ù?è
15、r=[0.35 0.8 0.8 0.5 3 2.5 1.5 1.5 0.5 0.35];%?òì?°???
?°°????a??ó|DDD?￡?μú?????a??D?￡?μúê????a???ò
16、c=sqrt(a.^2-b.^2);h=pi/18;h1=pi/10;f=pi/9;g=pi/8;g1=pi/
6;
17、aby=[h h 0;h1 h 0;h h 0;h h 0;h h 0;h h 0;h g 0;h h h;g 0
g;0 g1 0];%????1ìμà??????ò?éè??
18、
19、%D?DDD?′?ó?íáD??·éè??
20、runu=35:0.5:40;theta=(0:0.05*pi:2*pi)';runa=2.8:0.4:5.6
;
21、xx=cos(theta)*runu+20;
22、yy=0.9*sin(theta)*runu;
23、zz=-0.17*xx-0.17*yy; plot3(xx,yy,zz,':k'); %D?DD′?éè??
24、hx=cos(theta)*runa;hy=sin(theta)*runa;hz=-0.1*hx-0.2*hy
; %color={'y','m','b','m','r','c','b','b'};
25、
26、%ì???éè??
27、[X,Y,Z]=sphere(40);
28、surf(5*X,5*Y,5*Z);colormap(autumn) %éè??ì???
29、light ('position',[1 0 2],'style','infinite')
30、lighting phong
31、material shiny
32、t=0:0.01*pi:50*pi; t';
33、num=length(a);
34、
35、%1ìμà·?3ì
36、for n=1:num
37、x(:,n)=a(n)*cos(omga(n)*t)+c(n);
38、y(:,n)=b(n)*sin(omga(n)*t);
39、z(:,n)=0*t;
40、xuanz(:,:)=[1 0 0;0 cos(aby(n,1)) -sin(aby(n,1));0
sin(aby(n,1)) cos(aby(n,1))]*[cos(aby(n,2)) 0
sin(aby(n,2));0 1 0;-sin(aby(n,2)) 0
cos(aby(n,2))]*[cos(aby(n,3)) -sin(aby(n,3))
0;sin(aby(n,3)) cos(aby(n,3)) 0;0 0 1];
41、xyz(:,:)=[x(:,n) y(:,n) z(:,n)]*xuanz(:,:);
42、x(:,n)=xyz(:,1);
43、y(:,n)=xyz(:,2);
44、z(:,n)=xyz(:,3);
45、if n<=9
46、p(n)=surf(r(n)*X+x(1,n),r(n)*Y+y(1,n),r(n)*Z+z(1,n));sh
ading interp
47、plot3(x(:,n),y(:,n),z(:,n),'-k');%?-3??ùóD1ì?￡??
48、else
49、p(n)=surf(r(n)*X+x(1,n)+x(1,3),r(n)*Y+y(1,n)+y(1,3),r(n
)*Z+z(1,n)+z(1,3));shading interp
50、pmuw(n-9)=plot3(x(:,n)+x(1,3),y(:,n)+y(1,3),z(:,n)+z(1,
3),'-k');%???ò1ìμàéè??
51、end
52、end
53、
54、%íáD??·
55、for n=1:length(runa)
56、ph(n)=plot3(hx(:,n)+x(1,6),hy(:,n)+y(1,6),hz(:,n)+z(1,6
),'-c');
57、end
58、
59、%DDD???é?
60、set(p(1),'facecolor','y');set(p(2),'facecolor','m');set
(p(3),'facecolor','b');set(p(4),'facecolor','m');
61、set(p(5),'facecolor','r');set(p(6),'facecolor','c');set
(p(7),'facecolor','b');set(p(8),'facecolor','b');
62、set(p(9),'facecolor','r');set(p(10),'facecolor','k');
63、
64、%?-3?í?D?
65、for m=1:5000
66、for n=1:num
67、if n<=9
68、set(p(n),'xdata',r(n)*X+x(m,n),'ydata',r(n)*Y+y(m,n),'z
data',r(n)*Z+z(m,n));%?ùóDμ??′ê±????
69、else
70、set(p(n),'xdata',r(n)*X+x(m,n)+x(m,3),'ydata',r(n)*Y+y(
m,n)+y(m,3),'zdata',r(n)*Z+z(m,n)+z(m,3));
71、set(pmuw(n-9),'xdata',x(:,n)+x(m,3),'ydata',y(:,n)+y(m,
3),'zdata',z(:,n)+z(m,3));
72、end
73、end
74、
75、for n=1:length(runa)
76、set(ph(n),'xdata',hx(:,n)+x(m,6),'ydata',hy(:,n)+y(m,6)
,'zdata',hz(:,n)+z(m,6));%1a?·?′ê±????
77、end
78、
79、pause(pausetime); %?Yí￡ò??á
80、drawnow
81、end
2、分块诠释
(1)
这是初始的设置，包括一些必要数据的设置以及
后面要用到的参数的设置。
(2)
这是小行星带、土星带以及太阳的设置
(3)
这是土星环的设置。
(4)
这是对行星颜色的设置。
(5)后面两块分别是轨道方程的设置以及最终图
形的绘制及显示。
3、图形展示
五、实验总结
这次实验由于是第一次MATLAB的相关实验，我们对这个软件都不太熟悉，导致在过程中出现了很多的问题。当然对软件的不熟悉也导致了组内合作中出现了不好分工，出现争议的问题。还好，在组内成员的耐心研究，团结合作下，完成了此次实验。虽然这次没有做到很完美，但确实也付出了我们所有人不少的心血，希望下次会做的更好。
matlab通信系统仿真
通信系统基于（15,11）汉明编码的matlab仿真 clear m=4; n=2^m-1; k=11; msg1=randint(1,6000*k,[0,1])%产生信息序列 code0=vec2mat(msg1,k); code1=encode(msg1,n,k,'hamming/binary');%进行汉明纠错编码 code2=reshape(code1,90,1000)'; code3=zeros(1000,90);%设置零矩阵，以便储存交织后的序列 for i=1:1000 temp=code2(i,:); temp1=reshape(temp,15,6);%按15*6交织 code3(i,:)=reshape(temp1',1,90); end [row,cloums]=size(code3); code4=code3'; [row1,cloums1]=size(code4); code5=reshape(code4,6,15000)'; %将交织后的序列转换为6行，15000列矩阵msg2=zeros(15000,6); for j=2:6, msg2(:,j)=xor(code5(:,j),code5(:,(j-1))); end msg2(:,1)=code5(:,1);%此时得到的msg2为gray映射后的信号序列 code8=msg2';%转置，便于比特符号转换 A=[32,16,8,4,2,1]; %为比特符合转换所设的序列 msg3=A*code8;%生成符号序列 msg4=qammod(msg3,64);%将符号序列进行64QAM调制 [row2,cloums2]=size(msg4); dB=0:1:20; for k=1:length(dB), snr=10.^(dB(k)./10); %信噪比 sgma=sqrt(63./(6*snr));%标准差 b1=real(msg4)+sgma*randn(row2,cloums2);%分路叠加噪声 b2=imag(msg4)+sgma*randn(row2,cloums2);%分路叠加噪声 rx=complex(b1,b2); %————————量化判决——————————% for m=1:row2, for n=1:cloums2, if ((b1(m,n)<-6)) b11(m,n)=-7;
基于Matlab的光学衍射仿真
基于Matlab的光学衍射实验仿真 摘要 光学试验中衍射实验是非常重要的实验. 光的衍射是指光在传播过程中遇到障碍物时能够绕过障碍物的边缘前进的现象, 光的衍射现象为光的波动说提供了有力的证据. 衍射系统一般有光源、衍射屏和接受屏组成, 按照它们相互距离的大小可将衍射分为两大类, 一类是衍射屏与光源和接受屏的距离都是无穷远时的衍射, 称为夫琅禾费衍射, 一类是衍射屏与光源或接受屏的距离为有限远时的衍射称为菲涅尔衍射。 本文用Matlab软件对典型的衍射现象建立了数学模型，对衍射光强分布进行了编程运算，对衍射实验进行了仿真。最后创建了交互式GUI界面，用户可以通过改变输入参数模拟不同条件下的衍射条纹。 本文对于衍射概念、区别、原理及光强分布编程做了详细全面的介绍 关键字：Matlab；衍射；仿真;GUI界面;光学实验
Matlab-based Simulation of Optical Diffraction Experiment Abstract Optical diffraction experiment is a very important experiment. is the diffraction of light propagation of light in the obstacles encountered in the process to bypass the obstacles when the forward edge of the phenomenon of light diffraction phenomenon of the wave theory of light provides a strong Evidence. diffraction systems generally have light, diffraction screen and accept the screen composition, size according to their distance from each other diffraction can be divided into two categories, one is the diffraction screen and the light source and the receiving screen is infinity when the distance between the diffraction Known as Fraunhofer diffraction, one is diffraction screen and the light source or accept a limited away from the screen when the diffraction is called Fresnel diffraction. In this paper, Matlab software on a typical phenomenon of a mathematical model of diffraction, the diffraction intensity distribution of the programming operation, the diffraction experiment is simulated. Finally, create an interactive GUI interface, users can change the input parameters to simulate different conditions of the diffraction pattern. This concept of the diffraction, difference, intensity distribution of programming principles and a detailed comprehensive description Key word: matlab；diffraction; simulation; gui interface; optical experiment
基于matlab干涉系统仿真_
《工程光学》综合性练习一题目：基于matlab的干涉系统仿真 学院精密仪器与光电子工程学院 专业测控技术与仪器
综合练习大作业一 一、要求 3-4人组成小组，对下面给出的各题目利用Matlab等工具进行仿真。 二、仿真题目 1、对于杨氏双缝干涉，改变双缝的缝宽和缝间距，观察干涉图样变化 ①原理图 图中参数 光线波长：lam=500纳米； 双缝距离：d=0.1毫米；（可调） 双缝距接收屏距离：D=1米； 接收屏范围：xs：-0.005~0.005 ys：-0.005~0.005 光源振幅：AI=A2=1; （单位振幅，可调） ②matlab代码： clear; lam=500e-9; %设定波长lam（500纳米） d=0.5e-3; %设定两缝之间距离d（0.5毫米） D=1; %双缝到接收屏距离D（1米） A1=1; %初始两光源均为单位振幅 A2=1; xm=0.005; ym=xm; %接受屏的范围ym，xm(0.01*0.01矩形) n=1001; xs=linspace(-xm,xm,n); %用线性采样法生成两个一位数组xs，ys %（n为总点数） ys=linspace(-ym,ym,n); L1=sqrt((xs-d/2).^2+ys.^2+D^2);%光屏上点（xs,ys）距光源1距离r1 L2=sqrt((xs+d/2).^2+ys.^2+D^2);%光屏上点（xs,ys）距光源2距离r2 E1=A1./sqrt(L1).*exp(1i*L1*2*pi/lam);%光源1在接受屏上复振幅E1 E2=A2./sqrt(L2).*exp(1i*L2*2*pi/lam);%光源2在接受屏上复振幅E2 E=E1+E2; %复振幅叠加为合成振幅E
基于MATLAB的物理光学实验仿真平台构建
毕业设计(论文)开题报告题目：基于Matlab的物理光学实验仿真平台构建 院（系）光电工程学院 专业光信息科学与技术 班级120110 姓名闫武娟 学号120110127 导师刘王云 年月日
开题报告填写要求 1．开题报告作为毕业设计（论文）答辩委员会对学生答辩资格审查的依据材料之一。 此报告应在指导教师指导下，由学生在毕业设计（论文）工作前期内完成。2．开题报告内容必须按教务处统一设计的电子文档标准格式（可从教务处网页上下载）填写并打印（禁止打印在其它纸上后剪贴），完成后应及时交给指导教师审阅。3．开题报告字数应在1500字以上，参考文献应不少于15篇（不包括辞典、手册，其中外文文献至少3篇），文中引用参考文献处应标出文献序号，“参考文献”应按附件中《参考文献“注释格式”》的要求书写。 4．年、月、日的日期一律用阿拉伯数字书写，例：“2005年11月26日”。
这些仿真平台的使用不仅方便了教学，而且也使学生更容易理解物理光实验的基本原理，加深对理论知识的理解与记忆。 2.课题研究的主要内容和拟采用的研究方案、研究方法 2.1课题研究的主要内容 (1). 在光的干涉基本理论基础上，实现两束平面波、球面波的干涉实验，杨氏双缝和杨氏双孔干涉实验，平行平板的等倾干涉实验，楔形平板的等厚干涉实验，牛顿环干涉实验，迈克尔逊干涉实验以及平行平板的多光束干涉实验。 (2). 在菲涅尔衍射及夫琅和费衍射基本理论基础上，实现矩孔、单缝、圆孔、双缝、多缝、平面光栅及闪耀光栅的衍射实验。 2.2 研究方法及方案 物理光学实验可分为两大类：干涉与衍射。光的干涉有光源、干涉装置和干涉图形三个基本要素；衍射分为菲涅尔衍射和夫琅禾费衍射。光学领域的大部分图像及曲线分布都可以用MATLAB 软件加以计算和实现[16]， 以杨氏双缝干涉为例，简述实验方案 杨氏双缝干涉模型是典型的分波面干涉，其干涉装置图如图所示，用一个单缝与一个双缝，从同一波面上分出两个同相位的单色光，进而获得相干光源并观察分析干涉图样。 图1.1杨氏双缝干涉实验装置图 2.2.1数学建模 根据干涉的基本原理，点光源S 发出的光波经双缝分解为次波源S 1、S 2，这两个次波源发出的光波在空间相干叠加，继而在其后的接收屏形成一系列明暗相间的干涉条纹。 设入射光波波长为λ，两个次波源的强度相同，且间距为d （1）位相差的计算: 221)2 (y d x r ++ =222)2 - (y d x r +=（2.1） )(*12r r n -=?（2.2）
工程光学matlab仿真设计
工程光学仿真实验报告 1、氏双缝干涉实验 （1）氏干涉模型 氏双缝干涉实验装置如图1所示: S 发出的光 波射到光屏上的两个小孔S1 和S2 , S1 和S2 相 距很近,且到S 等距;从S1 和S2 分别发散出的光 波是由同一光波分出来的,所以是相干光波,它们在距离光屏为D 的屏幕上叠加,形成一定的干涉图 样。 图1.1 氏双缝干涉 假设S 是单色点光源,考察屏幕上某一点P ,从S1 和S2 发出的光波在该点叠加 产生的光强度为: I = I1 + I2 + 2 I1 I2 cos δ （1-1） 式中, I1 和I2 分别是两光波在屏幕上的光强度, 若实验装置中S1 和S2 两个缝 大小相等, 则有 I1 = I2 =I0 （1-2） δ= 2π（r2 - r1）/λ（1-3） （1-3） 2221)2/(D y d x r +++= （1-4） 2222)2/(D y d x r ++-= （1-5） 可得 xd r r 22 122=- （1-6） 因此光程差：12r r -=? （1-7） 则可以得到条纹的强度变化规律- 强度分布公式： ]/)([cos 1220λπd r r I I -= （1-8） （2）仿真程序 clear; Lambda=650; %设定波长,以Lambda 表示波长 Lambda=Lambda*1e-9; d=input('输入两个缝的间距 )'); %设定两缝之间的距离,以d 表示两缝之间距离 d=d*0.001; Z=0.5; %设定从缝到屏幕之间的距离,用Z 表示
yMax=5*Lambda*Z/d;xs=yMax; %设定y方向和x方向的围 Ny=101;ys=linspace(-yMax,yMax,Ny);%产生一个一维数组ys,Ny是此次采样总点数 %采样的围从- ymax到ymax,采样的数组命名为ys %此数组装的是屏幕上的采样点的纵坐标 for i=1:Ny %对屏幕上的全部点进行循环计算,则要进行Ny次计算L1=sqrt((ys(i)-d/2).^2+Z^2); L2=sqrt((ys(i)+d/2).^2+Z^2); %屏上没一点到双缝的距离L1和L2 Phi=2*pi*(L2-L1)/Lambda; %计算相位差 B(i,:)=4*cos(Phi/2).^2; %建立一个二维数组,用来装该点的光强的值 end%结束循环 NCLevels=255; %确定使用的灰度等级为255级 Br=(B/4.0)*NCLevels; %定标:使最大光强(4. 0)对应于最大灰度级(白色) subplot(1,4,1),image(xs,ys,Br); %用subplot创建和控制多坐标轴 colormap(gray(NCLevels)); %用灰度级颜色图设置色图和明暗 subplot(1,4,2),plot(B(:),ys); %把当前窗口对象分成2块矩形区域 %在第2块区域创建新的坐标轴 %把这个坐标轴设定为当前坐标轴 %然后绘制以( b (: ) , ys)为坐标相连的线title('氏双缝干涉'); （3）仿真图样及分析 a)双缝间距2mm b)双缝间距4mm
MATLAB实现通信系统仿真实例
补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样
MATLAB仿真教程
一、设计目的 通过运用MATLAB对函数进行Z域分析和单边带信号的调制与解调，使我们进一步加深对MATLAB的认识和运用，以实现以下目的： 1.本次试验进一步熟悉了MATLAB软件的使用方法及相关的操作。 2.对Z变换及其反变换函数在MATLAB中的调用有了掌握。 3.理论与实际的仿真相结合，更直观的看到结果。 4.观察了单边带信号调制与解调后的图像，加深认识。 二、设计原理 MATLAB是The MathWorks公司在1984年推出的一种商品化软件，它提供了大量丰富的应用函数，并且具有扩充的开放性结构。目前，该软件包涵盖了控制系统应用、数字信号处理、数字图像处理、通讯、神经网络、小波理论分析、优化与统计、偏微分方程、动态系统实时仿真等多学科专业领域。 其中单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同，产生单边带调制信号的方法有：滤波和相移法。 由于滤波法在技术上比较难实现所以在此我们将用相移法对单边带调制与解调系统进行讨论与设计。 三、设计内容和MATLAB图像
1、数字系统的响应 源代码如下： b=[0 1 2 1 0]; a=[1 -0.5 0 0.3 -0.005]; subplot(421);zplane(b,a); title('系统的零极点图'); subplot(422);impz(b,a,21); title('单位脉冲响应'); subplot(423);stepz(b,a,21); title('单位阶跃响应');
N=21;n=0:N-1; x=exp(-n); x0=zeros(1,N); y0=[1,-1]; xi=filtic(b,a,y0); y1=filter(b,a,x0,xi); xi0=filtic(b,a,0); y2=filter(b,a,x,xi0); y3=filter(b,a,x,xi); [h w]=freqz(b,a,21); subplot(424);stem(n,y1); title('零输入响应');grid on; subplot(425);stem(n,y2); title('零状态响应');grid on; subplot(426);stem(n,y3); title('系统的全响应');grid on; subplot(427);plot(w,abs(h)); title('幅频特性曲线');grid on; subplot(428);plot(w,angle(h)); title('相频特性曲线');grid on;
利用MATLAB模拟光学简单空间滤波系统
利用MATLAB 模拟光学简单空间滤波系统 摘要：阿贝成像原理是第一步在透镜的后焦面上得到物的空间频谱分布，第二步成像则是合频的过程，实则是两次傅立叶变换。利用阿贝-波特实验装置和空间滤波系统，可以对一幅光学图像进行光学信息处理。通过MATLAB 环境编写程序完成阿贝-波特实验和空间滤波的物理模型的构建并进行计算机模拟。 关键词：MATLAB ；阿贝成像原理；空间滤波；计算机模拟 引言： 早在1873年,阿贝（E ．Abbe,1840—1905）在德国蔡司光学器械公司研究如何提高显微镜的分辨本领问题时,就认识到相干成像的原理。空间滤波的主要目的是通过有意识地改变像的频谱,使像实现所希望的变化。光学信息处理是一个更为广阔的领域，它是基于光学频谱分析，利用傅里叶综合技术，通过空域或频域调制，借助空间滤波技术对光学信息进行处理的过程。阿贝提出的二次成像理论和20世纪初的阿贝—波特实验,已经为光学信息处理打下了一定的理论基础。 在阿贝成像理论的教学中,单纯依靠数学推演来讲解,效果不好,特别是空间频率、空间滤波等概念的形成有一定的困难。虽然可以通过演示阿贝- 波特实验来加强教学效果,但由于在普通教室难以完成演示实验,在实验室又受仪器、场地等方面的限制,实验现象不太理想。为此,我们设计出计算机模拟实验, 获得较好的模拟效果。在学习了解了阿贝成像原理的基础上，我们可以通过MATLAB 完成对阿贝-波特实验和空间滤波系统的计算机模拟，观察各种物体的空间频谱分布，设计各种不同的空间滤波器。 1.阿贝成像原理 在相干平行光照明下,显微镜的物镜成像可以分成两步：第一步即分频过程，由入射光经过物平面1P 发生衍射在物镜的后焦面2P 上形成夫琅禾费衍射图样;第二步称为合频或频谱综合过程，衍射图样作为新的子波源发出的球面波在像平面上相干叠加成像。相干光的成像过程本质上是两次傅立叶变换,第一次是将光场空间分布变成频率分布，第二次则是傅立叶逆变换,即将各频谱分量复合为像。如下图所示，为阿贝成像原理图。 L
matlab仿真模拟
中国地质大学长城学院 本科课程设计题目:双手协调机器人 系别信息工程系 学生姓名 专业电气工程及其自动化 学号 指导教师王密香 职称研究生 2015年12 月22日
双手协调机器人 摘要 多机器人的协同作业是制造业发展的必然要求，双臂机器人就是适应这一要求而开发出的一种新型机器人，相对于单臂机器人它可以大大增强机器人对复杂装配任务的适应性，同时可以提高工作空间的利用效率。当前大多数工业机器人的应用是为单臂机器人独自工作的能力准备的。一般地，单臂机器人只适合于刚性工件的操作，并受制于环境，随着现代工业的发展和科学技术的进步，对于许多任务而言单臂操作是不够的。为了适应任务的复杂性、智能性的不断提高以及系统柔顺性的要求而扩展为双手协调控制。即由两个单臂机器人相互协调、相互配合的去完成某种作业，但由于组成双手协调控制系统的是两个机器人它们不可能是两个单手机器人的简单组合，除了它们各自共同目标的控制实现外，它们相互间的协调控制以及对环境的适应性就成为组合的关键，这样双手协调控制机器人系统的进一步应用就受到了限制。而双臂机器人能完成对于人来说易于实现的功能，它比双手协调机器人更具有实用价值，它的高自律性以及学习性，能够适应许多环境，使其在工业生产、危险处理、国防、航天航空等方面运用广泛采用了Matlab/Simulink 仿真软件，分别用模块法和程序法对双手协调机器人系统进行了PID 控制器的校正仿真设计。 关键词：双手协调机器人；Matlab/Simulink；PID 校正；仿真 Using the Matlab/Simulink simulation software, using method of module and the procedural law on hands coordinate robot system has carried on the correction for the simulation of PID controller design simulation results show that these two kinds of design method is not only convenient and quick, and the correction effect is satisfactory to people Keywords：Hands coordinate robot; Matlab/Simulink; PID correction; The simulation
圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学(20200607000913)
工程光学综合练习-----圆孔、矩孔的菲涅尔衍射模拟
圆孔和矩孔的菲涅尔衍射模拟 一、原理 由惠更斯-菲涅尔原理可知接收屏上的P点的复振幅可以表示为 其中为衍射屏上的复振幅分布， 为倾斜因子。根据基尔霍夫对此公式的 完善，有 设衍射屏上点的坐标为（x1, y1），接收屏上点的坐标为（x, y），衍射屏与接收屏间距离为z1，当满足菲涅尔近似条件时，即 此时可得到菲涅尔衍射的计算公式 把上式指数项中的二次项展开，并改写成傅里叶变换的形式，可以写成上式为菲涅尔衍射的傅里叶变换表达式，它表明除了积分号前面的一个与 x1、y1无关的振幅和相位因子外，菲涅尔衍射的复振幅分布是孔径平面的复振 幅分布和一个二次相位因子乘积的傅里叶变换。 相对于夫琅和费衍射而言，菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中，输入变量和输出变量分别为衍射孔径平面的光场分布和观察平面的光场 以及光强分布，考虑到这三个量都是二维分布,而且Matlab主要应用于矩阵数值运算,所以本程序选择用二维矩阵来存储衍射孔径平面和观察平面的场分布,并分别以矩阵的列数和行数来对应平面的直角坐标值(x, y)以及(x1, y1)。 二、圆孔菲涅尔衍射 用MATLAB分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同，考虑到运算量的要求，采样点数不能过多，所以每个屏的x和y方向各取200到300点进行运算。根据式（4），选取合适的衍射屏和接收屏尺寸和相距的
距离，模拟结果如下： 取典型的He-Ne激光器波长=632.8nm，固定衍射屏和接收屏尺寸和相距的 距离，分别取不同的圆孔半径，得到以下三组衍射图样，其圆孔半径分别为12mm，20mm，50mm 图 1（r=12mm） 图 2（r=20mm）
MATLAB与系统仿真
学习中心/函授站_ 成都学习中心 姓名赵洪学号7020140122093 西安电子科技大学网络与继续教育学院 2015学年上学期 《MATLAB与系统仿真》期末考试试题 （综合大作业） 考试说明： 1、大作业于2015年4月3日公布，2015年5月9日前在线提交； 2、考试必须独立完成，如发现抄袭、雷同、拷贝均按零分计。 3、程序设计题（三（8，10））要求写出完整的程序代码，并在matlab软件环境调试并运行通过，连同运行结果一并附上。 一、填空题（1? ×25=25?） 1、Matlab的全称为矩阵实验室。 2、在Matlab编辑器中运行程序的快捷键是：F5 。 3、Matlab的工作界面主要由以下五个部分组成，它们分别是：菜单栏、 工具栏、当前工作目录窗口、工作空间管理窗口和命令窗口。 4、在Matlab中inf表示：无穷大；clc表示：清空命令窗口中的显示内容；more表示：在命令窗口中控制其后每页的显示内容行数；who表示：查阅Matlad内存变量名；whos表示：列出当前工作空间所有变量。 5、在Matlab命令窗口中运行命令Simulink 可以打开Simulink模块库浏览器窗口。 6、求矩阵行列式的函数：det ；求矩阵特征值和特征向量的函数eig 。 7、Matlab预定义变量ans表示：没有指定输出变量名；eps表示：系统精度 ；nargin表示：函数输入参数的个数。 8、Matlab提供了两种方法进行程序分析和优化，分别为：通过Profiler工具优化和通过tic和toc函数进行优化。 9、建立结构数组或转换结构数组的函数为：struct ； 实现Fourier变换在Matlab中的对应函数为：fourier() ；Laplace变换的函数：Laplace() 。
工程光学matlab仿真
工程光学仿真实验报告1、杨氏双缝干涉实验 （1）杨氏干涉模型 屏 图 , 0（1-8） 2 1 （2）仿真程序 clear; Lambda=650; %设定波长,以Lambda表示波长 Lambda=Lambda*1e-9; d=input('输入两个缝的间距 )'); %设定两缝之间的距离,以d表示两缝之间距离 d=d*0.001; Z=0.5; %设定从缝到屏幕之间的距离,用Z表示 yMax=5*Lambda*Z/d;xs=yMax; %设定y方向和x方向的范围
Ny=101;ys=linspace(-yMax,yMax,Ny);%产生一个一维数组ys,Ny 是此次采样总点数 %采样的范围从- ymax 到ymax,采样的数组命名为ys %此数组装的是屏幕上的采样点的纵坐标 for i=1:Ny %对屏幕上的全部点进行循环计算,则要进行Ny 次计算 L1=sqrt((ys(i)-d/2).^2+Z^2); L2=sqrt((ys(i)+d/2).^2+Z^2); %屏上没一点到双缝的距离L1和L2 Phi=2*pi*(L2-L1)/Lambda; %计算相位差 B(i,:)=4*cos(Phi/2).^2; %建立一个二维数组,用来装该点的光强的值 end %结束循环 NCLevels=255; %确定使用的灰度等级为255级 Br=(B/4.0)*NCLevels; %定标:使最大光强(4. 0)对应于最大灰度级(白色) subplot(1,4,1),image(xs,ys,Br); %用subplot 创建和控制多坐标轴 colormap(gray(NCLevels)); %用灰度级颜色图设置色图和明暗 subplot(1,4,2),plot(B(:),ys); %把当前窗口对象分成2块矩形区域 %在第2块区域创建新的坐标轴 %把这个坐标轴设定为当前坐标轴 %然后绘制以( b (: ) , ys)为坐标相连的线 title('杨氏双缝干涉'); （3）仿真图样及分析 a)双缝间距2mm b)双缝间距4mm c)双缝间距6mm d)双缝间距8mm 图1.2改变双缝间距的条纹变化 由上面四幅图可以看出，随着双缝之间的距离增大，条纹边缘坐标减小，也就是条纹 间距减小，和理论公式d D e /λ=推导一致。如果增大双缝的缝宽，会使光强I 增加，能够 看到条纹变亮。 二、杨氏双孔干涉实验 1、杨氏双孔干涉 杨氏双孔干涉实验是两个点光源干涉实 验的典型代表。如图2所示。当光穿过这两个 离得很近小孔后在空间叠加后发生干涉, 并 在像屏上呈现出清晰的明暗相间的条纹。 由 于双孔发出的波是两组同频率同相位的球面 波, 故在双孔屏的光射空间会发生干涉。 于是, 在图2中两屏之间的空间里, 如果一点P 处于 两相干的球面波同时到达 波 峰 (或波谷)的位置, 叠加后振幅达到最高, 图2.1 杨氏双孔干涉 表现为干涉波的亮点; 反之, 当P 处处于一个球面波的波峰以及另一个球面波的波谷时候， 叠加后振幅为零，变现是暗纹。
MATLAB仿真实验全部
实验一 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.
Matlab数字衍射光学实验讲义(一)
实验注意事项（必读） 1．没有弄清楚实验内容者，禁止接触实验仪器。 2．注意激光安全。绝对不可用眼直视激光束，或借助有聚光性的光学组件观察激光束，以免损伤眼睛。 3．注意用电安全。He-Ne激光器电源有高压输出，严禁接触电源输出和激光头的输入端，避免触电。 4．注意保持卫生。严禁用手或其他物品接触所有光学元件（透镜、反射镜、分光镜等）的光学表面；特别是 在调整光路中，要避免手指碰到光学表面。 5．光学支架上的调整螺丝，只可微量调整。过度的调整，不仅损坏器材，且使防震功能大减。 6．实验完成后，将实验所用仪器摆放整齐，清理一下卫生。
Matlab数字衍射光学实验一 计算机仿真过程是以仿真程序的运行来实现的。仿真程序运行时，首先要对描述系统特性的模型设置一定的参数值，并让模型中的某些变量在指定的范围内变化，通过计算可以求得这种变量在不断变化的过程中，系统运动的具体情况及结果。仿真程序在运行过程中具有以下多种功能： 1）计算机可以显示出系统运动时的整个过程和在这个过程中所产生的各种现象和状态。具有观测方便，过程可控制等优点； 2）可减少系统外界条件对实验本身的限制，方便地设置不同的系统参数，便于研究和发现系统运动的特性； 3）借助计算机的高速运算能力，可以反复改变输入的实验条件、系统参数，大大提高实验效率。因此．计算机仿真具有良好的可控制性(参数可根据需要调整)、无破坏性(不会因为设计上的不合理导致器件的损坏或事故的发生)、可复现性(排除多种随机因素的影响，如温度、湿度等)、易观察性(能够观察某些在实际实验当中无法或者难以观察的现象和难以实现的测量，捕捉稍纵即逝的物理现象，可以记录物理过程的每一个细节)和经济性(不需要贵重的仪器设备)等特点。 Matlab是MathWorks公司于1982年推出的一套高性能的数值计算和可视化软件。它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便、界面友好的用户环境。它还包括了ToolBox(工具箱)的各类问题的求解工具，可用来求解特定学科的问题。Matlab的长处在于数值计算，能处理大量的数据，而且效率比较高。MathWorths公司在此基础上开拓了符号计算、文字处理、可视化建模和实时控制能力，增强了Matlab的市场竞争力，使Matlab成为市场主流的数值计算软件。Matlab产品族支持概念设计、算法开发、建模仿真、实时实现的理想的集成环境。其主要功能有：数据分析、数值和符号计算、工程与科学绘图、控制系统设计、数字图像信号处理、财务工程，建模、仿真、原型开发，应用开发，图形用户界面。 在光学仪器设计和优化过程中，计算机的数值仿真已经成为不可缺少的手段。通过仿真计算，可以大幅度节省实验所耗费的人力物力，特别是在一些重复实验工作强度较大且对实验器材、实验环境等要求较苛刻的情况下。如在大型激光仪器的建造过程中，结合基准实验的仿真计算结果可为大型激光器的设计和优化提供依据。仿真光学实验也可应用于基础光学教学。光学内容比较抽象，如不借助实验，很难理解，如光的干涉、菲涅耳衍射、夫琅禾费衍射等。传统的光学实验需要专门的实验仪器和实验环境。其操作比较烦琐，误差大现象也不明显，对改变参数多次观察现象也多有不便。MATLAB是当今国际上公认的在科技领域方面最为优秀的应用软件和开发环境。利用它对光学实验仿真可避免传统实验中的缺点，强大的功能使光学实验变得简便准确。基于MATLAB的科学可视化功能对光学仿真实验现象进行计算机模拟的效果更加准确明显。 1.实验目的: 掌握基本的Matlab编程语言，了解其编程特点；模拟几种常用函数，了解其编程过程及图像显示命令函数，掌握Matlab画图方法；通过设计制作一系列光学研究物体掌握其编程方法；掌握光波的matlab编程原理及方法，初步了解Matlab
MATLABsimulink系统仿真分析仿真报告
仿真报告 课程名称：自动化技术导论 报告题目：MATLAB/simulink系统仿真分析 班级 姓名 学号 xxxxxx自动化学院 2016年4月 软件版本：MATLAB R2010b MATLAB强处理能力 MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而且经过了各种优化和容错处理。在通常情况下，可以用它来代替底层编程语言，如C和C++ 。在计算要求相同的情况下，使用MATLAB的编程工作量会大大减少。MATLAB的这些函数集包括从最简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。 MATLAB图形处理 MATLAB自产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程绘图。新版本的MATLAB对整个图形处理功能作了很大的改进和完善，使它不仅在一般数据可视化软件都具有的功能（例如二维曲线和三维曲面的绘制和处理等）方面更加完善，而且对于一些其他软件所没有的功能（例如图形的光照处理、色度处理以及四维数据的表现等），MATLAB同样表现了出色的处理能力。同时对一些特殊的可视化要求，例如图形对话等，MATLAB也有相应的功能函数，保证了用户不同层次的要求。另外新版本的MATLAB还着重在图形用户界面（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足。 MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱。一般来说，它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。领域，诸如数据采集、数据库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信号处理、图像处理、系统辨识、控制系统设计、LMI控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系统开发、定点仿真、DSP与通讯、电力系统仿真等，都在工具箱（Toolbox）家族中有了自己的一席之地。MATLAB程序接口
MATLAB仿真程序
窗型选择仿真程序： clear,clc bw=3e6;% 信号带宽 T=1e-4;%信号脉冲宽度 A=2;%信号幅度 fs=4*bw; lfft=round(T*fs);%采样点数 lfft=2^nextpow2(lfft); dt=1/fs;%采样间隔 f0=1e6; t=(0:lfft-1)*dt;%时域采样点 q=(0:lfft-1)*2*pi/lfft; s=A*exp(j*2*pi*f0*t+j*pi*bw*t.*t/T);%产生线性调频信号 S=(fft(s));%线性调频信号的傅立叶变换fft H=conj(S);%匹配滤波器的频率响应 Y=S.*H;%线性调频信号的频域匹配滤波输出 y=fftshift(ifft(Y));%线性调频信号的时域匹配滤波输出 %对chirp信号进行时域加权 h1=(triang(lfft))';%三角窗函数 s1=s.*h1;S1=fft(s1);H1=conj(S1); Y1=S1.*H1; y1=fftshift(ifft(Y1));%加三角窗后的线性调频信号的时域匹配滤波输出h2=(hanning(lfft))';%汉宁窗函数 s2=s.*h2;S2=fft(s2);H2=conj(S2); Y2=S2.*H2; y2=fftshift(ifft(Y2));%加汉宁窗后的线性调频信号的时域匹配滤波输出h3=(hamming(lfft))';%海明窗函数 s3=s.*h3;S3=fft(s3);H3=conj(S3); Y3=S3.*H3; y3=fftshift(ifft(Y3));%加海明窗后的线性调频信号的时域匹配滤波输出 figure; subplot(3,1,1), plot(t,real(s)),title('chirp signal'); subplot(3,1,2), plot(q,abs(S)),title('线性调频信号幅度谱'); subplot(3,1,3), plot(q,angle(S)),title('线性调频信号相位谱');
信息光学matlab仿真
%圆孔的夫琅禾费衍射： N=512; r=3; %衍射圆孔的半径 I=zeros(N,N); [m,n]=meshgrid(linspace(-N/16,N/16-1,N)); D=(m.^2+n.^2).^(1/2); I(find(D<=r))=1; subplot(1,2,1),imshow(I); title('生成的衍射圆孔'); % 夫琅禾费衍射的实现过程 L=500; [X,Y]=meshgrid(linspace(-L/2,L/2,N)); lamda_1=630; % 输入衍射波长; lamda=lamda_1/1e6 k=2*pi/lamda; z=1000000; % 衍射屏距离衍射孔的距离h=exp(1j*k*z)*exp((1j*k*(X.^2+Y.^2))/(2*z))/(1j*lamda*z);%脉冲相应 H =fftshift(fft2(h)); %传递函数 B=fftshift(fft2(I)); %孔径频谱 G=fftshift(ifft2(H.*B)); subplot(1,2,2),imshow(log(1+abs(G)),[]); title('衍射后的图样'); figure meshz(X,Y,abs(G)); title('夫琅禾费衍射强度分布')
%单缝的夫琅禾费衍射： N=512; a=25; % 单缝的宽度 b=1000;% 单缝的长度 I=zeros(N,N); [m,n]=meshgrid(linspace(-N/4,N/4,N)); I(-a