普通方程和微分方程-在線瀏覽
2024-08-03 23:58本頁面
【正文】 100 x0=y。 n=n+1。方程組的解y=39。yfprintf(39。)。迭代次數(shù)n=39。n在mand window輸入：A=[10 3 1。1 3 10]b=[14 5 14]39。y=gaussseidel(A,b,x0)收斂速度比Jacobi更快SOR迭代法：%該函數(shù)用SOR迭代法求解線性方程組%用戶需要輸入3個參數(shù)%再設(shè)A*y=b，用戶輸入矩陣A和b%再輸入一個初始向量x0和松弛系數(shù)w%使用格式為y=sor(A,b,w,x0)function y=sor(A,b,w,x0)D=diag(diag(A))。L=tril(A,1)。f=(Dw*L)\b*w。n=1。n=1000 x0=y。 n=n+1。方程組的解y=39。yfprintf(39。)。迭代次數(shù)n=39。n在mand window輸入：A=[10 3 1。1 3 10]b=[14 5 14]39。y=sor(A,b,x0)y=sor(A,b,x0)y=sor(A,b,x0)w對收斂速度影響很大非線性方程組的解法fsolve函數(shù)x=fsolve(fun,x0)：以x0為初始矩陣來求解方程fun，fun接受輸入量x并返回一個向量（矩陣），使得f=fun(x)x=fsolve(fun,x0,options)：以options為選擇參數(shù)的輸入變量，詳見help fsolvex=fsolve(fun,x0,options,
物理相關(guān)推薦
鄂ICP備17016276號-1