24小时热门版块排行榜

返回列表

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

[求助] 求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序

求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序

回复此楼

» 猜你喜欢

招博士已经有4人回复
救命帖已经有7人回复
限项规定已经有6人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复
英文综述是否需要润色及查重已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
最失望的一年已经有18人回复

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

如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！已经有11人回复
求用牛顿迭代法解非线性方程组的matlab程序已经有4人回复
请教matlab如何增加迭代次数已经有11人回复
求用牛顿迭代法求非线性方程组已经有4人回复
matlab数值计算时迭代项分母为零的时候怎么办？已经有15人回复
matlab菜鸟问3个问题~~~~~~~~~~~ 已经有2人回复
求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序已经有5人回复
求助matlab画出螺旋线的程序已经有6人回复
【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】已经有15人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼 2011-10-15 21:36:58

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

锐利的碎片

木虫 (正式写手)

star watcher

应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
jjdg(金币+2): 感谢支持 2011-10-16 23:35:49
jjdg: 图真漂亮！牛！ 2011-10-16 23:36:10

CODE:

import numpy as np

import pylab as pl

from matplotlib import cm

def iter_point(c):

    z=c

    tol=1.0e-2

    for i in range(1,1000):

        err=-(z**3-1)/(3*z**2)

        if abs(err)
            break

        z=z-(z**3-1)/(3*z**2)

    return i

def draw_newton(cx,cy,d):

    x0,x1,y0,y1=cx-d,cx+d,cy-d,cy+d

    y,x=np.ogrid[y0:y1:200j,x0:x1:200j]

    c=x+y*1j

    newton=np.frompyfunc(iter_point,1,1)(c).astype(np.float)

    pl.imshow(newton,cmap=cm.Blues_r,extent=[x0,x1,y0,y1])

    pl.gca().set_axis_off()

    return

if __name__=='__main__':

    x,y=0.0,0.0

    draw_newton(x,y,0.5)

    pl.show()

赞一下(2人)

回复此楼

2楼2011-10-16 10:04:54

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

相关版块跳转我要订阅楼主我要飞的主题更新

返回列表