jingquejie.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
t=(0:0.01:12);
aa=[1/3^0.5 0.5*(2/3)^0.5;1/3^0.5 -(2/3)^0.5];
n=length(t);
qx=zeros(3,n);
qx(1,:)=t;
qx(2:3,:)=aa*[5*(1-cos(2^0.5*t))/3^0.5;2*2^0.5*(-1+cos(5^0.5*t))/3^0.5];
save('jqj','qx');
plot(qx(1,:),qx(2,:),qx(1,:),qx(3,:))