exe4_21.m

来自「MATLAB 简单编程 包括很基础的程序」· M 代码 · 共 19 行

function y=exe4_21;
N=10;
g=inline('x.^2./y-x.*cos(y)');
[x,y]=ode23(g,0,5,1);
figure(1);plot(x,y);hold on;
xi=0:0.1:5;
a1=polyfit(x,y,N);
yi=polyval(a1,xi);
plot(xi,yi,'r-');
l=length(a1);
s=[num2str(a1(1)) '*x.^' num2str(l-1) ] ;
for i=2:l-1
    s=[s '+(' num2str(a1(i)) ')*x.^' num2str(l-i) ];
end
s=[s '+(' num2str(a1(l)) ')'];
F=inline(s);
Q=quad(F,0,5);
Q