怎么写出下面问题的matlab模型

如何用Matlab求解下面这个模型：在问题补充里：求大神指点，写出求解过程~

确定目标函数没有写错？
这还是凸优化问题么？

此微分方程组可用ode45（）函数求得其数值解。
首先，建立自定义函数
function dy = rigid(t,y)
dy = zeros(3,1);
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
end
其二，用ode45（）函数求出x、y值
[x,y] = ode45(@rigid,[0 12],[0 1 1]);
其三，用plot（）函数绘出x—y，x—dy/dx，x—d²y/dx²曲线图
plot(x,y(:,1),'-',x,y(:,2),'-.',x,y(:,3),'.')
运行结果

该问题实际上就是求非线性规划问题。可以用matlab的fmincon（）非线性规划函数来求解。求解思路：

1、建立目标函数 fmincon_fun( )，即

min z=5000*sum(xi)+6500*sum(yi)+200*sum(zi);

2、建立约束函数 fmincon_con( )，即

ceq1 = x1+y1-z2-3000

ceq2 = x2+y2+z2-z3-4500

ceq3 = x3+y3+z3-z4-3500

ceq4 = x4+y4+z4-z5-4000

ceq5 = x5+y5+z5-z6-4000

ceq6 = x6+y6+z6-5000

3、建立执行函数fmincon_main( )，即

x0=ones(1,18)*0.5  %初值

lb=zeros(1,18);   %下界值

ub=[ones(1,6)*3000 ones(1,6)*1500 ones(1,6)*5000]; %上界值

[x,fval,exitflag]=fmincon(@(x) fmincon_fun(x),x0,[],[],[],[],lb,ub,@(x) fmincon_con(x))

4、根据上述思路编程后，运行得到

最低生产成本F=最低生产成本F=129107500元

x1=3000件，x2=3000件，x3=3000件，x4=3000件，x5=3000件，x6=3000件

y1=1件，y2=1500件，y3=500件，y4=1000件，y5=1500件，y6=1500件

z1=1件，z2=1件，z3=1件，z4=1件，z5=1件，z6=500件

---

府谷县13937637066： 怎么用matlab程序表示下面的问题？
嬴旺安神： N=[3 9 5 6]; % 设备编号 V=[3 9; 9 5; 9 6]; % Excel的连接关系 % 连接矩阵,初始化为0,然后根据连接关系对有连接的元素赋值 M=zeros(length(N)); for i=1:size(V,1) r=find(V(i,1)==N); c=find(V(i,2)==N); M(r,c)=1; end disp(M)

府谷县13937637066： 用matlab如何编程解下面的问题 - ？
嬴旺安神： z=h*r; for i=2:n X(i)=X(i-1)+z*X(i-1)*(1-X(i-1)/K); end

府谷县13937637066： 怎么用matlab解决下面问题,要程序 - ？
嬴旺安神： 该问题可以这样来解决.1、用循环语句,求出当θ1=0到2π时的θ2、θ3、θ5、θ6的对应值.2、利用meshgrid()函数,对θ1和e数据平面网格化3、利用ω6的计算公式,计算对应θ1和e的ω6值.4、利用mesh()函数,绘出其θ1、e和ω6三维曲面图.

府谷县13937637066： (实验三)下面的题目用matlab怎么写?考虑正弦信号x(t)=sin(ω0t),若x(t)用频率固定在ωs =2π/T rad/s 采样,那么离散时间信号x[n]=x(nT)就是:x[n]= sin(ω0nT... - ？
嬴旺安神：[答案] n=0:8192;t=n*1/8192; X=sin(2*pi*1000.*t); subplot(2,1,1)stem(X(1:50))xlabel('n')ylabel('X[n]')subplot(2,1,2)plot(t(1:50),sin(2*pi*1000.*t(1:50)))xlabel('t')ylabel('x(t)')---------------------------一...

府谷县13937637066： 用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数和状态方程 - ？
嬴旺安神： (1)状态方程模型参数:编写matlab程序如下>> num=[1 7 24 24]; >> den=[1 10 35 50 24]; >> [A B C D]=tf2ss(num,den)得到结果:A=,B=,C=,D=[0] 所以模型为: =X+u,y=X(2) 零极点增益:编写程序 >> num=[1 7 24 24]; >> den=[1 10 35 50 24]; ...

府谷县13937637066： Matlab建立模型:狐狸与野兔问题 - ？
嬴旺安神： 首先用hu_tu.m建立微分方程组:function dX=hu_tu(t,X) % X=[x;y] dX=zeros(2,1); dX(1)=4*X(1)-0.02*X(1)*X(2); dX(2)=0.001*X(1)*X(2)-0.9*X(2); 其次,解这个微分方程(这里用ode45),建立run.m%%%%%%%%%%%%%%%%%%%%%%%...

府谷县13937637066： 用matlab怎么解决下面的问题:(要代码)最好有分析？
嬴旺安神： 试着做了一下,希望有所帮助. 代码 % max 2*x1+x2 % 约束条件: % x1>=0 % x2>=0 % -x1+2*x2<=2 % x1+x2<=4 % x1<=3 clc; clear all; f = -[2 1]; A = [-1 2 1 1]; b = [2; 4]; Aeq = []; beq = []; lb = [0; 0]; ub = [3; 4]; [x, fval] = linprog(f, A, b, Aeq, beq, lb, ub); x -fval 结果 Optimization terminated. x =3.0000 1.0000ans =7.0000 >>

府谷县13937637066： 如何用matlab建立数学模型及求解.哪位高手给个模版. - ？
嬴旺安神： 使用linprog函数.matlab中有详细的关于linprog的参数、用法的说明,可以查看.简单来说,linprog的一个常用标准形式是x = linprog(f,A,b,Aeq,beq,lb,ub),它用来求解最小化问题 min f(x), s.t. Ax ≤ b Aeq*x = beq lb ≤ x ≤ ub.所以针对你的问题,相当于minz'=-2x1+3x2,代码是x=linprog(-[2;3],[0,1;4,2;1,1],[12;20;6],[],[],[0;0],[]).运行得到结果x=[0;6],即x1=0,x2=6

府谷县13937637066： 如何用matlab画下面的函数图,条件在下面写出,想要步骤 - ？
嬴旺安神： w = 10; x = 0:0.5:10; % x = 0 0.5 1 1.5 ... 10 y = exp(-0.5*x).*sin(w*x); plot(x,y);

府谷县13937637066： 写出求下列表达式matlab程序 - ？
嬴旺安神： 在命令窗口中,输入>> z1=2*sind(85)/(1+exp(2)) %回车 z1 = 0.2375 %计算结果