MATLAB与数值分析第一次实验报告
- 格式:doc
- 大小:176.00 KB
- 文档页数:6
MATLAB与数值分析
实验报告
报告人:朱泽智
学号:2015020904033
时间:第五周周二7:30
一、 实验目的
通过上机编程实验,使我们熟悉对MATLAB 软件平台的使用,使学生掌握MATLAB 的编程技巧.
二、 实验原理
1. 给定一个n 维向量[]
121,,,,n n x x x x -= x ,编写MATLAB 程序,实现以下编程实现以下计算,用户能输入不同的初值n 以及系数i 。 121111i i i i n n x x x x y x x x x -++-++++=⋅⋅
2. 编程实现奥运5环图,允许用户输入环的直径。
3. 绘制以下二元函数的三维图像,采用mesh 、meshc 和surf 等三维函数命令
[][]10,10,10,10z x y =∈-∈-
4. 使用一个以上函数(function )的方式实现一个动画图形显示:使蓝色的小球沿着正弦曲线运动。
根据matlab 的程序设计对其进行程序编辑实现上面的要求。
三、 实验方案
1. clear all
close all
clc
n=input;
x=input;
i1=input;
sum=0;
for ii=1:i1
sum=sum+x(ii);
end
plu=1;
for jj=(i1+1):n
plu=plu*x(jj);
end
y=sum/plu;
fprintf
2.figure;
grid on;
hold on;
axis equal;
r=input('key in the r\n');
d=r/3;
p=pi/2:0.01:7*pi/4;
y=r*sin(p);x=r*cos(p);
plot(x,y,'b','linewidth',6);
p=pi/4:0.01:5*pi/4;
y=r*sin(p);x=r*cos(p);
plot(x+r+r+2*d,y,'k','linewidth',6);
p=-pi/2:0.01:5*pi/4;
y=r*sin(p);x=r*cos(p);
plot(x+4*d+4*r,y,'r','linewidth',6);
p=0:0.01:2*pi;
y=r*sin(p);x=r*cos(p);
plot(x+r+1*d,y-r,'y','linewidth',6);
p=5*pi/4:0.01:7*pi/4;
y=r*sin(p);x=r*cos(p);
plot(x+r+r+2*d,y,'k','linewidth',6);
p=0:0.01:2*pi;
y=r*sin(p);x=r*cos(p);
plot(x+3*r+3*d,y-r,'g','linewidth',6);
p=-pi/4:0.01:pi/2;
y=r*sin(p);x=r*cos(p);
plot(x,y,'b','linewidth',6);
p=-pi/4:0.01:pi/4;
y=r*sin(p);x=r*cos(p);
plot(x+r+r+2*d,y,'k','linewidth',6);
p=5*pi/4:0.01:3*pi/2;
y=r*sin(p);x=r*cos(p);
plot(x+4*d+4*r,y,'r','linewidth',6);
3. clear all;
close all;
clc;
x=-10:0.1:10;
y=-10:0.1:10;
[X,Y]=meshgrid(x,y);
u=sqrt(X.^2+Y.^2);
z=sin(u)./u;
figure(1)
meshc(z)
figure(2)
mesh(z)
figure(3)
surf(z)
4.
t=0:pi/500:10*pi;
y=sin(t);
plot(t,y,'m');
h=line('color',[0,0,1],'marker','.','markersize',50); for i=1:length(t);
set(h,'xdata',t(i),'ydata',y(i));
M(i)=getframe;
end
movie(M);
四、实验结果
1.
2.
3.
4.
五、分析与讨论
自己设计程序并运行后发现很多的问题,首先在第三个问题中开始并
没有做到让三个图同时出现,在老师的指导下做到了用figure来实现
多个图同时出现的情况。其次在第四个问题中没有做到让小球运行的
轨迹显示出来。还有待完善。总体来说,本次实验较为成功的完成了
给出的四个问题。