《信息学院信息与计算科学专业实验报告格式(ver2011)(1)》由会员分享，可在线阅读，更多相关《信息学院信息与计算科学专业实验报告格式(ver2011)(1)（15页珍藏版）》请在装配图网上搜索。
1、评分签名日期湖南商学院实验报告 课程名称 MATLAB科学计算编程语言 实验名称 MATLAB运算基础 专业班级 信科1021 姓 名 王佩 学 号 100322040 实验日期 第十周星期一 1-2节 实验地点 E601 20112012学年度第 一 学期一、 实验目的1. 掌握数据统计和分析的方法。2. 掌握数值插值与曲线拟合的方法及其应用。3. 掌握多项是的常用运算。4. 掌握求数值导数和数值积分的方法。5. 掌握代数方程数值求解的方法。6. 掌握常微分方程数值求解的方法。二、实验环境 1、硬件配置：Pentium（R） Dual-Core CUP E6500 2.93GHz，1.96的
2、内存 2、软件环境：系统版本：Microsoft Windows XP Professional Service Pack 3 软件版本：MATLAB 7.0三、实验内容 (包括本实验要完成的实验问题及需要的相关知识简单概述)1.利用MATLAB提供的rand函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：（1） 均值和标准方差（2） 最大元素和最小元素（3） 大于0.5的随机数个数占总数的百分比2.将100个学生5门功课的成绩存入距阵P中，进行如下处理： （！） 分别求每门课的最高分，最低分及相应学生序号。（2） 分别求每门课的平均分和标准方差。（3） 5门课总分的最高分、最
3、低分及相应的学生序号。（4） 将5门课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh。3.有3个多项式P1(x)=x4+2*x3+4*x2+5，P2(x)=x+2，P3(x)=x2+2*x+3，试进行下列操作：（1） 求P（x）=P1(x)+P2(x)*P3(x)。（2） 求P(x)的根。（3） 当x取矩阵A的每一元素时，求P（x）的值。其中 A=-1，1.2，-1.4；0.75，2，3.5；0，5，2.5（4） 当以矩阵A为自变量时，求P（x）的值。其中A的值与（3）题相同。4. 求函数在指定点的数值导。F(x)=x，x2，x3；1，2*x，3*x2；0，6，2*x，x=1，2，3
4、。5. 用数值方法求定积分。（1） I1=quad(sqrt（cos(t.2)+4*sin(2*t).2+1)，0,2*pi)的近似值。（2） I2=quad（log（1+x）./(1+x.*x)，0，1）。四实验步骤 （包括实验原理、数据的准备、运行过程分析、源程序（代码）、图形图象界面等）注：本项可以增加页数1.(1) A=rand(1,30000);b=mean(A)std(A,0,2)b = 0.5006ans =0.2888(2) max(A)min(A)ans = 0.9999ans = 5.2801e-006(3) n=0;for i=1:30000if A(i)0.5n=n+1
5、;endendp=n/30000p =0.50122.(1) A=45+51*rand(100,5);Y,U=max(A)a,b=min(A)Y = 95.5761 95.5995 94.5836 95.0827 95.5736U = 44 19 25 54 90a = 45.9833 45.0408 45.6024 45.0365 45.1063b =97 74 69 62 73(2) m=mean(A)s=std(A)m = 68.0703 67.5230 70.1698 68.9722 69.6234s = 14.8390 13.9109 13.6724 14.1876 14.9388(
6、3) sum(A,2)Y,U=max(ans)a,b=min(ans)ans = 377.4243 342.4960 319.9122 369.0263 322.3655 316.2726 350.3343 380.6332 359.8834 354.7636 386.6063 343.1142 305.1376 341.3232 370.2835 364.6223 366.8240 362.2434 372.3862 351.1440 358.7508 347.2765 397.4010 300.2541 348.5861 325.2011 407.4349 353.7849 374.218
7、0 357.5508 342.3605 356.6098 338.6148 369.4695 345.5606 366.7071 365.1810 307.1481 347.4952 352.1757 323.1823 373.4069 360.2183 357.3025 300.4217 378.3015 303.2919 330.7524 316.6786 344.6667 312.2993 312.8430 266.3132 372.2004 308.3939 344.5447 356.7960 342.6655 354.1309 388.3488 374.1251 335.7883 3
8、53.0814 268.9134 350.9218 418.4755 346.0711 286.8756 299.9256 312.8021 365.9626 324.4192 287.2028 306.8935 291.7168 308.0746 420.1381 363.2507 379.3287 263.7036 357.3413 381.1539 342.9214 342.8388 329.9134 376.8161 455.4518 373.6933 320.4081 329.0542 335.6717 347.8795 349.3276 371.7100 329.9586 304.
9、8626 318.8759 285.7668 320.8355 338.3844Y = 455.4518U = 87a = 263.7036b =80(4) zcj,xsxh=sort(ans)zcj = 263.7036 266.3132 268.9134 285.7668 286.8756 287.2028 291.7168 299.9256 300.2541 300.4217 303.2919 304.8626 305.1376 306.8935 307.1481 308.0746 308.3939 312.2993 312.8021 312.8430 316.2726 316.6786
10、 318.8759 319.9122 320.4081 320.8355 322.3655 323.1823 324.4192 325.2011 329.0542 329.9134 329.9586 330.7524 335.6717 335.7883 338.3844 338.6148 341.3232 342.3605 342.4960 342.6655 342.8388 342.9214 343.1142 344.5447 344.6667 345.5606 346.0711 347.2765 347.4952 347.8795 348.5861 349.3276 350.3343 35
11、0.9218 351.1440 352.1757 353.0814 353.7849 354.1309 354.7636 356.6098 356.7960 357.3025 357.3413 357.5508 358.7508 359.8834 360.2183 362.2434 363.2507 364.6223 365.1810 365.9626 366.7071 366.8240 369.0263 369.4695 370.2835 371.7100 372.2004 372.3862 373.4069 373.6933 374.1251 374.2180 376.8161 377.4
12、243 378.3015 379.3287 380.6332 381.1539 386.6063 388.3488 397.4010 407.4349 418.4755 420.1381 455.4518xsxh = 80 53 64 98 68 73 75 69 24 45 47 96 13 74 38 76 55 51 70 52 6 49 97 3 89 99 5 41 72 26 90 85 95 48 91 62 100 33 14 31 2 58 84 83 12 56 50 35 67 22 39 92 25 93 7 65 20 40 63 28 59 10 32 57 44
13、81 30 21 9 43 18 78 16 37 71 36 17 4 34 15 94 54 19 42 88 61 29 86 1 46 79 8 82 11 60 23 27 66 77 873.(1) p1=1,2,4,0,5;p2=1,2;p3=1,2,3;p=p1+0,conv(p2,p3) p = 1 3 8 7 11(2) A=roots(p)A = -1.3840 + 1.8317i -1.3840 - 1.8317i -0.1160 + 1.4400i -0.1160 - 1.4400i(3) A=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5;polyva
14、l(p,A)ans = 1.0e+003 * 0.0100 0.0382 0.0125 0.0223 0.0970 0.41220.0110 1.2460 0.1644(4) polyvalm(p,A)ans = 1.0e+003 * 0.0076 -0.1281 -0.0775 0.1328 1.3900 1.16440.1824 1.7364 1.51984. f(x)=(x,x2,x3),(1,2x,3x2),(0,2,6x)X=1,2,3? f(x)=(x,x2,x3),(1,2x,3x2),(0,2,6x)Error: Incomplete or misformed expressi
15、on or statement.5.（1）I1=quad(sqrt(t.2)+4*sin(2*t),0,2*pi)I1 = 19.7392（2）I2=quad(log(1+x)./(1+x.*x),0,1)? Error using = inlineevalError in inline expression = log(1+x)./(1+x.*x)? Error: Unbalanced or misused parentheses or brackets.Error in = inline.subsref at 25 INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);Error in = quad at 62y = f(x, varargin:); I2=quad(log(1+x)./(1+x.*x),0,1)I2 = 0.2722五、实验小结： （包括主要实验问题的最终结果描述、详细的收获体会，待解决的问题等） 本次实验内容较多，遇到的问题也就自然而然的变多了，但每一个问题都会有它自己的解决办法，通过一段时间的钻研，问题也都迎刃而解。总的来说，我们还是应该多实践，更好的把握每一个实验。 手写签名：
信息学院信息与计算科学专业实验报告格式(ver2011)(1)
