實(shí)驗(yàn)非線性方程求解

1、94 實(shí)驗(yàn)?zāi)繉?shí)驗(yàn)內(nèi)分別用 fzero 和 94 實(shí)驗(yàn)?zāi)繉?shí)驗(yàn)內(nèi)分別用 fzero 和 fsolve 程序求方程 sinx-x2/2=0 的所有根，準(zhǔn)確到 10-去不同的初始值計(jì)算，輸出初值、根的近似解和迭代次數(shù)，分析不同根的收斂域；自己構(gòu)造某個(gè)迭代公式（如 x=(2sinx)1/2 等）用迭代法求解，并自己編寫法的程序進(jìn)行求解和比較分析設(shè)計(jì)2知有兩個(gè)根，分布在 x=0 x=1.4 1結(jié)2法xx0000000004031112121155656-8.86E-7.83E-0556583356668-6.02E-2.60E-06656755553.51E-1.16E-7867775.76E-9.89E

2、-767-048672.34E-1.73E-778-4.19E-1.61E-0568-5.08E-79-5.08E-結(jié)2法xx0000000004031112121155656-8.86E-7.83E-0556583356668-6.02E-2.60E-06656755553.51E-1.16E-7867775.76E-9.89E-767-048672.34E-1.73E-778-4.19E-1.61E-0568-5.08E-79-5.08E-分在題中以-100000 和 1000000 作為 分在題中以-100000 和 1000000 作為 fzero迭代法：自己構(gòu)造的迭代法存在明顯的弊端

3、，0 0 0 32結(jié) 參考程序opt= %input the initial value 3 x(2)=input(initial value x(2): ); % clearx(1)=input(pleaseinputtheinitialvaleu:); clearallx(1)=input(pleaseinputtheinitialvaleu:); 4給定4 種物質(zhì)對應(yīng)的參數(shù)ai,bi,ci 和交互作用矩陣Q 如下：a1 b1 c1 a2 給定4 種物質(zhì)對應(yīng)的參數(shù)ai,bi,ci 和交互作用矩陣Q 如下：a1 b1 c1 a2 b2 c2 a3 b3 c3 a4 b4 c4 Q 在壓強(qiáng)下，

4、為了形成均相共沸混合物，溫度和組分分別是少？請盡量找出所有可能的解。分析設(shè)計(jì)分析：組分控制為 0，觀察其他組分能否實(shí)現(xiàn)均相共沸混合物。結(jié)果：XT =0.4142 71.9657 Y = 1.0e-006 *-0.0009 -0.4428 -XT= -0.0000 Y= 1.0e-005 0.0056 -結(jié)論：71.965758.58%5參考程序functionfor i=1:n-1 for參考程序functionfor i=1:n-1 ford(i)= x*forf(i)=x(i)*(b(i)/(T+c(i)+log(x*Q(i,1:n)+dd*Q(1:n,i)-a(i)-1+p); 0.31

5、61.00.4770.3770.361.0 for i=1:100 7. 用迭代公式xk+1=a xk exp(-bx)計(jì)算序列x，分析其收斂性，其a分取 5,11,15b（0）任意，初值 x=1.觀察是否有混沌現(xiàn)象出現(xiàn)并找出幾個(gè)分岔點(diǎn)，觀察分岔點(diǎn)的極限趨勢是否符合Feigenbaum常數(shù)揭示的規(guī)律。分析設(shè)計(jì)結(jié)果：6矩陣X形成表格如下：A 取不同值的曲線圖如下k0111123456789矩陣X形成表格如下：A 取不同值的曲線圖如下k0111123456789分析：分析：結(jié)論：xk+1=axkexp(-x=axexp(-一個(gè)平（ |f (y |f (y f(x)= x= a x exp(-bx) y*。根據(jù)迭代法拒（ 1-ln(a) =1-（ （ a=5 x k8 a=11、15xkxk 。混沌現(xiàn)象（參考程序如下 a=11、15xkxk ?；煦绗F(xiàn)象（參考程序如下a，縱坐標(biāo)表示迭代序列的極限）（由該圖也可以看出收斂