datachuli.m

内燃机转子仿真

clear all;
clc;
close all;
format long;
load 'sjjl';
pinglvb=(0.00:0.01:2);
ni=length(pinglvb);
mpyy=zeros(2,ni);
for i=1:ni
    w=i*0.01*fn;
    pim1x(w);
    mpyy1=gonglp1;
    
    mpyy(:,i)=[w;mpyy1];
end
x(1,:)=mpyy(1,:);
y(1,:)=mpyy(2,:);
c=polyfit(x,y,5);
cc=c(1)*x.^5+c(2)*x.^4+c(3)*x.^3+c(4)*x.^2+c(5)*x+c(6);
figure(1);
plot(x,cc,x,y,'-.k');
figure(2);
plot(x,y);