Re: [問題] 請教如何跑出輸入狀態的圖形??已回收
試了一下 可以跑了 不過不知道正不正確
程式要跑滿久的
%--------------------------
function test
close all;
clear;
clc;
[t x] = ode45(@controller2,[0 9],[0.002,0,0]);
plot(t,x(:,3)) % x(:,3)就是ut
function dx = controller2(t,x)
d = 1;
k = 0.000378;
m = 0.0793;
R = 24;
L = 0.12829;
g = 9.8;
n = 0.004;
V = 0.5*(d^2+1)*x(1)^2+d*x(2)*x(1)+0.5*x(2)^2;
a = ((d^2+1)*x(1)+d*x(2))*x(2);
b = (d*x(1)+x(2))*-(k/m)/((x(1)+n)^2);
u =-(a+sqrt(a^2+b^4))/b;
%dot_V=((d^2+1)*x(1)+d*x(2))*x(2)+(d*x(1)+x(2))*-(k/m)/((x(1)+n)^2)*u;
ug=g*((x(1)+n)^2)/(k/m);
dx=zeros(3,1);
dx(1) = x(2);
dx(2) = -(k/m)*(u/((x(1)+n)^2));
dx(3) = ug+u;
%--------------------------
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.118.207.32
推
05/09 01:14, , 1F
05/09 01:14, 1F
→
05/09 01:22, , 2F
05/09 01:22, 2F
討論串 (同標題文章)
完整討論串 (本文為第 2 之 2 篇)：