查看: 311|回复: 3|关注: 0

[已解决] 请问带三角函数的多元非线性方程组有什么求解的方法

[复制链接]

新手

16 麦片

财富积分


050


6

主题

13

帖子

0

最佳答案
我用的是fsolve函数,但是初值选取一直得不到解答,请问这类方程还有其他的求解方法吗
代码如下:
H=6;                %%墙高
gama=18;            %%土体容重
c=10;               %%填土黏聚力
n1=tand(20);        %%填土内摩擦角
q0=10;              %%表面外加荷载
a=deg2rad(70);      %%墙背倾角
delta=deg2rad(10);  %%墙土摩擦角
b=deg2rad(10);      %%填土表面倾角
f=@(x)[x(1)*sin(a-delta)*0.7579+x(3)*sin(x(4)-b)-9.3913;
x(1)*cos(a-delta)*0.7579+x(3)*cos(x(4)-b)-52.9281-58.0033;
x(2)*sin(a-delta)*0.8257-x(3)*sin(x(4)-b)-14.1916+30.9309*sin(0.0935-b)*4.6582;
x(2)*cos(a-delta)*0.8257-x(3)*cos(x(4)-b)-53.4482+30.9309*cos(0.0935-b)*4.6582]
[x,fval,exitflag]=fsolve(f,[17.9777,42.5202,38.1228,0])


回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

5

主题

2338

帖子

706

最佳答案
  • 关注者: 196
发表于 2020-6-7 07:42:46 | 显示全部楼层 |此回复为最佳答案
可以试着把f改写成一个向量的模(直接在原有基础上加norm函数)然后对他进行最小值优化,如果方程有解,这个模就是最小值零,如果没法满足全部方程,只能逼近零。
这样可以求得一个近似解
   16.1583   34.7126  104.6563    0.1638
提问请:①准确描述问题②提出你的思考(等着抄作业的一律锁帖)③提供代码文本而非截图④及时反馈
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


6

主题

13

帖子

0

最佳答案
 楼主| 发表于 2020-6-7 09:53:03 | 显示全部楼层
TouAkira 发表于 2020-6-7 07:42
可以试着把f改写成一个向量的模(直接在原有基础上加norm函数)然后对他进行最小值优化,如果方程有解,这 ...

请问能给下具体代码吗,这个问题弄了好久没做出来,谢谢大佬
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

中级

1197 麦片

财富积分


5001500


0

主题

2762

帖子

257

最佳答案
  • 关注者: 171
发表于 2020-6-23 10:28:24 | 显示全部楼层
flyer63 发表于 2020-6-7 09:53
请问能给下具体代码吗,这个问题弄了好久没做出来,谢谢大佬

楼主的方程组应该没有精确的实数解,近似解不加范围约束的话也有好多:
1:
x1: 4.05107617024465
x2: 45.8256699384377
x3: 109.451256379942
x4: 0.236916132191357

2:
x1: 3.61186887601818
x2: 46.2288129787646
x3: 109.635693085524
x4: 0.239445798783332

3:
x1: 0.00446044746381751
x2: 49.5400093550759
x3: 111.176692236181
x4: 0.259902041872311
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表