北京交通大学matlab复习资料

北京交通大学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;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)');

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

clc;clear all;close all;

A=[99234012,95,73,88;

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,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('');

for i=1:5

[x1,y1,button]=ginput(1);

if(button~=1)

break;

end

plot(x1,y1,'o');

x=[x,x1];y=[y,y1];

text(x1,y1,num2str(x1));

line(x,y);

end

hold off

老师： 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

油桶

clc;

clear;

syms x y h;

w=ezplot('x^2/(2.5^2)+(y-2)^2/4-1=0'); %绘制椭圆.

set(w,'Color','k');

hold on;

axis([-2.55,2.55,0,4.05]); %设置坐标值.

set(gca,'XTick',[],'XColor','w'); %%%将坐标隐去.

set(gca,'YTick',[],'YColor','w');

line([0,0],[0,4],'Color','k','LineStyle','-') %画出中间刻度线.

hold on;

s=int(((100-25*(y-2).^2)/16).^0.5,0,h); %用int函数求不定积分.

for i=1:25

H(i)=solve(s-0.3125*i,h); %求得各个刻度的高度.