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

登录注册

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

24小时热门版块排行榜

>虫友互识 (4286)
>留学DIY (2240)
>论坛更新日志 (2227)
>留学生活 (1780)
>海外博后 (1066)
>育儿交流 (840)
>海归之家 (700)
>医药类考试 (700)
>信息科学 (560)
>材料综合 (559)
>人文社科 (558)
>复合材料 (478)
>地学 (460)
>转基因 (420)
>我与小木虫的故事 (420)
>分子模拟 (420)

返回列表

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？

回复此楼

» 猜你喜欢

如何从铁电相到顺电相。已经有1人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有273人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有42人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有129人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复

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

如何利用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 的主题更新

返回列表

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

[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +24	医学老男孩 2026-08-13	50/2500	2026-08-15 06:38 by 学员8dgXkO
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4	k0dTPqJtl0jt 2026-08-14	5/250	2026-08-15 04:45 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +11	Tide man 2026-08-10	12/600	2026-08-15 02:12 by home3163
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	k0dTPqJtl0jt 2026-08-14	4/200	2026-08-15 01:28 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	HFw0lei2R37i 2026-08-14	5/250	2026-08-15 00:16 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4	HFw0lei2R37i 2026-08-14	8/400	2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7	启萌科技 2026-08-12	22/1100	2026-08-14 23:45 by Noways
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +5	Ldrop2023 2026-08-13	5/250	2026-08-14 18:48 by ssxclkj
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8	布布和一二 2026-08-10	11/550	2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候，自然基金真的不能太认真（我的申报经验） +6	majunge000 2026-08-11	8/400	2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23	Lanmanbaby 2026-08-09	37/1850	2026-08-14 13:40 by Equinoxhua
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10	家与远方 2026-08-10	15/750	2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9	majunge000 2026-08-10	12/600	2026-08-13 23:20 by iwuli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3	mbygzh 2026-08-08	4/200	2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8	gjjjzhong 2026-08-10	21/1050	2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了，巨变 +3	WH3796 2026-08-12	4/200	2026-08-13 14:13 by 小木虫6752397
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11	yufeiwaner 2026-08-09	13/650	2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6	瞬息宇宙 2026-08-10	7/350	2026-08-11 19:25 by Tide man
[基金申请] 面上项目filecode邪修 +5	西山十月 2026-08-09	7/350	2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读？ +10	Tide man 2026-08-08	18/900	2026-08-09 12:55 by 仁砚薪传

24小时热门版块排行榜

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

» 猜你喜欢

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

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