资源描述:
《数学建模实验报告matlab的基本操作编程》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、MATLAB基本操作、编程一、实验目的及意义1.熟悉MATLAB软件的用户环境;2.了解MATLAB软件的一般命令;3.掌握MATLAB向量操作与矩阵运算函数;4.掌握MATLAB软件的基本符号运算命令;5.掌握MATLAB语言编程的循环、条件和选择结构。二、实验内容1.MATLAB软件的向量操作及矩阵运算练习;2.使用MATLAB软件进行符号运算练习;3.用MATLAB语言编写命令M-文件和函数M-文件;三、实验步骤1.在D盘建立一个自己的文件夹;2.开启软件平台——MATLAB,将你建立的文件夹加入到MATLAB的搜索路径中。3.利用帮助了解函数max,
2、min,sum,mean,sort,length,rand,size和diag的功能和用法。4.开启MATLAB编辑窗口,键入你编写的M文件(命令文件或函数文件);5.保存文件(注意将文件存入你自己的文件夹)并运行;6.若出现错误,修改、运行直到输出正确结果;7.写出实验报告,并浅谈学习心得体会。四、实验要求与任务 根据实验内容和步骤,完成以下具体实验,按要求写出实验报告。1.设有分块矩阵,其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵,试通过数值计算验证。2.某零售店有9种商品的单件进价(元)、售价(元)及一周的销量如下表,问哪种商品的利润最大,哪
3、种商品的利润最小;按收入由小到大,列出所有商品及其收入;求这一周该9种商品的总收入和总利润。货号123456789单件进价7.158.253.2010.306.6812.0316.8517.519.30单件售价11.1015.006.0016.259.9018.2520.8024.1515.50销量5681205753580395210415388106943.建立一个命令M-文件:求所有的“水仙花数”,所谓“水仙花数”是指一个三位数,其各位数字的立方和等于该数本身。例如,153是一个水仙花数,因为。4.编写函数M-文件sq.m:用迭代法求的值。平方根迭代公
4、式为迭代的终止条件为前后两次求出的的差的绝对值小于。5.求函数的极限、导数或积分:五.程序代码及运行结果(经调试后正确的源程序)1.设有分块矩阵,其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵,试通过数值计算验证。程序代码:clear;E=eye(3);R=rand(3,2);O=zeros(2,3);X=[3,4];S=diag(X);A=[E,R;O,S];%EQ1,EQ2分别表示待验证的等式左边和右边EQ1=A^2;EQ2=[E,R+R*S;O,S^2];EROSAEQ1EQ2ifEQ1==EQ2fprintf('左边=右边,该等式成立')e
5、lsefprintf('左边不等于右边,该等式不成立')end运行结果:E=100010001R=0.92180.40570.73820.93550.17630.9169O=000000S=3004A=1.0000000.92180.405701.000000.73820.9355001.00000.17630.91690003.0000000004.0000EQ1=1.0000003.68732.028501.000002.95284.6773001.00000.70514.58450009.00000000016.0000EQ2=1.0000003.
6、68732.028501.000002.95284.6773001.00000.70514.58450009.00000000016.0000左边=右边,该等式成立2.某零售店有9种商品的单件进价(元)、售价(元)及一周的销量如下表,问哪种商品的利润最大,哪种商品的利润最小;按收入由小到大,列出所有商品及其收入;求这一周该9种商品的总收入和总利润。货号123456789单件进价7.158.253.2010.306.6812.0316.8517.519.30单件售价11.1015.006.0016.259.9018.2520.8024.1515.50销量568
7、120575358039521041538810694程序代码:clear;inPrice=[7.15,8.25,3.20,10.30,6.68,12.03,16.85,17.51,9.30];%进价outPrice=[11.10,15.00,6.00,16.25,9.90,18.25,20.80,24.15,15.50];%售价sold=[568,1205,753,580,395,2104,1538,810,694];%销量inCome=outPrice.*sold;%收入profit=(outPrice-inPrice).*sold;%利润[max1,p
8、os1]=max(profit);[min1,pos