资源描述：
《matlab上机程序设计》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、*******************实践教学*******************兰州理工大学理学院2015年秋季学期Matlab上机程序设计专业班级：信息与计算科学（2）班姓名：学号：13540216指导教师：成绩：2015matlab上机题目练习二-----------科学计算1.对于，如果，，求解X，并求A的行列式、转置、逆矩阵、秩和条件数等运算A=[492;764;357];B=[37;26;28];X=ABC=det(A)R=rank(A)A'D=inv(A)A1=norm(A,1)A2=norm(A,2)A3=norm(
2、A,inf)A4=cond(A,1)A5=cond(A,2)A6=cond(A,inf)X=-0.51184.04271.3318C=-211.0000R=3ans=473965247D=-0.10430.2512-0.11370.1754-0.10430.0095-0.0806-0.03320.1848A1=20A2=16.0018A3=17A4=7.7725A5=5.4858A6=7.97631.，，观察a与b之间的六种关系运算的结果a=[125;36-4];b=[8-74;362];r1=abr3=a<=br4=a>
3、=br5=a==br6=a~=br1=100001r2=011000r3=100111r4=011110r5=000110r6=1110011.对多项式进行因式分解，并求的值symsx;f=x^4-5*x^3+5*x^2+5*x-6;factor(f)subs(f,2.5)ans=(x-1)*(x-2)*(x-3)*(x+1)ans=-21/164.，用符号微分求A=sym('[ax^21/x;exp(a*x)log(x)sin(x)]')diff(A)ans=[0,2*x,-1/x^2][a*exp(a*x),1/x,cos(x)]5
4、.求和symsx;f1=x^2*atan(x)f2=x-x^2int(f1,x)int(f2,x,0,1)ans=log(x^2+1)/6+(x^3*atan(x))/3-x^2/6ans=1/66.微分方程的通解.y=dsolve('D2y-2*Dy+5*y=exp(x)*cos(2*x)','x')y=(3*cos(2*x)*exp(x))/32+(cos(6*x)*exp(x))/32+sin(2*x)*exp(x)*(x/4+sin(4*x)/16)+C2*cos(2*x)*exp(x)+C3*sin(2*x)*exp(x)7.
5、求代数方程组关于x,y的解[x,y]=solve('a*x^2+b*y+c=0','x+y=0','x','y')x=(b+(b^2-4*a*c)^(1/2))/(2*a)(b-(b^2-4*a*c)^(1/2))/(2*a)y=-(b+(b^2-4*a*c)^(1/2))/(2*a)-(b-(b^2-4*a*c)^(1/2))/(2*a)8.求函数的一阶导数.并求symsabxyf=sin(a*x)*cos(b*x);d=diff(f,x)subs(d,1/(a+b))d=a*cos(a*x)*cos(b*x)-b*sin(a*x)*
6、sin(b*x)ans=a*cos(a/(a+b))*cos(b/(a+b))-b*sin(a/(a+b))*sin(b/(a+b))9.已知,计算和并在区间上作出函数的图形。symsx;f=(x^3-9*x)/(x^3-x);a1=limit(f,x,inf)a2=limit(f,x,1)x=-4:0.01:4;f=(x.^3-9.*x)./(x.^3-x);plot(x,f)a1=1a2=NaN练习三--------图形绘制1.符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t的变化范围为
7、[0,2p]ezplot('sin(3*t)*cos(t)','sin(3*t)*sin(t)',[0,2*pi])2.有一组测量数据满足，t的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线，并加入标题、X轴和Y轴的标记、网格显示和图例（用代码形式生成）；再用子图画出三条曲线，对每个子图命名如：“a=0.1”。t=0:0.1:10;y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5*t);plot(t,y1,'b*',t,y2,'r+',t,y3,'y.-'
8、);title('y=exp(-a*t),t从0到10');legend('y1','y2','y3');xlabel('x');ylabel('y');gridonsymsay=exp(-a*t);t=0
