MATLAB程序(牛顿法及线形方程组)【DOC精选】.doc

MATLAB程序 1、图示牛顿迭代法（M文件）文件名：newt_g function x = new_g(f_name,x0,xmin,xmax,n_points) clf,hold off % newton_method with graphic illustration del_x = 0.001; wid_x = xmax - xmin; dx = (xmax - xmin)/n_points; xp = xmin:dx:xmax; yp = feval(f_name,xp); plot(xp,yp);xlabel(x);ylabel(f(x)); title(newton iteration),hold on ymin = min(yp); ymax = max(yp); wid_y = ymax-ymin; yp = 0. * xp; plot(xp,yp) x = x0; xb = x+999; n=0; while abs(x-xb) 0.000001 if n 300 break; end y=feval(f_name,x); plot([x,x],[y,0]);plot(x,0,o) fprintf( n = % 3.0f, x = % 12.5e, y = % 12.5e \ n, n, x, y); xsc = (x-xmin)/wid_x; if n 4, text(x,wid_y/20,[num2str(n)]), end y_driv = (feval(f_name,x + del_x) - y)/del_x; xb = x; x = xb - y/y_driv; n = n+1; plot([xb,x],[y,0]) end plot([x x],[0.05 * wid_y 0.2 * wid_y]) text( x, 0.2 * wid_y, final solution) plot([ x ( x - wid_x * 0.004)], [0.01 * wid_y 0.09 * wid_y]) plot([ x ( x + wid_x * 0.004)], [0.01 * wid_y 0.09 * wid_y]) 传热问题 假设一个火炉是用厚度为0.05m的砖单层砌成的。炉内壁温度为T0=625K, 外壁温度为T1（未知）。由于对流和辐射造成了外壁的热量损失，温度T1由下式决定： 其中： k：炉壁的热传导系数，1.2W/mK : 发射率，0.8 T0：内壁温度，625K T1：外壁温度（未知），K ：环境温度，298K Tf：空气温度，298K H：热交换系数，20W/m2K :Stefan-Boltzmann常数，5.67*10-8W/m2K4 :炉壁厚度，0.05m 试用牛顿迭代法求解T1. 先根据待求解的方程式建立（M文件）文件名：wall_ht function f= wall_ht(t1) k = 1.2; e = 0.8; tinf = 298; tf = 298; h = 20; t0 = 625; sig = 5.67e-8; wall_thick = 0.05; f = k/wall_thick * (t1-t0) + e * sig * (t1.^4-tinf^4) + h * (t1 - tf); 程序的调用及计算结果： newt_g(wall_ht,550,400,600,500) n = 0, x = 5.50000e+002, y = 7.03301e+003 n = 1, x = 4.55199e+002, y = 6.58551e+002 n = 2, x = 4.44423e+002, y = 6.44623e+000 n = 3, x = 4.44316e+002, y = 6.27680e-004 n = 4, x = 4.44316e+002, y = 5.67525e-010 ans = 444.3157 平衡问题 一氧化碳与氢按以下反应生成甲醇 现有1molCO与2molH2的混合物，在温度t=5900C，压力p=3.04*107Pa下进行反应并达到平衡（Kf=1.393*10-15，Kr=0.43）,求CH3OH在平衡气中的摩尔分数。 达到平衡时的关系式： 其中 据待求函数建立（M文件）文件名：ping_hen function f= ping_heng(x) t = 590; p = 3.04e7; kf = 1.393e-15; k