北京交通大学matlab复习资料

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)

4.求解方程x5+6x4-3x2=10的5个根，并将其位置用五角星符号标记在复平面上，要求横纵坐标轴刻度等长，注明虚轴和实轴，在title位置上写出方程。

clc;closeall;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)');

6.某班同学成绩已经存放在矩阵A中，每行为某一位同学的数据，第1列为学号，第2列至第4列为其三门课程的成绩，试编程按照3门课平均成绩由小到大的顺序重新排列成绩表，并放在矩阵B中。

clc;clearall;close all;

A=[99234012,95,73,88;

99234033,84,77,80;

99234009,66,80,72;

99234067,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;clearall;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,i);

x=0:0.01:1;

y=0:0.02:2;

[x,y]=meshgrid(x,y);

z=1./sqrt((1-x.^2)+y.*2)+1./sqrt((1+x.^2)+y.^2);

mesh(x,y,z)

view(angle{i});title(ttl{i});

end

6. 用鼠标左键在图形窗口上取5个点，在每个点的位置处写出一个字符串来显式鼠标点的横坐标值，然后将这些点连成折线。

axis([0,5,0,5]);

hold on;

box on;

x=[];y=[];

gtext('');