计算物理课业.docVIP

Matlab上机程序代码及结果 MATLAB基础 1、A=reshape(1:12,3,4),r=[3 2 1]将数组A的第i列减去r(i),如何实现？ 程序代码： A=reshape(1:12,3,4) r=[3 2 1] for i=1:3 B(:,i)= A(:,i)-r(i); End disp(B=); disp(B); 运行结果： shuzu A = 1 4 7 10 2 5 8 11 3 6 9 12 r = 3 2 1 B= -2 2 6 2 -1 3 7 5 0 4 8 8 2、求解线性方程组AX=B 　　 1 1.5 2 9 7 3 　　 0 3.6 0.5 -4 4 -4 其中A= 7 10 -3 22 33 ， B= 20 　 3 7 8.5 21 6 5 　 3 8 0 90 -20 16 程序代码：A=[1,1.5,2,9,7;0,3.6,0.5,-4,4;7,10,-3,22,33;3,7,8.5,21,6;3,8,0,90,-20] B=[3;-4;20;5;16] X=A\B 运行结果： jiefc A = 1.0000 1.5000 2.0000 9.0000 7.0000 0 3.6000 0.5000 -4.0000 4.0000 7.0000 10.0000 -3.0000 22.0000 33.0000 3.0000 7.0000 8.5000 21.0000 6.0000 3.0000 8.0000 0 90.0000 -20.0000 B = 3 -4 20 5 16 X = 3.5653 -0.9255 -0.2695 0.1435 0.0101 3、计算函数值 x (x1) 2x-1 (1x10) 3x-11 (x=10) 程序代码： x=input(请输入一个数值：) if x1 y=x; elseif x1 x10 y=2*x-1; else y=3*x-11; end disp(函数值为：); disp(y); 运行结果： hsqiujie 请输入一个数值：11 x = 11 函数值为： 22 4、求一元二次方程ax2 +bx+c=0 程序代码： a=input(a=); b=input(b=); c=input(c=); d=b*b-4*a*c; if d0 disp(不存在实数根); else x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)]; disp([x1=,num2str(x(1)),,x2=,num2str(x(2))]); end 运行结果： yyecfcqiujie a=3 b=4 c=2 不存在实数根 yyecfcqiujie a=5 b=8 c=2 x1=-0.3101,x2=-1.2899 5、输入三角形的三条边，求面积 程序代码： a=input(a=); b=input(b=); c=input(c=); cos(c)=(a^2+b^2-c^2)/(2*a*b); sin(c)=sqrt(1-cos(c)^2); s=a*b*sin(c)/2; disp([三角形面积为：,num2str(s)]); 运行结果： jisuanmj a=3 b=4 c=5 三角形面积为：6 6、产生两个随机n阶方阵和一个n阶列向量，将两个n阶方阵进行四则运算，利用克拉莫法则将第1个方阵和列向量所构成的方程组求解 程序代码： n=input(请输入一个整数); A=randn(n) B=randn(n) C=randn(n,1) jia=A+B jian=A-B