mainf.m

内燃机转子仿真

%mainf
clear all;
clc;
close all;
m=[1.000 1.000 0.000;1.000 2.000 0.000;0.000 0.000 3.000];%3
c=[1.000 0.000 0.000;0.000 0.000 1.000;0.000 1.000 0.000];%3
k=[2.000 1.000 0.000;1.000 3.000 0.000;0.000 0.000 4.000];%3
F=[2.000;4.000;6.000];%3
n=length(F);
[M,C,K]=arryf(m,c,k,n);
x1=[-1;-3;-2];%3
x2=[1;2;3];%3
b=1.4;
tspan=6;
S=inv(M)*K;
[px pn]=eig(S);
pn=pn^0.5;
p=diag(pn);
w=1;
fmin=min(p);
if w==0
Tmax=1/2/pi/fmin;
else
Tmax=2*pi/w;
end    
h=Tmax/320.0000;
dn=rem(tspan,h);
dnn=tspan-dn;
n1=tspan/h;
f0=zeros(n,1);
xy1=wilson1((M,C,K,dK,f,x1,x2,n,b,n1,w);
xy2=wilson2(M,C,K,dK,f,x1,x2,b,n,n1,w);
figure(1);
plot(x,y);
figure(2);
plot(x,y);
figure(3);
plot(x,y);
pyy=gonglp(t,yy,n1);