利用MATLAB求方程x-e^(-x)=0的解

x=0:0.1:10;%x的取值范围,可以为任意的.y=exp(x);%e指数函数p=polyfit(x,y,5)%五次多项式去拟合e指数函数,得到的p是五次多项式的系数.

用corrcoef函数设a1,b1,c1,d1,a2,b2,c2,d2分别为f(x)和g（x）的系数x=[a1,b1,c1,d1];y=[a2,b2,c2,d2];z=corrcoef(x,y)

先建立二分法的fun.m文件,代码如下：functionfun(a,b,e)%f是自定义的函数%a为隔根区间左端点,b为隔根区间右端点,e为绝对误差限ifnargin==2e=1.0e-6;elsei

以下仅作参考,若大家有更好的方法,欢迎交流1.若手动挑选方程的解,可以这样写fx=inline('5*x.^2.*sin(x)-exp(-x)');x0=fsolve(fx,0:10)y=subs(f

symsx>>fun=inline(cos(x)-x.*exp(x)','x')fun=Inlinefunction:fun(x)=cos(x)-x.*exp(conj(x))>>fzero(fun,

symsaf=((2.72.^a-a-1).^(-1/2)-1);g=int(f,a,a,3)结果：g=int(1/((68/25)^a-a-1)^(1/2)-1,a=a..3)

程序如下：clear,clc;a=0;%a=input('inputa:');b=1;%b=input('inputb:');err=10^-5;y1=a*exp(a)

源代码:在matlab中保存为:demimethod.mfunction[x,k]=demimethod(a,b,f,emg)%a,b:求解区间的两个端点%f:所求方程的函数名%emg：精度指标%x：

试试我这个程序,在附件中

程序如下：clear,clc;a=0;%a=input('inputa:');b=1;%b=input('inputb:');err=10^-5;y1=a*exp(a)-1;y2=b*exp(b)-1

令f(x)=e^x+x-2;f(0)=-1f(1)=e+1-2>0f'(x)=e^x+1>0=>所以f(x)单调存在0

f=@(x)5*x+exp(-x)；x=fsolve(f,0)

symsaxyx=[1223555766];y=[25.634789497105];a=polyfit(x,y,2);%这是求拟合代数式的系数a=polyfit(x,y,1)系数为：1.129233.

这个定义一个函数.f=@(x)x^3-3*x^2-x+3;的意思是说定义一个函数f(x),它只有一个自变量.使用时,直接可以用f(1)它就是x=1时的函数值.另外多参数可以是：f=@(x,y)sqrt

请到我的百度空间,标题为TOchenyan52288

function y=wtxy(beta,x)       %保存为 wtxy.m 文件y=beta

用solve函数可以解,symsabcs=solve(a*exp(b)+c-9.846659,a*exp(4*b)+c-7.03417,a*exp(8*b)+c-4.489171)s=a:[6x1sy

f=@(x)x.*sin(x)-1/2;x=fsolve(f,[-10:1:10])x=Columns1through6-9.3714-9.3714-6.3619-6.3619-2.9726-2.97

不定积分举例>>symsx>>int(x^2)ans=1/3*x^3>>int(cos(x))ans=sin(x)

你这个有点像有限元里面的东西.下面的K是我乱写的.你如果给定的不是前三个位移后三个应力的话,你自己得适当地调整,实际上就是调整下面的K1和K2,你一个命令地看,应当可以OK的.K=[2,1,0,0,0