24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

study3040

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 31.3
红花: 1
帖子: 34
在线: 35.6小时
虫号: 493782
注册: 2008-01-10
性别: GG
专业: 机构学与机器人

[求助] 亲们，求助一下，这个非线性方程组怎么解？已有1人参与

各位亲们，想解这个非线性方程组，含有积分项，且上下界中含有变量。
求助各位了，怎么用matlab解呢？结果是多少，求解求解~~非常3ks
或者有其他解的办法没，都可以

求解非线性方程组.jpg

回复此楼

» 猜你喜欢

最近几年招的学生写论文不引自己组发的文章已经有5人回复
职称评审没过，求安慰已经有54人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

一个关于有约束非线性方程组的求解问题已经有7人回复
求助，帮忙求解下面方程组，谢谢！已经有13人回复
关于非线性方程组的问题已经有9人回复
求助解复杂非线性方程组的好的方法已经有24人回复
用matlab求解一个非线性方程组的解已经有4人回复
请教一个matlab求解非线性方程组的问题已经有9人回复
用fortran程序遗传算法解非线性方程组已经有7人回复
非线性方程求解问题，求助各位兄弟姐妹们~~ 已经有7人回复
MATLAB求解非线性方程组已经有5人回复
【求助】非线性方程组的求解问题已经有6人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
求助matlab---fsolve解非线性方程组已经有6人回复
matlab求解非线性方程组已经有16人回复
请教Matlab高手解多元非线性方程组中的参数已经有6人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】求助解决这样一个非线性方程组已经有8人回复
【求助】非线性方程组求解已经有3人回复
【求助】求助，解非线性方程组，急已经有5人回复
【求助】求教matlab解非线性方程组已经有9人回复

被上帝咬过的苹果

1楼 2014-06-20 12:31:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

study3040

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 31.3
红花: 1
帖子: 34
在线: 35.6小时
虫号: 493782
注册: 2008-01-10
性别: GG
专业: 机构学与机器人

引用回帖:

2楼: Originally posted by dingd at 2014-06-20 14:13:29
楼主图方便自己贴个图片省事，难道还要帮助的人照着图片一个个输入？应先学会懂得如何求助。

小弟不常发贴，不懂规矩，请见谅哈！程序如下。
求解时会出出现如下问题：
Warning: Infinite or Not-a-Number function value encountered.
查了下，有人说是存在第二类间断点，不知何解。
1、不知我写的这段程序是否可行？
2、是否有其他别的方法求解含积分的方程组？

function F = myfun2(x)
syms p theta theta1 q g alfa;
s0=0.055;    %%已知参数1，0.055
r=0.04;       %%已知参数2，0.04
theta1=x(1); %%待求量1初值，初值范围0<theta1<1.57
alfa=x(2);    %%待求量2赋初值，初值范围0<alfa<1.57
p=x(3);       %%待求量3赋初值，初值范围-10<p<10
q=x(4);       %%待求量4赋初值，初值范围-10<q<10
c=1000;       %%常数

g=-p*cos(alfa)-q*sin(alfa); %计算中间变量g
fun_cal1=@(aaa)1./sqrt(abs(2*(p*cos(aaa)+q*sin(aaa)+g)/c)); %积分函数1
fun_cal2=@(bbb)cos(bbb)./sqrt(abs(2*(p*cos(bbb)+q*sin(bbb)+g)/c)); %积分函数2
fun_cal3=@(ccc)sin(ccc)./sqrt(abs(2*(p*cos(ccc)+q*sin(ccc)+g)/c)); %积分函数3

%%待求方程组如下：（4个方程，想求4个变量theta1,alfa,p,q）
F = [s0-theta1*r-quad(fun_cal1,theta1,alfa)+quad(fun_cal1,alfa,pi/2);
r*(1-sin(theta1))-quad(fun_cal2,theta1,alfa)+quad(fun_cal2,alfa,pi/2);
r*cos(theta1)-quad(fun_cal3,theta1,alfa)+quad(fun_cal3,alfa,pi/2);
r-1/sqrt(abs(2*(p*cos(theta1)+q*sin(theta1)+g)/c))];

% options=optimset('Display','iter');
% fsolve(@myfun2,[0.4,1.2,1.2,1.3]) %%初值