[繪圖] 內插法練習

看板MATLAB作者 (文化裘德洛)時間7年前 (2017/03/02 23:33), 7年前編輯推噓1(100)
留言1則, 1人參與, 最新討論串1/1
先貼我的程式碼 r = [9.37 12.32 12.46 9.51 5.92 4.33 4.29 3.88 4.51 10.92 27.50 49.67 69.59 81.73 88.19 86.05]; w = linspace(400,700,16); plot(w,r,'bo'); for i=3:13 a0 = r(i); a1 = (2*r(i-2)-16*r(i-1)+16*r(i+1)-2*r(i+2))/24; a2 = (-r(i-2)+16*r(i-1)-30*r(i)+16*r(i+1)-r(i+2))/24; a3 = (-9*r(i-2)+39*r(i-1)-70*r(i)+66*r(i+1)-33*r(i+2)+7*r(i+3))/24; a4 = (13*r(i-2)+64*r(i-1)+126*r(i)-124*r(i+1)+61*r(i+2)-12*r(i+3))/24; a5 = (-5*r(i-2)+25*r(i-1)-50*r(i)+50*r(i+1)-25*r(i+2)+5*r(i+3))/24; wave = 400+(i-1)*20; w = linspace(wave+2,wave+18,19); x = linspace(0.1,0.9,19); y = a0+a1*x+a2*x.^2+a3*x.^3+a4*x.^4+a5*x.^5; hold on plot(w,y,'ro') end 小弟想要把圖畫成這樣 http://imgur.com/a/4wIxA 但我打出來後會變這樣 http://imgur.com/a/VkTIA 我的程式碼是按照書上打的,也檢查很多次都沒錯 請問各位神人可以幫我看看哪邊出錯嗎 拜託了 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 49.217.178.122 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1488468811.A.A17.html

03/03 01:28, , 1F
a4的係數是-64還是+64?檢查好多次都沒錯呢!
03/03 01:28, 1F
!!我真的眼殘 謝謝k大 ※ 編輯: johnny1229 (49.217.178.122), 03/03/2017 02:04:52
文章代碼(AID): #1Ok3jBeN (MATLAB)
文章代碼(AID): #1Ok3jBeN (MATLAB)