24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lzht1986

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2462.1
散金: 50
红花: 2
帖子: 93
在线: 60.9小时
虫号: 2592348
注册: 2013-08-11
性别: GG
专业: 机械摩擦学与表面技术

[求助] 求助一个简单的方程，用MATLAB求解已有2人参与

方程为图片中的方程式
R为常量
S为已知量，且S为一组数据
求不同S值时对应的X值
我用MATLAB编码总是提示我有问题，还请高人帮我看看，谢谢啦

1R)HI@X$%IP{[Q4I%90`K7M.png

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有234人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

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

关于一个看着挺简单的非线性方程组，被搞得晕乎乎。。。。 T-T 特来求助大神们已经有5人回复
COMSOL新手，关于COMSOL的简单程序求助，求大神随手帮忙已经有3人回复
求助很简单的Lyapunov方程的LMI解法与过程已经有9人回复
求助一个基础二价铁，三价铁浓度比计算的问题已经有0人回复
comsol菜鸟求助，一个简单的偏微分方程，不知道边界条件如何设置已经有2人回复
金币全部求助，求二次型，矩阵，不等式化简方面的高手，只是二阶方程跟三阶矩阵已经有2人回复
求助：matlab求解一个方程组，其中一个方程有三角函数已经有10人回复
向各位求助：一个简单的仿真错误已经有19人回复
求助matlab解非线性方程组并作图的问题，需要具体程序已经有1人回复
求助大神：要做参数回归，fortran和matlab，该选择哪个？已经有14人回复
急切求助！一个mathematica求解一个简单一元非线性方程的问题已经有4人回复
Mathematica 简单四元一次方程组，求解问题~ 已经有9人回复
求助，一个很简单的二元一次微分方程组的解析解已经有2人回复
一个简单的问题，关于微分方程，泣血求助！已经有1人回复
【求助】末学愚钝，请教一个Poisson方程的问题已经有2人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】关于数据拟合，请教数学达人已经有4人回复
【求助】急求，方程的求助已经有5人回复
急求，求助方程的计算已经有5人回复
【讨论】【求助】一个微分方程的解已经有7人回复

1楼 2016-04-19 14:33:10

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

8楼: Originally posted by lzht1986 at 2016-04-20 15:32:54
大哥，谢谢你的指导。数据我带进去了代码可行，但是出来的数据有点问题，能否麻烦您再帮我看看，谢谢啦！！！为什么中间有一多半的数据解出来是0.100000呢？
ans =

Columns 1 through 2

0.0020892490 ...

因为初值设置就是0.1，有些方程用0.1的初值没有得到解，把初值修改成其他数值。

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

lzht1986

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

9楼2016-04-20 15:35:57

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

查看全部 16 个回答

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lzht1986: 金币+30, ★★★很有帮助, 谢谢了，但是提示还是不行啊 2016-04-20 14:13:26

CODE:

function jie_fangcheng

clear all;clc

format long

global s

s_data=[0.1 0.2 0.3 0.4];

for i=1:length(s_data)

    s=s_data(i);

    [x(i) fval(i)]=fsolve(@obj,0.1);

end

[s_data' x' fval']

function f=obj(x)

global s

R=1;

f=R^2*asin(x/R)-x*sqrt(R^2-x^2)-s;

R值自行设定。
上述代码以S=[0.1 0.2 0.3 0.4]为例。
计算结果：
S x 方程残差
0.100000000000000 0.515718393818021 0.000000000076744
0.200000000000000 0.637417584593469 0.000000003345383
0.300000000000000 0.717211724786695 0.000000001460760
0.400000000000000 0.776601502032687 0.000000376769810

赞一下

回复此楼