“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