求助四元数具体用法,有matlab程序最好
四元数是简单的超复数。
复数是由实数加上虚数单位 i 组成,其中i^2 = -1。
相似地,四元数都是由实数加上三个虚数单位 i、j、k 组成,而且它们有如下的关系: i^2 = j^2 = k^2 = -1, i^0 = j^0 = k^0 = 1 , 每个四元数都是 1、i、j 和 k 的线性组合,即是四元数一般可表示为a + bk+ cj + di,其中a、b、c 、d是实数。
对于i、j、k本身的几何意义可以理解为一种旋转,其中i旋转代表X轴与Y轴相交平面中X轴正向向Y轴正向的旋转,j旋转代表Z轴与X轴相交平面中Z轴正向向X轴正向的旋转,k旋转代表Y轴与Z轴相交平面中Y轴正向向Z轴正向的旋转,-i、-j、-k分别代表i、j、k旋转的反向旋转。
轨迹发生器是什么?
是一种程序,通过一些列数值模拟如导弹等一些飞行物体的轨迹,模拟输出载体在飞行过程中的位置、速度、姿态等导航参数,同时模拟给出在各飞行时刻载体上的惯导元件所测得的真实比例和姿态角速度信息将该比例和姿态角速度输入给捷联惯导导航解算软件,加上各种误差,即可模拟捷联惯导导航解算输出,包括位置、速度及姿态信息。
(部分信息摘自《机载武器轨迹发生器设计与仿真》)
matlab的几个干涉实验模拟
展开全部 function photo_diffraction%光学衍射仿真 矩形lmda=632.8e-9;% 波长xmax=0.05;% 观察屏所取范围a=1e-3;f=1;ymax=xmax;def=0.0001;x=-xmax:def:xmax;y=-ymax:def:ymax;lenm=length(x);lenn=length(y);for m=1:lenmfor n=1:lennalpha=pi*x(m)*a/(lmda*f);%0.5*k*l*a;beta=pi*y(n)*a/(lmda*f);%bb=0.5*k*w*b;I(m,n)=((sin(alpha))/(alpha))^2*((sin(beta))/(beta))^2;endendI=I/(max(max(I)));[X,Y]=meshgrid(x,y);figuremesh(X,Y,I);xlabel('x');ylabel('y');zlabel(' 光强');rotate3Dhold onfigureimshow(255*I);xlabel('x');ylabel('y'); %---------光学衍射仿真 单缝------------clearlam=500e-9;a= 1e-3; f=1;xm= 3*lam*f/a;nx= 51;xs=linspace(-xm,xm,nx);np=51;xp=linspace(0,a,np);for i=1:nxsinphi= xs(i)/f;alpha=2*pi*xp*sinphi/lam;sumcos=sum(cos(alpha));sumsin=sum(sin(alpha));B(i,:)=(sumcos^2+sumsin^2)/np^2;endN=255;Br=(B/max(B))*N;figure;subplot(1,2,1)image(xm,xs,Br);colormap(gray(N));subplot(1,2,2)plot(B,xs);%--------光学衍射仿真 多缝--------------clearlam=500e-9;N=2;a= 2e-4;z=5;d=5*a;xm=2*lam*z/a;y0=xm;n=1001;x0=linspace(-xm,xm,n);for i= 1: nsinphi=x0(i)/z;alpha=pi*a*sinphi/lam;beta=pi*d*sinphi/lam;B(i,:)=(sin(alpha)./alpha).^2.*(sin(N*beta)./sin(beta)).^2;B1=B/max(B);endNC=255;Br=(B/max(B))*NC;figure;subplot(1,2,1)image(y0,x0,Br);colormap( gray(NC) );subplot(1,2,2)plot(B1,x0);%--------光学衍射仿真 圆孔--------------clearlam=500;% '请输入光的波长:')lam=lam*1e-9;a=2e-3;f=1;m=200;ym=2000*lam*f;ys=linspace(-ym,ym,m);xs=ys;n=255;for i=1:mr=xs(i)^2+ys.^2;sinth=sqrt(r./(r+f^2));x=2*pi*a*sinth./lam;hh=(2*BESSELJ(1,x)).^2./x.^2;b(:,i)=(hh)'.*5000;endfiguresubplot(1,2,1)image(xs,ys,b)colormap(gray(n))subplot(1,2,2)b(:,m/2);plot(ys,b(:,m/2))
转载请注明出处51数据库 » 用matlab软件模拟捷联惯导的飞行模拟