1樓:半條渴死的小魚
function newton(x0,e,n)%輸入xo為估計的迭代初值,e為規定的誤差,n為最大迭代次數.
%輸出x,y為最後迭代的兩個近似根,k為迭代次數.
clcformat long;
disp('迭代次數 近似根')k=0;
x1=0;
x2=x0;
while (abs(x2-x1))>e
x1=x2;
x2=x1-f(x1)./df(x1);
k=k+1;
if k>n
return;
end%%%%%%記錄並輸出%%%%%%%%%o1=sprintf('%3d',k);
o2=sprintf('%3.8f',x2);
ol=[o1,' ' o2];
disp(ol);
y(k)=x2;
end%%%%畫圖%%%%%%%
i=1:k;
figure(2)
plot(i,y,'rd-')
grid on
xlabel('迭代次數')
ylabel('近似根')
title(['牛頓法求出的該方程的近似根 x^*=', num2str(x2,9)])
function y=f(x)
y=x^2/2-sin(x)-1;
function y=df(x)
y=x-cos(x);
2樓:晰美酒窩
%求方程x^3-sinx-12x+1=0的根。
n =1000;
x0=0;
x1=x0^3-sin(x0)-12*x0+1;
e=1.0e-6;%f=x*x*x+4*x*x-10=0,f1為f=x*x*x+4*x*x-10的一階導
k=1;
用matlab做,牛頓迭代法
3樓:sky不用太多
function [ a ] = cal( a,b,v )%a,b表示區間,v是精度
i=1;
x = (a+b)/2;
a=[i x];
t = x-(x^3-x-1)/(3*x^2-1);%迭代函式
while(abs(t-x)>v)
i=i+1;
x = t;
a = [a;i x];
t = x-(x^3-x-1)/(3*x^2-1);%迭代函式
enda = [a;i+1 t];
end執行結果:
>> format long;
>> cal(1,2,0.00001)
ans =
1.000000000000000 1.500000000000000
2.000000000000000 1.347826086956522
3.000000000000000 1.325200398950907
4.000000000000000 1.324718173999054
5.000000000000000 1.324717957244790
牛頓迭代法(newton's method)又稱為牛頓-拉夫遜(拉弗森)方法(newton-raphson method),它是牛頓在17世紀提出的一種在實數域和複數域上近似求解方程的方法。多數方程不存在求根公式,因此求精確根非常困難,甚至不可能,從而尋找方程的近似根就顯得特別重要。方法使用函式f(x)的泰勒級數的前面幾項來尋找方程f(x) = 0的根。
牛頓迭代法是求方程根的重要方法之一,其最大優點是在方程f(x) = 0的單根附近具有平方收斂,而且該法還可以用來求方程的重根、復根,此時線性收斂,但是可通過一些方法變成超線性收斂。另外該方法廣泛用於計算機程式設計中。
已經證明,如果是連續的,並且待求的零點是孤立的,那麼在零點周圍存在一個區域,只要初始值位於這個鄰近區域內,那麼牛頓法必定收斂。 並且,如果不為0, 那麼牛頓法將具有平方收斂的效能. 粗略的說,這意味著每迭代一次,牛頓法結果的有效數字將增加一倍。
[1]迭代法也稱輾轉法,是一種不斷用變數的舊值遞推新值的過程,跟迭代法相對應的是直接法(或者稱為一次解法),即一次性解決問題。迭代演算法是用計算機解決問題的一種基本方法。它利用計算機運算速度快、適合做重複性操作的特點,讓計算機對一組指令(或一定步驟)重複執行,在每次執行這組指令(或這些步驟)時,都從變數的原值推出它的一個新值。
利用迭代演算法解決問題,需要做好以下三個方面的工作:
一、確定迭代變數
在可以用迭代演算法解決的問題中,至少存在一個可直接或間接地不斷由舊值遞推出新值的變數,這個變數就是迭代變數。
二、建立迭代關係式
所謂迭代關係式,指如何從變數的前一個值推出其下一個值的公式(或關係)。迭代關係式的建立是解決迭代問題的關鍵,通常可以使用遞推或倒推的方法來完成。
三、對迭代過程進行控制
在什麼時候結束迭代過程?這是編寫迭代程式必須考慮的問題。不能讓迭代過程無休止地執行下去。
迭代過程的控制通常可分為兩種情況:一種是所需的迭代次數是個確定的值,可以計算出來;另一種是所需的迭代次數無法確定。對於前一種情況,可以構建一個固定次數的迴圈來實現對迭代過程的控制;對於後一種情況,需要進一步分析得出可用來結束迭代過程的條件。
4樓:匿名使用者
syms x
f=x^x-10;
df=diff(f,x);
eps=1e-6;
x0=10;
cnt=0;
maxcnt=200; %最大迴圈次數
while cnt if (abs(x1-x0) break; endx0=x1; cnt=cnt+1; endif cnt==maxcnt disp '不收斂' else vpa(x1,8)end 用matlab如何編寫牛頓迭代法問題,謝謝 5樓:chenbd文庫 看一本書叫數值計算, x=x^3-1或者 x=sqrt3(x+1) 不停的迭代即可 其中一個應該滿足迭代的收斂原則,具體的收斂條件請看書。 6樓:匿名使用者 function [ a ] = cal( a,b,v )%a,b表示區間,v是精度 i=1; x = (a+b)/2; a=[i x]; t = x-(x^3-x-1)/(3*x^2-1);%迭代函式 while(abs(t-x)>v) i=i+1; x = t; a = [a;i x]; t = x-(x^3-x-1)/(3*x^2-1);%迭代函式 enda = [a;i+1 t]; end執行結果: >> format long; >> cal(1,2,0.00001) ans = 1.000000000000000 1.500000000000000 2.000000000000000 1.347826086956522 3.000000000000000 1.325200398950907 4.000000000000000 1.324718173999054 5.000000000000000 1.324717957244790 include int main printf d r n s return 0 c語言程式設計1 2 3 4 5 6 99 100 1 2 3 4 5 6 99 100 1 2 3 4 5 6 99 100 1 100 2 這裡指一共有50組的意思,將100個數分成50組,每組兩個數相減 50 為... 1 電bai力系統 matab的simulink工具箱裡面du,有simpowersystem模組,專zhi 門為電力系dao統提供專各種 特點就是簡單屬,很多模型都有。2 程式設計計算。設計演算法的,都可以用matlab程式設計計算,方便快速。多用於simulink 用途不多 matlab對電氣工... k menu please choose a number 1 4 switch k case 1 a 10 b 10 result a b a rand 1 case 4 a 0 b 1 result a sqrt b randn 1 enddisp result 看看du 這個題,zhi吧,或許...求12345699100用c語言編出來
matlab在電氣專業用多嗎,Matlab對電氣工程及其自動化專業的人具體用途
在matlab中switch怎麼用