作者diva ( )
看板MATLAB
标题Re: [绘图] 方程式绘图问题
时间Sat Mar 10 22:44:07 2012
你要的是不是类似这种结果
假设你的方程式是
y=x^3+2*x^2-10*x+10
图形为
http://imgur.com/BysVz
要用迭代的方法找y=0的解
opts=optimset('display','iter','tolx',0.2) ;
^^^^^^^^^^^^^^^ ^^^^^^^^^
显示迭代结果 tolerance <=0.2
fzero('x^3+2*x^2-10*x+10',10,opts) ;
^^
从10开始找
结果为
Search for an interval around 10 containing a sign change:
Func-count a f(a) b f(b) Procedure
1 10 1110 10 1110 initial
3 9.71716 1019.2 10.2828 1205.92 search
5 9.6 983.056 10.4 1247.18 search
7 9.43431 933.383 10.5657 1307.1 search
9 9.2 865.968 10.8 1394.99 search
11 8.86863 776.159 11.1314 1525.76 search
13 8.4 659.824 11.6 1724.02 search
15 7.73726 515.55 12.2627 2032.13 search
17 6.8 348.912 13.2 2526.45 search
19 5.47452 179.269 14.5255 3351.45 search
21 3.6 46.576 16.4 4794.86 search
23 0.949033 3.16576 19.051 7459.71 search
25 -2.8 31.728 22.8 12674 search
26 -8.10193 -309.52 22.8 12674 search
Search for a zero in the interval [-8.10193, 22.8]:
Func-count x f(x) Procedure
26 -8.10193 -309.52 initial
27 -4.86116 -9.00013 interpolation
28 -4.86116 -9.00013 interpolation
Zero found in the interval [-8.10193, 22.8]
这是你要的吗？
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.103.220
※ 编辑: diva 来自: 140.112.103.220 (03/10 22:44)
1F：推 luckykey:大大你好强! 我欠你一杯饮料 85度C乙张兑换卷送你! 03/10 22:48
2F：→ luckykey:我会研究一下你的写法! 方程式只是举例! 03/10 22:49
3F：推 luckykey:不过你的方程式是多项式 我要找x=x.^-2+1/x-1之类的 03/10 22:56
4F：→ diva:x=x.^-2+1/x-1 这也是多项式啊 不然改成 y=x^-2+1/x-1-x 03/11 09:26
5F：→ diva:不过不一定会跑出结果 不一定所有的多项式都有解 03/11 09:27
6F：→ diva:我指的是实数解 03/11 09:27