计算物理matlab简介中的编程.docVIP

例1：计算 》2*(300-2*i)^2+1/sqrt(5+2*i) ans = 1.7999e+005 -2.4001e+003i 例2：解三角方程cos 2x + sin x = 1 》s = solve(cos(2*x)+sin(x)=1) s = pi 0 1/6*pi 5/6*pi 例3：解三次方程x3-2x2 = x-1 》s=solve( x ^ 3 - 2 * x ^ 2 = x - 1) s = 1/6*(28+84*i*3^(1/2))^(1/3)+14/3/(28+84*i*3^(1/2))^(1/3)+2/3 -1/12*(28+84*i*3^(1/2))^(1/3)-7/3/(28+84*i*3^(1/2))^(1/3)+2/3+1/2*i*3^(1/2)*(1/6*(28+84*i*3^(1/2))^(1/3)-14/3/(28+84*i*3^(1/2))^(1/3)) -1/12*(28+84*i*3^(1/2))^(1/3)-7/3/(28+84*i*3^(1/2))^(1/3)+2/3-1/2*i*3^(1/2)*(1/6*(28+84*i*3^(1/2))^(1/3)-14/3/(28+84*i*3^(1/2))^(1/3)) 由于solve将所有的数值都当作符号计算，需简化，用指令vpa（S,D）S表示要替换的对象，D表示要用几位数字。 vpa(s,5) ans = 2.2470+.1e-4*i -.80192+0.*i .55498+0.*i 例4：解方程ax2 + bx + c = 0 syms a b c x f=a*x^2+b*x+c; solve(f) ans = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) solve默认x作为求解的变量，但也可指定变量求解 solve(f,b) ans = -(a*x^2+c)/x 例5：解非线性方程组（课本中的电阻题） syms x y z r1 r2 r3 [x,y,z] = solve(1/x+1/(y+z)-1/r1,1/y+1/(z+x)-1/r2,1/z+1/(x+y)-1/r3) x = -1/2*(r3^2-2*r2*r1-2*r1*r3-2*r2*r3+r2^2+r1^2)/(-r1+r3+r2) y = -1/2*(r3^2-2*r2*r1-2*r1*r3-2*r2*r3+r2^2+r1^2)/(r3+r1-r2) z = 1/2*(r3^2-2*r2*r1-2*r1*r3-2*r2*r3+r2^2+r1^2)/(r3-r1-r2) 例6：解常微分方程 dsolve(Dx=-a*x) ans = C1*exp(-a*t) 例7：画图 x=0:0.1:6; A=[x;4*x]; B=sin(A); plot(x,B(1,:),x,B(2,:)) 将B(1)中所有的负值换成正值画图 a=find(B(1,:)0); B(1,a)=abs(B(1,a)); plot(x,B(1,:),r-.*) 也可将Plot改成comet comet(x,B(1,:)) 例8：运行程序 ballode funtool taylortool demo 例9：数据显示不同 a=23.456; A=23.568974123; format short format long 例10：数值计算误差 1/3连续乘300次，选用不同的数据显示，结果不同 A=0;for k=1:300,A=A+1/3;end format short,A A = 100.0000 format long,A A = 99.99999999999966 将A再连续减1/3共300次，看结果 for k=1:300,A=A-1/3;end A A = -6.994405055138486e-015 B=sqrt(A) B = 0 +8.363255977870393e-008i 例1：向量的生成 a=[3, 5, 7, 8] a=[3 5 7 8] a=[3; 5; 7; 8]