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

《数学实验》报告 学 院： 电子与信息学院 专业班级： 通信工程4班 学 号： 201130301443 姓 名： 李腾辉 实验名称： Matlab基础知识 实验日期： 2013.02.28 Matlab基础知识 目的： 熟悉MATLAB的具体操作与操作键。 掌握MATLAB中的常用函数与变量、表达式的定义方法。 熟悉MATLAB常用的工作方式M文件的编程工作方式 掌握MATLAB语言中的程序结构。 任务 建立一个M文件，求所有的水仙花数。所谓的水仙花数是指一个三位数，其各位数字的立方和等于该数本身。例如153是一个水仙花数，因为153=13+53+33 。 用subplot分别在不同的坐标系下画出下列四条曲线，为每幅图形加上标题： 概率曲线 四叶玫瑰曲线 叶形线 正弦曲线 实验过程 建立一个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 运行结果： 用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(正弦曲线); 运行结果：