MATLAB(陈怀琛)20111018.ppt

* * 线性方程组的求解 线性方程组的形式为A*X=B；其中A至少行满秩。 X=linsolve(A,B) 输出方程的特解X。 例如：键入 A=sym([cos(t),sin(t);sin(t),cos(t)]); B=sym([1;1]); c=linsolve(A,B) c = [ 1/(sin(t)+cos(t))] [ 1/(sin(t)+cos(t))] * * 代数方程的求解 solve(P,v) 对方程P中的指定变量v求解。v可省略。 solve(P1,P2,…，Pn,v1,v2,…，vn) 对方程P1，P2，…Pn中的指定变量v1, v2…vn求解。 * * 例：可输入 solve(p+sin(x)=r) 得： ans =-asin(p-r) 例输入：P1=x^2+x*y+y=3;P2=x^2-4*x+3=0; [x,y]=solve(P1,P2) 得： x = [ 1; 3] y = [1; -3/2] 可输入：P1=a+u^2+v^2=0;P2=u-v=1; [u,v]=solve(P1,P2,u,v) 得： u = [ 1/2+1/2*(-1-2*a)^(1/2); 1/2-1/2*(-1-2*a)^(1/2)] v =[ -1/2+1/2*(-1-2*a)^(1/2) ;-1/2-1/2*(-1-2*a)^(1/2)] * * 对于有些无法求出解析解的非线性方程组，MATLAB只给出一个数值解。这一点可以 从表示解的数字不被方括号括住而确定。例如：键入： [x,y]=solve(sin(x+y)-exp(x)*y=0,x^2-y=2) 得： x = -6.0173272500593065641097297117905 y = 34.208227234306296508646214438330 由于这两个数字没有被[ ]括住，所以它们是数值解。 * * 另外，可利用solve来解线性方程组的通解。例如：键入 P1=2*x1+7*x2+3*x3+x4=6; P2=3*x1+5*x2+2*x3+2*x4=4; P3=9*x1+4*x2+x3+7*x4=2; u=solve(P1,P2,P3,x1,x2,x3,x4) ? Warning: 3 equations in 4 variables. u = x1: [1x1 sym] x2: [1x1 sym] x3: [1x1 sym] x4: [1x1 sym] * * 可以看到：屏幕提示“有3个方程4个变量”，意为解不唯一。（有时会提示解不唯一）且输出的是解的结构形式。为进一步得到解，可输入： u.x1,u.x2,u.x3,u.x4, 得： ans = x1 ans = -5*x1-4*x4 ans = 11*x1+9*x4+2 ans = x4 这样就得到了原方程组的通解。 * * 解符号微分方程 dsolve(eq1,eq2,…) 其中eq表示相互独立的常微分方程、初始条件或 指定的自变量。默认的自变量为t。 如果输入的初始条件少于方程的个数，则在输出结果中出现常数c1,c2,等字符。 关于微分方程的表达式有如下的约 定： 字母y表式函数，Dy表示y对t的一阶导数； Dny表示y对t的n阶导数。 * * 例如： 求 的解。可键入： [x,y]=dsolve(Dx=y,Dy=-x) 得x = cos(t)*C1+sin(t)*C2 y = -sin(t)*C1+cos(t)*C2 * * 例如求 , , f(0)=0 , g(0)=1 的解。 可输入指令： P=Df=3*f+4*g,Dg=-4*f+3*g; v=f(0)=0,g(0)=1; [f,g]=dsolve(P,v) f = exp(3*t)*sin(4*t) g = exp(3*t)*cos(4*t) 注意：微分方程表达式中字母D必须大写 *