主程式如此,報錯的就是(x:,i+1)那一句,不知道怎么改,
function [x,t]=LONGE1(x0,t0,a,b,h)
n=(b-a)/h;
x(1)=x0;
t(1)=t0;
for i=1:n
K1=fun(x(i),t(i));
K2=fun(x(i)+h/2*K1,t(i)+h/2);
K3=fun(x(i)+h/2*K2,t(i)+h/2);
K4=fun(x(i)+h*K3,t(i)+h);
x(:,i+1)=x(:,i)+h*(K1+2*K2+2*K3+K4)/6;
t(i+1)=t(i)+h;
end
end
函式子程式如下
function z=fun(x,t)
ss=[-10 -100 0 0;1 0 0 0;0 0 -5 -25;0 0 1 0];
F=10;
z=ss*x+F*sin(10*t);
就是函式里多了個矩陣就無法執行了
轉載請註明出處,本文鏈接:https://www.uj5u.com/qukuanlian/35120.html
標籤:區塊鏈技術
