MATLAB第五课 代数方程求解.ppt

《MATLAB第五课 代数方程求解.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第五章：代数方程求解和其他符号函数1、解基本代数方程2、二次方程求解3、符号方程绘图4、高阶方程求解5、解方程组6、方程展开与合并7、使用指数和对数函数求解方程8、函数的级数表示9、习题1解基本代数方程>>x+3=0???x+3=0

2、Error:Theexpressiontotheleftoftheequalssignisnotavalidtargetforanassignment.>>x=3x=3用solve命令解方程>>x=solve('x+3=0')x=-3>>x=solve('x+3')x=-3当你传递给solve函数x+3时，MA

3、TLAB假设你传递的就是x+3=0提问：如果你求解的方程有时候可能含有多个符号，如ax+5=0，这个时候怎么求解？提问：如果你在MATLAB中输入方程ax+5=0，MATLAB只认为求解x。>>solve('a*x+5')ans=-5/a>>solve('a*x+5','a')ans=-5/x提问：上面研究了简单的一次方程，那么遇到二阶、三阶等方程怎么办呢？2二次方程求解例如求解：x2-6x-12=0？>>s=solve('x^2-6*x-12=0')s=3+21^(1/2)3-21^(1/2)>>y=3+s(1)y=6+21^(1/2)>

4、>s(1)-s(2)ans=2*21^(1/2)>>d=‘x^2+9*x-7=0’；>>solve(d)ans=-9/2+1/2*109^(1/2)-9/2-1/2*109^(1/2)把方程赋给一个变量，然后把它传递给solve也完全是可以的提问：上面解决了方程组求解问题，那么怎么作图呢？3符号方程绘图用ezplot（d），自动显示标题>>d='x^2+9*x-7';>>ezplot(d)>>d='x^2+9*x-7';>>ezplot(d,[-2,8])使用ezplot注意的地方>>ezplot('x+3=0')???Errorusing

5、==>inlineevalErrorininlineexpression==>x+3=0???Error:Theexpressiontotheleftoftheequalssignisnotavalidtargetforanassignment.Errorin==>inline.fevalat34INLINE_OUT_=inlineeval(INLINE_INPUTS_,INLINE_OBJ_.inputExpr,INLINE_OBJ_.expr);Errorin==>specgraphprivateezplotfevalat54z=f

6、eval(f,x(1));Errorin==>ezplot>ezplot1at448[y,f,loopflag]=ezplotfeval(f,x);Errorin==>ezplotat148[hp,cax]=ezplot1(cax,f{1},vars,labels,args{:});>>ezplot('x+3')如果：-4>eq='x^2+x-2^(1/2)';>>s=solve(eq)s=-1/2+1/2*(1+4*2^(1/2))^(1/2)-1/2-1/2*(1+4*2^(1/2))^(1

7、/2)>>x=double(s(1))x=0.7900>>ezplot(eq)4高阶方程求解求（x+1）2（x-2）=0>>eq='(x+1)^2*(x-2)';>>solve(eq)ans=2-1-1练习：求x4-5x3+4x2-5x+6=0，并绘制-10>eq1='x^4-5*x^3+4*x^2-5*x+6';>>s=solve(eq1);>>a=s(1)a=5/4+1/12*3^(1/2)*((43*(8900+12*549093^(1/2))^(1/3)+2*(8900+12*549093^(1/2))^(2

8、/3)+104)/(8900+12*549093^(1/2))^(1/3))^(1/2)+1/12*((258*(8900+12*549093^(1/2))^(1/3)*((43*(8900+12*549093^(1/2))^(1/3)+2*(8900+12*549093^(1/2))^(2/3)+104)/(8900+12*549093^(1/2))^(1/3))^(1/2)-6*((43*(8900+12*549093^(1/2))^(1/3)+2*(8900+12*549093^(1/2))^(2/3)+104)/(8900+12*5

9、49093^(1/2))^(1/3))^(1/2)*(8900+12*549093^(1/2))^(2/3)-312*((43*(8900+12*549093^(1/2))^(1