[問題] 解非線性問題?已回收
以下是我的程式，遇到的問題是[-15,1,1]是x的初值代入時跑不出初
始狀態圖形??
怎麼會這樣??
如果在[0,0,0]時也是沒有圖形的??
可以跟我說是怎麼回事嗎??
該怎麼解決呢?
幫幫忙!! 謝謝
function new
clc;
clear;
close all;
[t x] = ode45(@controller2,[0 1],[-1.5,1,1])
plot(t,x);
function dx = controller2(t,x)
k = 0.000378;
m = 0.0793;
R = 1;
L = 128.29;
g = 9.8;
n = 0.005;
d = 5;
u1 = (0.5*(d^2+1)*x(1)+0.5*d*x(2)*x(3)+x(1)*x(2)^2+x(1)*x(3)^2)*x(2);
u2 = (0.5*d*x(1)*x(3)+0.5*x(2)+x(2)*x(3)^2+x(1)*x(2)^2)*
-(k/m)*(x(3)^2/(x(1)+n)^2);
u3 = (0.5*d*x(1)*x(2)+0.5*x(3)+x(3)*x(2)^2+x(3)*x(1)^2)*-(R/L)* x(3);
b = (0.5*d*x(1)*x(2)+0.5*x(3)+x(3)*x(2)^2+x(3)*x(1)^2)*(1/L);
a = u1+u2+u3;
u = -(a+sqrt(a^2+b^4))/b;
dx=zeros(3,1);
dx(1) = x(2);
dx(2) = -(k/m)*(x(3)^2/(x(1)+n)^2);
dx(3) = -(R/L)* x(3)+(1/L)*u;
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.121.136.17
※ 編輯: agersir 來自: 140.121.136.17 (02/23 18:45)
※ 編輯: agersir 來自: 220.136.20.17 (02/23 23:32)