jingque.m

内燃机转子仿真

clear all
clc
close all
t=40;
h=0.00281;
dn=rem(t,h);
ddn=t-dn;
n1=ddn/h;
i=1;
A=[1.000/3.000^0.5 0.500*(2.000/3.000)^0.5;1.000/3.000^0.5 -(2.000/3.000)^0.5];
B=[5.000*(1.000-cos(2.000^0.500*i*h))/3^0.5;2.000*(2.000/3.000)^0.5*(-1.000+cos(5.000^0.5*i*h))]
xyz2=zeros(2,n1);
x12=zeros(2,1);
x12=A*B;
xyz2=x12;
for i=2:n1
    B=[5.000*(1.000-cos(2.000^0.500*i*h))/3^0.5;2.000*(2.000/3.000)^0.5*(-1.000+cos(5.000^0.5*i*h))];
    x12=A*B;
    xyz2(:,i)=x12;
end