MATLAB与数值分析实验报告一

MATLAB与数值分析实验报告

报告人：秦旸照

学号： 2015020901033 时间： 2016.4.8

电子科技大学电子工程学院

一、实验目的

实验一：MATLAB软件平台与程序设计实验

二、实验原理

1.熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序）

2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序）

3. 掌握Matlab函数的编写规范。

4.掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释）

5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。

三、实验方案

1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以

x，y为坐标显示图像

x(n+1) = a*x(n)-b*(y(n)-x(n)^2);

y(n+1) = b*x(n)+a*(y(n)-x(n)^2)

2. 编程实现奥运5环图，允许用户输入环的直径。

3. 实现对输入任意长度向量元素的冒泡排序的升序排列。

不允许使用sort函数。

四、实验结果

1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以

x，y为坐标显示图像

x(n+1) = a*x(n)-b*(y(n)-x(n)^2);

y(n+1) = b*x(n)+a*(y(n)-x(n)^2)

clc

clear all

n=input('please input n');

i=input('please input i');

x=10*rand(1,n);

x1=0;

x2=1;

for j=1:i

x1=x(j)+x1;

end

for j=i+1:n

x2=x2*x(j)

end

y=x1/x2

2. 编程实现奥运5环图，允许用户输入环的直径。

clear;

clc;

t=0:0.01:2*pi;

d=input('d=');

r=d/2

x=2*r*cos(t);

y=2*r*sin(t);

plot(x,y,'b');

hold on;

plot(x+3*r,y,'k');

plot(x+6*r,y,'r');

plot(x+1.5*r,y-2.5*r,'y');

plot(x+4.5*r,y-2.5*r,'g');

plot off;

axis equal;

3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort函数。

clc

clear

x=-10.0:0.1:10.0;

y=x'

u=ones(size(y))*x;

v=y*ones(size(x));

r=sqrt(u.^2+v.^2)+eps;

z=sin(r)./r;

mesh(z)

figure;meshc(z)

figure;surf(z)

shading interp

xlabel('x');ylabel('y'); zlabel('z');