北京交通大学matlab预习复习资料.doc

- P105 1. 主程序： [s1,kk]=quad(f1,0,2); [s2,kk]=quad(f2,2,4); s=s1+s2 函数： function y=f1(x) y=log(x.^2); function y=f2(x) y=log(16)./(2+sin((x+1).*pi)); f1=inline(log(x.^2),x); f2=inline(log(16)./(2+sin((x+1).*pi)),x); [s1,kk]=quad(f1,0,2); [s2,kk]=quad(f2,2,4); s=s1+s2 3. syms x a tau s=int(log(x.^2),0.5,tau) solve(s==10) 求解方程x5+6x4-3x2=10的5个根，并将其位置用五角星符号标记在复平面上，要求横纵坐标轴刻度等长，注明虚轴和实轴，在title位置上写出方程。 clc;close all;clear all; factor=[1 6 0 -3 0 -10]; ro=roots(factor); disp(ro); scatter(real(ro),imag(ro),p); axis equal; title(x^5+6x^4-3x^2=10); xlabel(real axis); ylabel(imagine axis(j)); 某班同学成绩已经存放在矩阵A中，每行为某一位同学的数据，第1列为学号，第2列至第4列为其三门课程的成绩，试编程按照3门课平均成绩由小到大的顺序重新排列成绩表，并放在矩阵B中。 clc;clear all;close all; A=95,73,88; 84,77,80; 66,80,72; 92,93,59]; [y,index]=sort(mean(A(:,2:end),2)); B=A(index,:); disp(B); 9. figure (1) [u,v,b]=ginput(5); line(u,v); for i=1:5; str=[(,num2str(u(i)),)]; text(u(i),v(i),str); end 11.解如下微分方程，并绘制出图像y(x). （1）function ydot= f(x,y) ydot=zeros(2,1); ydot(1)=y(2); ydot(2)=2*(y(2).^2); t0=0;tf=1; x0=[0,-1]; [t,x]=ode45(f,[t0,tf],x0); plot(t,x); title(y=f(x)); xlabel(x);ylabel(y); （2）t0=1;tf=2; x0=[1,0]; [t,x]=ode45(f,[t0,tf],x0); plot(t,x); title(y=f(x)); xlabel(x);ylabel(y); function xdot= f(t,x) xdot=zeros(2,1); xdot(1)=x(2); xdot(2)=-1./(x(1).^3); p87 3. 在同一坐标系下绘制t2，- t2，t2sint在t[0,2pi]内的曲线图；选择合适的“曲线线型”，“曲线颜色”、“标记符号”选项，并在图形上加注坐标轴名和图名。 close all;clear all;clc; t=linspace(0,2*pi,1024); y1=t.^2; y2=-t.^2; y3=(t.^2).*sin(t); figure(2); plot(t,y1,:y,t,y2,--g,t,y3,-r); legend(y=t^2,y=-t^2,y=t^2*sin(t)); title(三个函数图形);xlabel(t);ylabel(f(t)); 4. theta=linspace(0,16*pi,2048); subplot(2,2,1);polar(theta,1.0013*theta.^2); subplot(2,2,2);polar(theta,cos(3.5*theta)); subplot(2,2,3);polar(theta,sin(theta)./theta); subplot(2,2,4);polar(theta,1-(cos(7*theta)).^3); 5. ttl={主视图,左视图,俯视图,三维图}; angle={[0,0],[-90,0],[0 90],[-37.5,30]}; for i=1:4 subplot(2,2,