| 查看: 1579 | 回复: 7 | ||
notdingy铁虫 (正式写手)
|
[求助]
利用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？ |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有64人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2015-11-17 22:51:58
凌风7388
铁虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1153.2
- 散金: 72
- 红花: 5
- 帖子: 573
- 在线: 160.6小时
- 虫号: 759985
- 注册: 2009-04-29
- 专业: 化学反应工程
3楼2015-11-18 09:36:02
notdingy
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3963.2
- 散金: 27
- 帖子: 680
- 在线: 31.8小时
- 虫号: 3429502
- 注册: 2014-09-20
- 性别: MM
- 专业: 光学
4楼2015-11-18 10:39:19
notdingy
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3963.2
- 散金: 27
- 帖子: 680
- 在线: 31.8小时
- 虫号: 3429502
- 注册: 2014-09-20
- 性别: MM
- 专业: 光学
5楼2015-11-18 10:40:19
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2015-11-18 11:44:10
1314168apple
金虫 (知名作家)
- 应助: 68 (初中生)
- 金币: 677
- 红花: 12
- 帖子: 6872
- 在线: 1462.3小时
- 虫号: 287760
- 注册: 2006-10-21
- 专业: 色谱分析
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
notdingy: 金币+10, ★★★很有帮助 2015-11-18 17:38:30
感谢参与，应助指数 +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
凌风7388
铁虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1153.2
- 散金: 72
- 红花: 5
- 帖子: 573
- 在线: 160.6小时
- 虫号: 759985
- 注册: 2009-04-29
- 专业: 化学反应工程
|
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
