版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

notdingy

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 3968.2
散金: 27
帖子: 681
在线: 31.8小时
虫号: 3429502
注册: 2014-09-20
性别: MM
专业: 光学

[求助] 利用MATLAB计算如下方程组，求大神算出结果。解析法显示为空，可能要利用数值法已有3人参与

L1='0.8334=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.14-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.14-1/z)^2)))*u)^2-(2*pi*y*(1/1544.14-1/z)^2))';
L2='0.8799=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.145-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.145-1/z)^2)))*u)^2-(2*pi*y*(1/1544.145-1/z)^2))';
L3='0.91092=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.15-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.15-1/z)^2)))*u)^2-(2*pi*y*(1/1544.15-1/z)^2))';
L4='0.92915=(x^2*(sinh(sqrt(x^2-(2*pi*y*(1/1544.155-1/z)^2)))*u)^2)/(x^2*(cosh(sqrt(x^2-(2*pi*y*(1/1544.155-1/z)^2)))*u)^2-(2*pi*y*(1/1544.155-1/z)^2))';
solve(L1,L2,L3,L4,x,y,z,u)

x，y，z，u为未知数，四个方程形式都一样。利用fsolve？

回复此楼

» 猜你喜欢

基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有100人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有13人回复
散金币祈福已经有94人回复
青基已中已经有1人回复
散金币祈福已经有66人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复

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

如何利用MATLAB模拟建筑物，用射线跟踪法计算接收机状态已经有7人回复
怎么利用matlab计算如下函数的积分已经有8人回复
matlab计算数据的邻接矩阵已经有0人回复
MATLAB计算问题已经有0人回复
matlab计算晶相含量已经有3人回复
求一份利用蒙特卡罗方法，并结合MATLAB计算光在胶体中散射系数的代码已经有1人回复
求教如何利用matlab 模拟计算两块异性磁体之间磁感应强度分布？已经有0人回复
如何利用matlab 模拟计算两块异性磁体之间磁感应强度分布？已经有0人回复
利用MATLAB计算热管轴向参数的问题已经有1人回复
怎么利用matlab计算出非线性函数模型的系数？急求！！！！！已经有4人回复
利用Matlab模拟空间滤波实验已经有43人回复
计算数据，利用matlab画三维曲面图已经有2人回复
如何利用Matlab编程计算热加工图中功率耗散和失稳值（利用Origin画图）已经有0人回复
【求助】如何利用MATLAB计算N元N次方程已经有9人回复
【求助】如何利用matlab计算这个无穷级数的和已经有12人回复

1楼 2015-11-17 19:55:08

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
notdingy: 金币+2, ★★★很有帮助 2015-11-18 17:39:05

数值解的话用1stOpt试试，很简单的。

[ 发自手机版 http://muchong.com/3g ]

赞一下(1人)

回复此楼

2楼2015-11-17 22:51:58

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

凌风7388

铁虫 (正式写手)

应助: 4 (幼儿园)
金币: 1153.2
散金: 72
红花: 5
帖子: 573
在线: 160.6小时
虫号: 759985
注册: 2009-04-29
专业: 化学反应工程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
notdingy: 金币+8, ★★★很有帮助 2015-11-18 17:38:48

用fsolve可以解，但是解之间需要给定初值，初值不同得到的结果就不同。
比如假设x=y=z=u=1 1 1 1，得到的结果依次为：1.8634 0.3599 1.9406 1.1542
比如假设x=z=1，y=u=2，得到的结果依次为：1.6796 2.0456 0.8471 2.4493

赞一下(1人)

回复此楼

3楼2015-11-18 09:36:02

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

notdingy

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 3968.2
散金: 27
帖子: 681
在线: 31.8小时
虫号: 3429502
注册: 2014-09-20
性别: MM
专业: 光学

引用回帖:

3楼: Originally posted by 凌风7388 at 2015-11-18 09:36:02
用fsolve可以解，但是解之间需要给定初值，初值不同得到的结果就不同。
比如假设x=y=z=u=1 1 1 1，得到的结果依次为：1.8634 0.3599 1.9406 1.1542
比如假设x=z=1，y=u=2，得到的结果依次为：1.6796 ...

具体的程序要怎么写呢，没用过fsolve

赞一下

回复此楼

4楼2015-11-18 10:39:19

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

notdingy

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 3968.2
散金: 27
帖子: 681
在线: 31.8小时
虫号: 3429502
注册: 2014-09-20
性别: MM
专业: 光学

引用回帖:

2楼: Originally posted by dingd at 2015-11-17 22:51:58
数值解的话用1stOpt试试，很简单的。

网上找了很多1stopt，都安装不成功

赞一下

回复此楼

5楼2015-11-18 10:40:19

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

引用回帖:

5楼: Originally posted by notdingy at 2015-11-18 10:40:19
网上找了很多1stopt，都安装不成功...

坛子里就有啊

回复此楼

6楼2015-11-18 11:44:10

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

1314168apple

金虫 (知名作家)

应助: 68 (初中生)
金币: 677
红花: 12
帖子: 6872
在线: 1462.3小时
虫号: 287760
注册: 2006-10-21
专业: 色谱分析

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
notdingy: 金币+10, ★★★很有帮助 2015-11-18 17:38:30

fun=@(x)[0.8334-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));0.8799-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.145-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.145-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.145-1/x(3))^2));0.91092-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.15-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.15-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.15-1/x(3))^2));0.92915-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.155-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.155-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.155-1/x(3))^2))];

TT=fsolve(fun,randn(4,1));

都没有精确的解，可能不存在？检查检查。
误差(好些解得到的差不多)：

fun(TT)

-0.0549
-0.0084
0.0226
0.0408

赞一下

回复此楼

为了最终理解你所不理解的，你必须经历一条愚昧无知的道路。为了占有你从未占有的东西，你必须经历被剥夺的道路。为了达到你现在所不在的名位，你必须经历那...

7楼2015-11-18 12:19:55

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

凌风7388

铁虫 (正式写手)

应助: 4 (幼儿园)
金币: 1153.2
散金: 72
红花: 5
帖子: 573
在线: 160.6小时
虫号: 759985
注册: 2009-04-29
专业: 化学反应工程

引用回帖:

4楼: Originally posted by notdingy at 2015-11-18 10:39:19
具体的程序要怎么写呢，没用过fsolve...

function xFsolve
clear all;clc
x0=[1 1 1 1]';x1=fsolve(@nonlineqs,x0)
x0=[1 2 1 2]';x2=fsolve(@nonlineqs,x0)

function f=nonlineqs(x)
f(1)=0.8334-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));
f(2)=0.8799-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));
f(3)=0.91092-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));
f(4)=0.92915-(x(1)^2*(sinh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2)/(x(1)^2*(cosh(sqrt(x(1)^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2)))*x(4))^2-(2*pi*x(2)*(1/1544.14-1/x(3))^2));

赞一下

回复此楼

8楼2015-11-18 13:28:44

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

相关版块跳转我要订阅楼主 notdingy 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5	ASdOkHsho7FD 2026-08-28	11/550	2026-08-29 20:18 by 63GmZBFV29mG
[基金申请] 投票: 有多少人是今天查系统知道结果的？ +16	爱看书的可乐 2026-08-26	18/900	2026-08-29 20:07 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15	晴天加油 2026-08-26	16/800	2026-08-29 18:28 by symmetry
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3	ASdOkHsho7FD 2026-08-28	4/200	2026-08-29 14:13 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4	ASdOkHsho7FD 2026-08-28	6/300	2026-08-29 11:40 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗？ +3	349506619 2026-08-28	3/150	2026-08-29 10:49 by 木风9012
[基金申请] 2026年叶企孙基金 +4	bud_bud 2026-08-27	7/350	2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12	Kittylucky 2026-08-27	13/650	2026-08-29 00:04 by superceng
[基金申请] 基金未中，这种答复是模板吗？ +5	zhaosm1982 2026-08-27	6/300	2026-08-27 16:00 by lfy8008
[基金申请] 为什么国际(地区)合作与交流项目没有放榜？ 10+3	majunge000 2026-08-26	11/550	2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15	Edward_pc 2026-08-26	23/1150	2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23	koalala 2026-08-24	26/1300	2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9	trojank 2026-08-26	9/450	2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4	bjdxyxy 2026-08-26	4/200	2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18	Ldrop2023 2026-08-26	18/900	2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4	yanglien 2026-08-26	5/250	2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3	paperzjh 2026-08-26	3/150	2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19	angus9576 2026-08-25	23/1150	2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6	ziyangfang 2026-08-25	9/450	2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8	scalable 2026-08-24	8/400	2026-08-25 12:52 by jnhyjjm

24小时热门版块排行榜

[求助] 利用MATLAB计算如下方程组，求大神算出结果。解析法显示为空，可能要利用数值法 已有3人参与

» 猜你喜欢

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

[求助] 利用MATLAB计算如下方程组，求大神算出结果。解析法显示为空，可能要利用数值法已有3人参与