怎样用MATLAB求方程的解ln
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/14 07:22:27
symsxy=x*sin(x);dy=diff(y);%一阶导数d2y=diff(dy);%二阶导数x1=double(solve(dy));%找出导数为0的点,就是极点x2=subs(d2y,x1)
将方程的三对解分别代入上述方程,化简后可得到三个新的方程,联立该方程组,就可得到另外三个未知系数.但是由于分母中有个指数函数,且指数中有未知数,所以不一定能直接求出,需要变换形式.再问:求问怎么变换呢
eqn=@(x)4*x^4-4*x^2;fzero(eqn,0)fzero(eqn,-1.1)fzero(eqn,-1)fzero(eqn,-0.9)fzero(eqn,0.0)fzero(eqn,1
矩阵有两种乘法:点乘和插乘.比如矩阵A乘以矩阵B,matlab语言表示如下:点乘:A.*B插乘:A*B点乘为两个矩阵的对应项相乘,插乘为一般的矩阵乘法.
初始条件呢?告诉初始条件就可求解了.好像是导弹追击问题啊!比较复杂,调用maple了:maple('dsolve','')结果:ans=[,],[,]结果没有问题,请看例子:>>maple
由|sin(t)|>1/16可得:n*PI+arcsin(1/16)
%%定义xx为最终解x=solve('x*(10+30*sqrt(1+(20*x)^2))-y');i=1;foryy=0:0.01:1xx(:,i)=subs(x,'y',yy);i=i+1;end
用solve命令就可以了,具体用法你查看一下MATLAB的帮助文档吧比如:solve('2*x+y-z=1','2*y+3*z=-2','x-z=-1');上面的方程是我胡乱写的,意思是有几个方程就往
limit(F,x,a)计算当x→a时符号表达式F=F(x)的极限值limit(F,a)用函数findsym(F)确定F中的自变量x,再计算当x→a时F=F(x)的极限值limit(F)用函数find
ank(A)即可再问:能具体点吗?从来没用过.再答:输入矩阵A=[12;34]回车输入rank(A)回车出结果>>a=[12;34]a=1234>>rank(a)ans=2
输入下面4行命令即可.>>h=29.8149;b=68.4982*0.5;>>symsx>>123.5;hr=h/b;xr=x/b;A=1/((hr^2+xr^2)^0.5);B=hr/sqrt(1+
symsaxyy=int(x^2+a,x,0,a)然后运行solve('1/3*a^3+a^2=3','a');vpa(a)这结果出来了,三个.自己根据需要选择合适解.一般是第一个,后面俩个是复数解
先用数值拟合求出其时域方程,在用ss2tf求出其传输函数.
%bydynamic%2009.2.%积分Matlab中提供了数值和解析解法两种1.数值方法>>fun=@(x)738.5*x.^2.*exp(-0.00105*x.^2);>>quadl(fun,0
编程观察了下F1(t)/F0(t)-t函数,图像如下:也就是说a/K需要是个正数才行.你给的数据a是正的,K是负的,无解.给定a,K,计算yita的程序是:%%示例a=1;K=2;Fn=@(n,yit
f1=sym('x1+x01=sqrt(x2^2+y2^2)*cos(a+atan(y2/x2))+qrt(x01^2+y01^2)*cos(a+atan(y01/x01))');f2=sym('y1
给出方程看看吧.如果方程没有什么特殊性的话,那就基本等于靠猜了;而且一般都是无穷多解,不知道你要的“多组解”有什么要求没有.
symsxyy=x^3+4*x^2+5*x+6;solve('y=x^3+4*x^2+5*x+6','x')ans=1/6*(-440+108*y+12*(1344-660*y+81*y^2)^(1/
参考代码:>> syms w kc real>> Tjw = 7.8*kc/(7.8*kc-0.11*w^
对于非线性方程,可用nlinfit()函数来求解.求解步骤如下:x=[.]y=[.]z=[.]X=[xy]func=inline('a(1).*exp(a(2)*x(1))*(1+exp(-a(2)*