华工数学实验-作业1-Matlab基础知识

《数学实验》报告

学院：电子与信息学院

专业班级：通信工程4班

学号：201130301443

姓名：李腾辉

实验名称：Matlab基础知识

实验日期：2013.02.28

Matlab 基础知识

1． 目的：

- 熟悉MATLAB 的具体操作与操作键。

- 掌握MATLAB 中的常用函数与变量、表达式的定义方法。 - 熟悉MATLAB 常用的工作方式M 文件的编程工作方式

- 掌握MATLAB 语言中的程序结构。

2． 任务

（1） 建立一个M 文件，求所有的水仙花数。所谓的水仙花数是指一

个三位数，其各位数字的立方和等于该数本身。例如153是一个水仙花数，因为153=13+53+33 。

（2） 用subplot 分别在不同的坐标系下画出下列四条曲线，为每幅

图形加上标题：

概率曲线2

x y e -=

四叶玫瑰曲线sin 2ρθ

= 叶形线3

233131t x t t y t ?=??+??=?+?

正弦曲线sin y x

=

3．实验过程

（1）建立一个M文件，求所有的水仙花数。所谓的水仙花数是指一个三位数，其各位数字的立方和等于该数本身。例如153是一个水仙花数，因为153=13+53+33。

思路如下：

取一个数num从100到999之间进行循环递加，循环体中将num 分解成ge,shi,bai三个位数，并且进行判断，如果num ==

ge^3 + shi^3 + bai^3 则输出这个num，否则继续进行循环。

程序如下：

i = 0;

for num = 100:999

ge = mod(num,10);

shi = floor(mod(num,100)/10);

bai = floor(num/100);

if num ==( bai^3 + shi^3 + ge^3 )

i = i + 1;

tips = sprintf('第%d个水仙花数是',i);

disp(tips);

disp(num);

end

end

运行结果：

（2）用subplot分别在不同的坐标系下画出下列四条曲线

思路：

使用 plot和polar进行直角坐标和圆坐标下的绘图

程序：

x1 = -10:0.01:10; %第一幅图概率曲线

y1 = exp(-x1.^2);

%第二幅图四叶玫瑰曲线

theta2 = 0:pi/100:8*pi; % θ的取值范围为0-8∏

r2 = sin(2*theta2);

x3 = @(t)(3*t./(1+t.^3)); %第三幅图叶形线

y3 = @(t)(3*t.^2./(1+t.^3));

x4 = -10:0.01:10; %第四幅图正弦曲线

y4 = sin(x4);

subplot(2,2,1), plot(x1,y1), title('概率曲线');

subplot(2,2,2), polar(theta,r),title('四叶玫瑰曲线');

subplot(2,2,3), plot(x3(-100:0.2:100),

y3(-100:0.2:100)), title('叶形线');

subplot(2,2,4), plot(x4,y4), title('正弦曲线');

运行结果：