MATLAB实验报告3(1)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
四、实验内容与步骤:
1.绘制下列曲线.
(1) y=x-(x^3)/6
程序输入如下:
fplot('x-(x^3)/6',[-5,5],'r.');
程序输出:
(2) x^2+2*y^2=64
程序输入如下:
ezplot(' x^2+2*y^2-64',[-8,8])
程序输出:
2.设
y=1/(1+exp(-t)) –pi<=t<=pi
在同一图形窗口采用子图的形式绘制条形图阶梯图杆图和对数坐标图等不同图形,并对不同图形加标注说明.
程序输入如下:
t=-pi:pi/10:pi;
y=1./(1+exp(-t));
subplot(2,2,1);bar(t,y,'r');
title('条形图');axis([-4,4,0,1]);
subplot(2,2,2);stairs(t,y,'b');
title('阶梯图');axis([-4,4,0,1]);
subplot(2,2,3);stem(t,y,'g');
title('杆图');axis([-4,4,0,1]);
subplot(2,2,4);semilogx(t,y,'k');
title('对数坐标图');axis([-4,4,0,1]);
程序输出:
3.绘制下列极坐标图.
(1) y=5*cos(x)+4
(2) y=(5*sin(x)*sin(x))/cos(x) (1)程序输入:
x=0:pi/50:2*pi;
y=5*cos(x)+4;
polar(x,y,'-*');
程序输出:
(2)程序输入:
x=-pi/3:pi/50:pi/3;
y=(5.*sin(x).*sin(x))./cos(x);
polar(x,y,'-*');
程序输出:
4.绘制下列三维图形
(1)x=exp(-t/20).*cos(t)
y=exp(-t/20).*sin(t)
z=t
0<=t<=2*pi
(2)z=5
abs(x)<=5
abs(y)<=5
要求应用插值着色处理(1)程序输入:
t=0:pi/10:2*pi;
x=exp(-t/20).*cos(t);
y=exp(-t/20).*sin(t);
z=t;
plot3(x,y,z);
title('三维图形4-1');
xlabel('x');ylabel('y');zlabel('z'); grid on;
程序输出:
(2)程序输入:
[x,y]=meshgrid(-5:0.5:5); z=0*(x-y)+5;
surf(x,y,z);
shading interp;
title('三维图形4-2');
程序输出:
五、实验总结:
2.绘制下列曲线,
(1) y=exp(-x*x/2)/(2*pi)
程序输入:
fplot('exp(-x*x/2)/(2*pi)',[0,5],'r.')
程序输出:
(2) x=t*sin(t)
y=t*cos(t)
程序输入:
t=0:0.1:2*pi;
x=t.*sin(t);
y=t.*cos(t);
plot(x,y);
程序输出:
3.在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点
(1) y=2*x-0.5
(2) x=sin(3*t).*cos(t)
y= sin(3*t).*sin(t)
0<=t<=pi
程序输入:
t=0:pi/100:pi;
x=sin(3*t).*cos(t);
y2=sin(3*t).*sin(t);
y1=2*x-0.5;
plot(x,y1,'m',x,y2,'g');
hold on
k=find(abs(y2-y1)<1e-4);
x1=x(k);
y3=2*x1-0.5;
plot(x1,y3,'bp');
程序输出:
4.分别用plot和fplot函数绘制函数y=sin(1/x)的曲线,分析两曲线的差别程序输入:
x=-1:pi/100:1;
y=sin(1./x);
subplot(2,1,1);plot(x,y,'g');
subplot(2,1,2);fplot('sin(1./x)',[-1,1],'m');
程序输出:
两曲线的差别plot函数在取数据点时一般都是等间隔采样,fplot函数可自适应地对函数进行采样,能更好的反应函数的变化规律
6.绘制曲面图形
(1)x=3*u*sin(v)
y=2*u*cos(v)
z=4*u*u
程序输入:
[u,v]=meshgrid(0:pi/100:2*pi);
x=3*u.*sin(v);
y=2*u.*cos(v);
z=4*u.*u;
mesh(x,y,z);
程序输出:
严重觉得对细节方面很重要,,差一个点就能导致整个程序的不能运行