资源描述：
《matlab实验(西北农林科技大学)》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、matlab实验(西北农林科技大学)课程实验报告学年学期２０11５--62021学年第22学期课程名称MATLAＢ与科学计算实验名称实验一MATLAＢ环境与数值计算实验室测量测绘实验中心计算机室专业年级土木313级级学生姓名学生学号22０130提交时间22０１66－０66--02成绩任课教师王增红水利与建筑工程学院实验一MAＴＬAAＢ环境与数值计算1、目得与要求(１）熟练掌握MATＬAB得启动与退出、MATLAB得命令窗口；常用选单与工具栏,通过实例初步认识对MATLAB得命令与程序运行等。(2)
2、熟练掌握ＭATLＡB变量得使用；(３)熟练掌握矩阵得创建、掌握MAＴLＡB得矩阵与数组得运算;(４)熟悉MAＴLAB多项式得运用,掌握ＭATLＡＢ得拟合与插值方法、积分运算与统计运算。2、内容与步骤参见教材实验一、二。3、实验报告提交要求(1)用MALTAＢ命令与M程序两种方式求解三元一次方程组。答:以求解三元一次方程组{2+33a+b5c=54ab+c=9为例:MALTＡＢ命令a=[２-13;31-5;４-11];ｂ＝［5;5;9];x=ａｂｘ=2-１0Ｍ程序a=[2-13;31－5；4-11
3、];b=［5;5；9];x=ab运行结果：x=２-１0(2)用ＭＡLTAＢ命令与M程序两种方式求解33复数矩阵得模与相角。答:MALTＡＢ命令a=ｆix(20＊ｒand(３,3）)ａ=１６185１８121０2119＞b=maｇic(3);A=ａ+ｂＩ；ｒ=ａｂs(A)r=17、８88５18、０２787、810218、２48313、0００012、20６６4、47２19、０5541９、1050ｔheta＝angle(A)1８0/ｐithｅta＝２６、５6513、1７9８50、１944９、4６2322
4、、6１9934、9９2０６3、４3４983、65986、0090MM程序a=ｆix(２0rand(３,3）)b＝mａgiｃ(３）;A=a+bi;r=ａbｓ(Ａ）tｈeta=ａｎgｌe(A）180/pi运行结果：a=15013１９１615１31814ｒ=17、00001、000０14、31７819、2３5416、７63116、5５2913、601520、12４614、14２1ｔheta=28、０72590、000024、77518、9７261７、3540２5、0169１７、1０272６、5６5１8、
5、13０１(3）创建一个４╳４矩阵,求逆矩阵、产生对角矩阵、求矩阵得秩等。＞ａ=fix(2０＊raｎd（4，4）)a=1261010１1３31５9151218０65＞ni=ｉnv(a)nｉ=－0、00５40、0２５1-0、0２940、0664-0、007９0、086３０、0074－0、0536-0、1３７5-0、1１740、１74０-0、07２20、１846０、０507－０、1０290、０４７4＞dｕijiao=ｄｉaｇ(a)duijiao=111１5５zｈi=rａnk(ａ)zｈｉ=4（4)创建两个
6、4╳4矩阵进行加、减、乘法与左除运算。＞＞a=ｆiｘ(1０rand(4,4)+1)ａ=372389２９３１03366９10＞b=fiｘ(10ｒａnｄ(4，4）＋１)b=45682410439367681＞ａ-ｂａｎｓ=a+ｂａｎs=7１28１１10131213６1９６９1３121711ａ＊bａns＝53７91１86711９1４８21６1２1６21００15185133195２03１36＞＞a＼banｓ=14、0000－５、０00０20、000028、0０00－2、13701、64３８-3、0８22-
7、３、73975、4１1０-0、93１56、2466-１２-４-565-85０10-3-10199、2１92-11、287７３、45２１－14、９726－２2、75３4(5)创建一个一元四次多项式,进行多项式求值,求根计算。解:以函数y＝５x4+4x3+3x2+2x+1为例：p=［5,４,3,2,1]ｐ=5432１＞p＝[5,4,3,2,1]；pｏlyvａl（p,3)%计算ｘ=3时多项式得值anｓ=547roots（ｐ)％计算多项式得根ans=０、1378+0、67８2i０、1３７8-０、6７82i
8、-０、5378+0、３5８３i-0、５378-0、3583i(６)创建一个一元三次多项式,取十个不同自变量值，计算多项式得值，并分别进行一阶、二阶拟合,图示原始数据与拟合结果解:＞ｐ＝fｉx（10＊rand(1，４)＋1）ｐ=77８５＞x=1:2:20；y=polyval（ｐ，x);ｐ1=pｏｌyfit(ｘ,y,1);p2＝pｏlｙｆit(x，y，2）；y1=polyｖal(ｐ1,x）;ｙ2＝poｌｙval(p2，ｘ);plｏt(x,y１,"ｒ")hoｌdon＞＞plot(ｘ，y２
