24小时热门版块排行榜

>论坛更新日志 (3561)
>虫友互识 (437)
>休闲灌水 (219)
>导师招生 (199)
>文献求助 (157)
>考博 (60)
>论文道贺祈福 (56)
>博后之家 (53)
>论文投稿 (53)
>硕博家园 (48)
>考研 (44)
>基金申请 (34)
>公派出国 (29)
>招聘信息布告栏 (28)
>找工作 (19)
>教师之家 (15)

返回列表

我要飞

铁虫 (正式写手)

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

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

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

回复此楼

» 猜你喜欢

过年走亲戚时感受到了所开私家车的鄙视链已经有10人回复
今年春晚有几个节目很不错，点赞！已经有12人回复
情人节自我反思：在爱情中有过遗憾吗？已经有13人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有12人回复
基金正文30页指的是报告正文还是整个申请书已经有5人回复

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

如何用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的回帖

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

返回列表