7. subs (f,’x’,a):将a赋值给符号表达式f中的x，并计算出值。
【实验环境】
MatlabR2010b
二、实验内容
实验1设 ， ，分别绘出 的图形.
1.分析问题
本题很明显先利用for语句求出函数在个点的值，再通过plot函数分别绘出图形。
2.问题求解
x=-4:0.01:4;
gx=exp(x);
3.结果
4.结论及分析
通过实验，结论正确，证明分析无误。
三、实验小结
通过本次matlab实验我发现matlab并不是别人说的那么无聊，你可以在实验中找到很多乐趣，如用plot函数、surf函数绘制各种奇异的图形，并且你还会在完成试验后有一种成就感。
g2 =
(13^(1/2)/12 - x + 5/12)*((5*13^(1/2))/6 - 12*(13^(1/2)/12 + 5/12)^2 + 19/6) + 13^(1/2)/12 - 5*(13^(1/2)
4.结论及分析
通过实验，结论正确，证明分析无误。
实验3判断平面 、球 、椭球
之间的位置关系.
⑵求 ，并求出 的点 ；
⑶分别求出过 ，斜率分别为 、 的两条直线，并作出题目的图像；
⑷把上述图像在同一个坐标系下显示。
1.分析问题:
本题是根据给出的函数方程绘制图形，再用diff函数求方程的倒数，求出满足方程 的x1、x2再求出在x1、x2的y1、y2绘制经过他们的直线。
2.问题求解:
symsx
fx=4*x.^3-5*x.^2+x-2;
f=diff(fx)
ezplot(fx,[0,1])
b=solve(f-(subs(fx,'x',1)-subs(fx,'x',0))/(1-0));
x1=b(1)
x2=b(2)
y1=subs(fx,'x',x1);
y2=subs(fx,'x',x2);
f1=subs(f,'x',x1)
f2=subs(f,'x',x2)
数学与信息科学系实验报告
实验名称Matlab作图
所属课程数学软件与实验
实验类型综合型实验
专业信息与计算科学
班级
学号
姓名
指导教师
一、实验概述
【实验目的】掌握matlab环境下怎样绘图，了解用plot函数、surf函数绘制二维图、三维图，用diff函数求导数以及用solve函数求解，进一步巩固for语句和abs函数。
g1=f1*(x-x1)+y1
holdon
ezplot(g1,[0,1])
g2=f2*(x-x2)+y2
holdon
ezplot(g2,[0 1])
axis([0 1 -2.4 -1.9])
holdoff
3.结果
f =
12*x^2 - 10*x + 1
x1 =
5/12 - 13^(1/2)/12
x2 =
1.分析问题:
本题可以根据三个函数方程作出三个图行来判断三个图行间的位置关系。
2.问题求解:
x=-60:1:60;
y=-60:1:60;
[x,y]=meshgrid(x,y);
z=1/3*(2*x-y+11);
surf(x,y,z)
holdon
surf(x2,y2,z2)
[u,v]=meshgrid(0:0.1:2*pi);
x1=60*cos(u).*sin(v)+3;
y1=60*sin(u).*sin(v)-5;
z1=60*cos(v)-2;
surf(x1,y1,z1)
holdon
x2=40*cos(u).*sin(v)-30;
y2=100*sin(u).*sin(v)+20;
z2=10*cos(v)+100;
surf(x2,y2,fgx(i)=0;
end
end
plot(x,fx,':r')
holdon
plot(x,gx,'b')
holdon
plot(gx,fgx,'--y')
axis([-6,10,-5,10])
3.结果
4.结论及分析
通过实验，结论正确，证明分析无误。
实验2已知函数
⑴做出 的图像；
fx=[];
fgx=[];
fori=1:length(x)
ifabs(x(i))<1
fx(i)=1;
elseifabs(x(i))>1
fx(i)=-1;
else
fx(i)=0;
end
end
fori=1:length(gx)
ifabs(gx(i))<1
fgx(i)=1;
elseifabs(gx(i))>1
13^(1/2)/12 + 5/12
f1 =
(5*13^(1/2))/6 + 12*(13^(1/2)/12 - 5/12)^2 - 19/6
f2 =
12*(13^(1/2)/12 + 5/12)^2 - (5*13^(1/2))/6 - 19/6
g1 =
(x + 13^(1/2)/12 - 5/12)*((5*13^(1/2))/6 + 12*(13^(1/2)/12 - 5/12)^2 - 19/6) - 5*(13^(1/2)/12 - 5/12)^2 - 4*(13^(1/2)/12 - 5/12)^3 - 13^(1/2)/12 - 19/12
【实验原理】
1.for循环：for循环控制变量=表达式1：表达式2：表达:3
语句
end
2. Plot(x,y):求二维图形。
3.solve(‘表达式’)：求表达式的解。
4. diff(f):对独立变量求导数，f为符号表达式。
5. abs():求绝对值函数。
6. axis（xmin xmax ymin ymax）:指定x、y的取值范围。