matlab松壹元叁次方程,违反掉落的邑是副数根。

时间:2018-09-22      作者:admin

  关于你昨天的阿谁方程,是典型的匪线性方程,我修改了以次,用fsolve函数到来松,以次不才面。先新建壹个m文件,又把以下整顿个以次直接骈制进m文件中,运转即却,不要在主窗口中运转,用m文件到来写以次,比直接在主窗口中操干便宜,也算是壹种好习惯吧。

  由计算结实却知,松更适宜还愿,同时没拥有拥有出产即兴虚数。

  function solution_nonl

  clear all;clc

  global aa bb tt pp

  %-------------根据计算结实,修改迭代初值------------------------------------

  xx0=[1]; %迭代初值,假设计算结实不快宜还愿,却修改,普畅通由pV=nRT给出产初始的V,此雕刻边我直接取1了

  t=[350.1500 348.9500 347.1500 345.4500 342.3500 340.7500 338.9500];

  p=[101.3300 101.3300 101.3300 101.3300 101.3300 101.3300 101.3300 101.3300];

  a=[0.9098 0.9365 0.9763 1.0223 1.0549 1.0968 1.1340 1.1773];

  b=[0.0573 0.0565 0.0551 0.0535 0.0523 0.0507 0.0491 0.0470];

  v=[];

  for n=1:7

  aa=a(n);

  bb=b(n);

  tt=t(n);

  pp=p(n);

  x(n)=fsolve(@nonlin,xx0);

  end

  disp(x)

  function f=nonlin(x)

  global aa bb tt pp

  f=8.314*tt/(x-bb)-aa/(x*(x+bb))-pp;

  计算结实:

  28.7864 28.6871 28.5380 28.3969 28.1413 28.0084 27.8591,

标签: