1. 首頁 > 生活 > >

“matlab求線性規劃最大值“matlab如何求兩面相交線

生活百科matlab線性擬合

matlab求！速度
給出的題1是Van der Pol方程，x"(t) (x^2-1)x'(t) x=0，x(0)=0，x'(0)=0.25，可下列來求解 。
1、用ode23函解，其過程為
func=@(t,x)[x(2);-(x(1)-1)*x(2)-x(1)] %自微分方程函數
x0=[0,0.25]; %初始條件
tspan=[0,20];%時間區間
[t,x]=ode23(func,tspan,x0)
plot(t,x)%用圖像來表示其方程的數值解，即t—x曲線
2、用dsolve函數求解，其過程為
syms x
D2x=diff(x,2)%x"(t)—x(t)的二階導數
Dx=diff(x,1)%x'(t)—x(t)的一階導數
x=dsolve(D2x (x^2-1)*Dx x==0,x(0)==0,Dx==0.25)
運行該代碼，我們可以看到，得到警告信息（‘無法找到顯式解’）該方程沒有解析值 。




matlab求線性規劃最大值
matlab求線性規劃最大值可以按下列方法求解：
1、首先應根據問題，創建線性規劃極值目標，例如：
function maxZ=fmincon_fun(k);
maxZ=-（ 。。。。。。）; %目標函數表達式
end
2、第二據問題，創建線性規劃極值約束條件函數，例如：
function [g,h]=fmincon_con(k);
g= 。。。。。。; %約束條件不等式表達式
h= 。。。。。。; %約束條件等式表達式
end
3、確定k的初值，即k0=[ 。。。]
4、確定k的上下限，lb=[ 。。。],ub=[ 。。。]
5、使用fmincon函數求其極值，即
fun=@(k) fmincon_fun(k);
nonlcon=@(k) fmincon_con(k);
[k,fval,exitflag] = fmincon(fun,k0,A,b,Aeq,beq,lb,ub,nonlcon);
x1=k(1),x2=k(2), 。。。最優化值
maxZ=-fval最

求問各位大佬，這個matlab程序哪里錯了，為啥是個直線
題主給這個matlab程序運行結果為啥直線？，是錯了于題主給出Lorenz模狀態方程，我們知道該方程完全依賴于時間tspan條件，即tspan=0～t，當t比較小時，就會出現一條直線，但隨著t不斷增加，其圖形就會出現題主給的相位空間圖形（蝴蝶圖形） 。
分析題主的代碼，其主要錯誤是：
1、t值偏小，建議取 t=100 以上
2、rk4（）函數的自定義不太明確，一般應這樣定義（網上可以搜索） 。即
[t,x] =rk4@(t,x)func(t,x),x0,h,a,b)
這里，func是自定義微分方程組函數，x0是初始條件，a是初始時間，b是結束時間，h是步長
x0=[1,2,3]
a=0
b=100
【“matlab求線性規劃最大值“matlab如何求兩面相交線】h=0.01