| 查看: 2291 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
亲们，求助一下，这个非线性方程组怎么解？ 已有1人参与
|
||
|
各位亲们，想解这个非线性方程组，含有积分项，且上下界中含有变量。 求助各位了，怎么用matlab解呢？结果是多少，求解求解~~非常3ks 或者有其他解的办法没，都可以  求解非线性方程组.jpg |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个关于有约束非线性方程组的求解问题
已经有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人回复
arising2010
铜虫 (正式写手)
- 应助: 32 (小学生)
- 金币: 1433.4
- 红花: 11
- 帖子: 996
- 在线: 412.6小时
- 虫号: 2464689
- 注册: 2013-05-15
- 专业: 可再生与替代能源利用中的
3楼2014-06-20 15:49:57
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼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]) %%初值 |
4楼2014-06-20 16:48:40
5楼2014-06-20 16:50:24
40