matlab 用solve解方程组如何只保留全部未知数都大于0的一组解_matlab吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：191,253贴子：
matlab 用solve解方程组如何只保留全部未知数都大于0的一组解收藏
比如 [V1,V2,V3,Vges]=solve('4=0.36+0.083932*Vges^2+0.^2 ','4=0.66+0.083932*Vges^2+0.016782*(Vges-V1)^2+0.^2 ','
4=0.96+0.083932*Vges^2+0.016782*(Vges-V1)^2+0.^2 ','Vges=V1+V2+V3','Real',true) 解出来这样的，但是只有第二组解才是符合要求全部大于0的，怎样才能使它只显示这一组解呢？？还有怎样保留4位小数呢？急求V1 =
-2.9181214
3.7274971 -2.7836328 -3.2889327 -3.7274971V2 =
-1.0599859
1.0599859 -2.5221083
2.2311552 -2.2892777 -2.2311552
2.2892777V3 =
-1.4836024
-2.221335 -2.3329928
2.221335 -2.6913728
2.3329928Vges =
-5.4617096
3.793453 -1.8947734
-3.793453 -2.4581971
企业版MATLAB中文官方试用下载,让工程师和科学家更高效地进行科学计算.
您好，您最后是怎么解决的呢？
楼主，还有代码吗？可以给我发一份吗
楼主，还有代码吗？可以给我发一份吗
登录百度帐号matlab解含变量方程组问题_百度知道
matlab解含变量方程组问题
求matlab高手解答：
要计算一个坐标轴旋转变换，使坐标系中向量A2=[a2x, a2y, a2z];变成A1=[a1x,0,0];（a2x^2+a2y^2+a2z^2=a1x^2）
首先设置变量
写出坐标轴旋转矩阵（先绕y轴旋转b角度，再绕z...
0, a2z].5;A1=[(a2x^2+a2y^2+a2z^2)^0换成A2=[a2x, a2y
我有更好的答案
看你的程序中，有7个未知量，用5个未知量表示2个未知量，当然不好解，且有三角函数式当算式中有三角函数时，往往solve不出结果，更是这样，解是很不确定的。因此，特别是未知数有2个以上，要尽量减少未知量的个数
如果我只需要在[0,pi)之间求解是否可以呢？
采纳率：77%
来自团队：
为您推荐：
其他类似问题
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。matlab中等式两边都有未知量,怎么求出这个未知量？_百度知道
matlab中等式两边都有未知量,怎么求出这个未知量？
我有更好的答案
f1(x)=f2(x)问题；转化为两条曲线y=f1(x)和y=f2(x)的交点。只需要在一定范围内遍取x画图就可以了。也可以通过判据abs(f1(x)-f2(x))&eps 判断
采纳率：87%
为您推荐：
其他类似问题
雪地胎的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。当前位置：
&关于matlab求解非线性指数方程组出现问题
关于matlab求解非线性指数方程组出现问题
作者 blbailei
我现在用Matlab 求解如下非线性指数方程组：
& && && &该非线性指数方程组，是为了求解极大对数似然估计，利用diff()对各变量取偏导数求得。
& & function f=nle(xvert)
& & b11=xvert(1);
& & b12=xvert(2);
& & b13=xvert(3);
& & b21=xvert(4);
& & b22=xvert(5);
& & b23=xvert(6);
& & b31=xvert(7);
& & b32=xvert(8);
& & b33=xvert(9);
& & f(1)= 2/(b11 + b12/2 + b13/2) - 2/(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2) - (28*exp(-7*b12)*exp(-7*b13)*exp(-14*b11))/(exp(-7*b12)*exp(-7*b13)*exp(-14*b11) - exp(-7*b22)*exp(-7*b23)*exp(-14*b21)) - 28 ;
& & f(2)= 2/(b21 + b22/2 + b23/2) - (28*exp(- 14*b21 - 7*b22 - 7*b23))/(exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33)) - (2*(((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)^2 - (exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) + (14*exp(- 14*b21 - 7*b22 - 7*b23)*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)))/(((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) + 1) - 2/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) + (2*(((exp(- 14*b11 - 7*b12 - 7*b13) - exp(- 14*b21 - 7*b22 - 7*b23))*(b11 + b12/2 + b13/2))/(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2)^2 + (14*exp(- 14*b21 - 7*b22 - 7*b23)*(b11 + b12/2 + b13/2))/(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2))*(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2))/((exp(- 14*b11 - 7*b12 - 7*b13) - exp(- 14*b21 - 7*b22 - 7*b23))*(b11 + b12/2 + b13/2)) - 28 ;
& & f(3)= 1/(b11 + b12/2 + b13/2) - 1/(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2) - (14*exp(-7*b12)*exp(-7*b13)*exp(-14*b11))/(exp(-7*b12)*exp(-7*b13)*exp(-14*b11) - exp(-7*b22)*exp(-7*b23)*exp(-14*b21)) - 14 ;
& & f(4)= 1/(b21 + b22/2 + b23/2) - (14*exp(- 14*b21 - 7*b22 - 7*b23))/(exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33)) - (2*(((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(2*(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)^2) - (exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))/(2*(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)) + (7*exp(- 14*b21 - 7*b22 - 7*b23)*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)))/(((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) + 1) - 1/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) + (2*(((exp(- 14*b11 - 7*b12 - 7*b13) - exp(- 14*b21 - 7*b22 - 7*b23))*(b11 + b12/2 + b13/2))/(2*(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2)^2) + (7*exp(- 14*b21 - 7*b22 - 7*b23)*(b11 + b12/2 + b13/2))/(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2))*(b11 + b12/2 + b13/2 - b21 - b22/2 - b23/2))/((exp(- 14*b11 - 7*b12 - 7*b13) - exp(- 14*b21 - 7*b22 - 7*b23))*(b11 + b12/2 + b13/2)) - 14 ;
& & f(5)= (28*exp(- 14*b31 - 7*b32 - 7*b33))/(exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33)) + (2*(((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)^2 + (14*exp(- 14*b31 - 7*b32 - 7*b33)*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)))/(((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) + 1) + 2/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2) - 28 ;
& & f(6)= ((28*exp(- 14*b31 - 7*b32 - 7*b33)*(b21 + b22/2 + b23/2))/(2*b21 + b22 + b23 - 2*b31 - b32 - b33) + ((exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(b21 + b22/2 + b23/2 - b31 - b32/2 - b33/2)^2)/((2*(exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33))*(b21 + b22/2 + b23/2))/(2*b21 + b22 + b23 - 2*b31 - b32 - b33) + 1) + 2/(2*b21 + b22 + b23 - 2*b31 - b32 - b33) + (14*exp(- 14*b31 - 7*b32 - 7*b33))/(exp(- 14*b21 - 7*b22 - 7*b23) - exp(- 14*b31 - 7*b32 - 7*b33)) - 14 ;
主函数代码：
& & x0=[0.01;0.002;0.003;0.004;0.005;0.006;0.007;0.008;0.009];
& & [x,FVAL]=fsolve(@nle3,x0,optimset( 'PlotFcns', @optimplotfval));
使用fsolve 函数时，显示：
Warning: Trust-region-dogleg algorithm of FSOLVE cannot handle non- using Levenberg-Marquardt algorithm instead.
& In fsolve at 309
&&In main at 10
No solution found.
fsolve stopped because the last step was ineffective. However, the vector of function
values is not near zero, as measured by the default value of the function tolerance.
&stopping criteria details&
问题是，如何使用设置optimset，来使用 Levenberg-Marquardt algorithm 呢？
用1stOpt跑了下，似乎的确无解。好好检查下方程，尤其是最后三个。
我建议你对原始的exp模型进行转换，将其转化成对数模型进行求解，我估计最大的问题在于数据计算中产生的误差，
24小时热帖
下载小木虫APP
与700万科研达人随时交流扫二维码下载作业帮
3亿+用户的选择
下载作业帮安装包
扫二维码下载作业帮
3亿+用户的选择
matlab怎么解物理中含有参数的方程组?如要用 M-T*R=k1*b1; T*r=k2*b2; b1*R=r*b2 算出b1 b2,其中b1 b2 T是未知数,结果用已知量表示.这题用笔一下就可以算出 我只是举个例子 用solve解这方程组怎么含有未知数呀 肿么回事
作业帮用户
扫二维码下载作业帮
3亿+用户的选择
syms M T R r k1 k2 b1 b2f = solve(M-T*R-k1*b1,T*r-k2*b2,b1*R-r*b2,b1,b2,T);
为您推荐：
其他类似问题
扫描下载二维码}