版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

我要飛

鐵蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1515.3
散金: 144
紅花: 5
帖子: 458
在線: 591.2小時
蟲號: 671362
注冊: 2008-12-09
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

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

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

回復(fù)此樓

» 猜你喜歡

26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工（0856）304求 B區(qū) 調(diào)劑已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有8人回復(fù)
265求調(diào)劑已經(jīng)有14人回復(fù)
350 本科985求調(diào)劑，求老登收留已經(jīng)有3人回復(fù)
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有9人回復(fù)
初試 317 已經(jīng)有3人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有15人回復(fù)
279求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

如何用matlab將這個方程代表的圖畫出來，最好能把畫圖的代碼貼出來，謝謝��！已經(jīng)有11人回復(fù)
求用牛頓迭代法解非線性方程組的matlab程序已經(jīng)有4人回復(fù)
請教matlab如何增加迭代次數(shù) 已經(jīng)有11人回復(fù)
求用牛頓迭代法求非線性方程組已經(jīng)有4人回復(fù)
matlab數(shù)值計(jì)算時迭代項(xiàng)分母為零的時候怎么辦？已經(jīng)有15人回復(fù)
matlab菜鳥問3個問題~~~~~~~~~~~ 已經(jīng)有2人回復(fù)
求助matlab里面畫出x^3-1=0的Newton迭代的吸引域的程序已經(jīng)有5人回復(fù)
求助matlab畫出螺旋線的程序已經(jīng)有6人回復(fù)
【求助】求助：vb編程中用牛頓迭代解三次方程為什么只得到一個根？【已完成】已經(jīng)有15人回復(fù)
【求助】求一個用牛頓迭代法解非線性方程組的程序已經(jīng)有10人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

銳利的碎片

木蟲 (正式寫手)

star watcher

應(yīng)助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時
蟲號: 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
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人)

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主我要飛的主題更新