24小时热门版块排行榜    

CyRhmU.jpeg
查看: 937  |  回复: 4

麦穗和稻谷

新虫 (初入文坛)

[交流] 自己的写的matlab跑不动,哪位大神帮我用牛顿迭代编一下已有4人参与


现在已知n=73(有73个题),π =0.6,a,b,c见下图 我要求解θ,我因为不会编写MATLAB,就把73组 abc值都代入公式,发现MATLAB跑不了,就是这个意思。
方程只有一个未知数θ, a,b,c,n, π都已知。能不能编程?

题号        a        b        c
A01        1.03        0.15        0.36
A02        0.87        -1.21        0.20
A03        0.48        -1.96        0.25
A04        0.79        -0.07        0.28
A05        0.47        -1.96        0.28
A06        0.52        -2.74        0.24
A07        0.41        -0.29        0.16
A08        0.29        -1.69        0.28
A09        0.57        -2.11        0.25
A10        0.61        -1.28        0.30
A11        0.63        -3.30        0.21
A12        0.52        -3.05        0.20
A13        0.41        -3.68        0.21
A14        0.44        0.18        0.25
A15        0.94        0.01        0.45
B01        1.18        -0.20        0.22
B02        1.27        -0.34        0.28
B03        1.39        -0.53        0.29
B04        1.12        -0.07        0.20
B05        0.73        0.02        0.15
B06        1.01        -0.05        0.22
B07        1.00        -0.50        0.16
B08        0.85        0.45        0.21
B09        1.02        -0.60        0.24
B10        0.96        -0.25        0.20
B11        0.95        0.47        0.10
B12        0.60        -0.86        0.13
B13        1.01        -0.32        0.10
B14        0.98        3.12        0.36
B15        0.57        -0.35        0.25
B16        0.86        -1.62        0.12
B17        0.87        0.47        0.14
B18        0.52        -2.40        0.23
B19        0.41        -1.63        0.16
B20        0.46        -1.51        0.16
B21        0.89        3.88        0.40
B22        1.16        -0.40        0.15
B23        0.55        -0.14        0.19
B24        0.84        -1.09        0.15
B25        0.92        -0.79        0.09
B26        0.89        -0.33        0.15
B27        0.77        -1.28        0.29
B28        1.16        1.09        0.11
B29        0.51        -0.76        0.17
B30        0.57        1.42        0.25
B31        0.63        -0.91        0.23
B32        0.93        2.10        0.21
B33        1.26        -0.91        0.08
B34        0.65        -2.98        0.17
B35        0.60        -0.94        0.10
B36        0.55        -1.65        0.19
B37        0.38        -0.17        0.16
B38        0.75        -2.11        0.16
B39        1.02        0.17        0.11
B40        0.39        -1.76        0.16
C01        1.03        0.71        0.50
C02        0.62        1.78        0.08
C03        0.49        2.91        0.22
C04        0.84        0.75        0.41
C05        1.05        4.39        0.01
C06        0.71        0.71        0.15
C07        1.03        0.26        0.30
C08        0.46        -1.32        0.21
C09        0.86        -0.11        0.21
C10        0.70        2.01        0.06
C11        0.91        0.76        0.13
C12        0.42        0.13        0.16
C13        0.90        1.80        0.15
C14        0.65        -0.96        0.09
C15        0.57        2.11        0.06
C16        0.68        0.91        0.13
C17        0.35        -0.31        0.18
C18        0.42        -1.10        0.20
自己的写的matlab跑不动,哪位大神帮我用牛顿迭代编一下
QQ图片20140916185605.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

仙木映月

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你查一下matlab如何求解隐式方程,会方便一点。至于方程能不能编程实现,那是一定可以的。
如果不会用隐式方程求解,还有一个办法,作图求解,你秩序要以theta做自变量作图就能找到交点,放大之后横坐标就是你要的解了。
Hardtosay.
2楼2014-09-17 08:45:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

解:
θ=1.03261
3楼2014-09-17 10:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gctpaper

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这个确实要编程,且要解方程。
而录入那些数据就很累了,耗时严重……
天外飞仙,自在悠然~
4楼2014-09-20 11:02:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suntree4152

铁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
将你的原始数字存入一个excel文档,然后数组a, b, c,可以通过matlab的xlsread命令读入,然后编写你图中所示的目标函数,记为obj(theta; a, b,c), 最后用matlab的solve命令来求解obj(theta; a, b, c)=0即可。
5楼2014-09-20 12:03:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 麦穗和稻谷 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见