91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 1KB
    文件類型: .rar
    金幣: 2
    下載: 1 次
    發布日期: 2021-08-08
  • 語言: Matlab
  • 標簽: 對分法??MATLAB??

資源簡介

數值分析中,用對分法求方程根的完整MATLAB代碼。首先是對分法的完整代碼,,文件中舉了一個函數的例子函數可以更改。是數值分析課程很好的實驗案例

資源截圖

代碼片段和文件信息

function?bisec_g(f_nameacxminxmaxn_points)
% ac:end?points?of??initial??interval
% tolerance :?tolerance
% it_limit:limit?of?iteration?number
% Y_a?Y_c?:?y?values?of?the?current?end?points?
clfhold?off
clear?Y_aclear?Y_c
wid_x=xmax-xmin;
dx=(xmax-xmin)/n_points;
xp=xmin:dx:xmax;
yp=feval(f_namexp);
plot(xpyp‘r‘);????????????????????????%?y=?f(x)
xlabel(‘x‘);ylabel(‘f(x)‘);?title(‘Bisection?Method‘)hold?on
ymin=min(yp);ymax=max(yp);wid_y=ymax-ymin;
yp=0.0*xp;
plot(xpyp)?????????????????????????????%?x?axes?
fprintf(‘Bisection?Scheme?\n\n‘);
fprintf(‘It ?a ?b ?c fa=?f(a)‘);
fprintf(‘ fc=?f(c)?abs(fc-fa) abs(c-a)/2\n‘);
tolerance=0.000001;
it_limit?=?30;
it=0;
Y_a=feval(f_namea);
Y_c=feval(f_namec);
plot([aa][Y_a0]‘black‘);text(a-0.1*wid_y‘x=a‘)???%line
plot([cc][Y_c0]‘black‘);text(c-0.3*wid_y‘x=c‘)???%line
if(Y_a*Y_c>0)
????fprintf(‘f(a)f(c)>0\n‘);
else
????while?1
????????it=it+1;
????????b=(a+c)/2;
????????Y_b=feval(f_nameb);
????????plot([bb][Y_b0]‘:‘);???%line
????????plot(b0‘o‘)???????????????%o?point
?if?it<4
?????text(bwid_y/20[num2str(it)])
?end
????fprintf(‘%3.0f?%10.6f?%10.6‘itab);
????fprintf(‘%10.6f?%10.6f?%10.6‘cY_aY_c);
????fprintf(‘%12.3e?%12.3e\n‘abs(Y_c-Y_a)abs(c-a)/2);
????if(abs(c-a)/2<=tolerance)
????????fprintf(‘Tolerance?is?satisfied.\n‘);
????????break
????end
???if(it>it_limit)
???????fprintf(‘Iteration?limit?exceeded.\n‘);
???????break
???end
???if??(Y_a*Y_b<=0)
???????c=b;Y_c=Y_b;
???else
???????a=b;Y_a=Y_b;
???end
????end
fprintf(‘Final?result;Root=%12.6f\n‘b)
end
x=b;
plot([xx][0.05*wid_y0.2*wid_y]‘g‘);text(x0.25*wid_y‘Final?solutin‘)
plot([x(x-wid_x*0.004)][0.05*wid_y0.09*wid_y]‘r‘)????%?arrow?line
plot([x(x+wid_x*0.004)][0.05*wid_y0.09*wid_y]‘r‘)???%?arrow?line
????????????
????????
????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1872??2017-12-06?10:36??對分法求根\bisec_g.m

?????文件?????????41??2017-12-06?10:01??對分法求根\ex_5x2.m

?????文件?????????29??2017-12-06?09:59??對分法求根\main.m

?????目錄??????????0??2017-12-06?10:37??對分法求根

-----------?---------??----------?-----??----

?????????????????1942????????????????????4


評論

共有 條評論