数学实验报告.docVIP

数学实验报告 实验1（定积分） 实验目的 加深对定积分概念的理解。 学会用MATLAB计算定积分。 实验内容 用数值计算和图形展示相结合的方法研究函数的积分和随分割细度的变化趋势。 用数值方法和符号演算计算定积分。 求由抛物线 y=2*x2 与y=x2+3所围成的图形面积A。 解： x=-3:0.1:3; y1=2*x.^2; y2=x.^2+3; plot(x,y1,x,y2) syms x f=(x.^2+3)-2*x.^2; A=int(f,x,-sqrt(3),sqrt(3)) 结果为 A = 4*3^(1/2) 实验2（常微分方程） 实验目的 了解常微分方程的基本概念。 了解常微分方程的解析解。 了解常微分方程的数值解。 学习、掌握MATLAB软件有关的命令。 实验内容 常微分模型的建立及求解。 求解 解： function z=dy1(x,y); z=[y(2);2*y(2)-3*y(1)] y0=[1;0]; [x,y]=ode23s(dy1,[0,30],y0); plot(x,y(:,1),r:,x,y(:,2),k-); xlabel(x) ylabel(y) 结果为 y = C1*exp(3*x)+C2*exp(-x) y = 1/4*exp(3*x)+3/4*exp(-x) 实验3（线性规划） 实验目的 了解线性规划求解的基本方法。 学习、掌握用MATLAB求解线性规划的命令。 实验内容 线性规划的求解。 求解模型 解： c=[6 3 4]; A=[0 1 0]; b=[50]; Aeq=[1 1 1]; beq=[120]; vlb=[30;0;20]; vub=[ ]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 结果为 x = 30.0000 50.0000 40.0000 fval = 490.0000 实验4（非线性规划） 实验目的 了解非线性规划求解的基本方法。 学习、掌握用MATLAB求解非线性规划的命令。 实验内容 非线性规划的求解。 S.t. 解： 化为标准形 S.t. H=[1 0;0 1]; c=[-1,-2]; A=[2 3 ;1 4]; b=[6;5]; aeq=[]; beq=[]; vlb=[0;0]; vub=[]; [x,fval]=quadprog(H,c,A,b,aeq,beq,vlb,vub) 结果为 x = 0.7647 1.0588 fval = -2.0294 实验5（数据统计分析） 实验目的 加深对统计基本概念的理解。 参数估计。 假设检验。 学习、掌握用MATLAB进行参数估计和假设检验的命令。 实验内容 参数估计和假设。 随机抽查某专业60名同学“高等代数”课程的期末考试成绩，成绩如下： 93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55 试统计该专业学生“高等代数”课程成绩的分布规律。 解： x=[|93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55] T=[mean(x),median(x),var(x),std(x),kurtosis(x),skewness(x)] T = 80.1000 80.5000 94.2949 9.7106 3.1529 -0.4682 该专业学生“高等代数”课程成绩的平均分为80.1分，中位数为80.5分，方差为94.2949，标准差为9.7106，峰度为3.1529，偏度系数为-0.4682。 hist(x,5) normplot(x) [muhat,sigmahat,muci,sigmaci]=normfit(x) muhat = 80.1000 sigmahat = 9.7106 muci = 77.5915 82.6085 sigmaci = 表明该专业学生“高等代数”课程成绩的平均分为80.1分，方差为9.7106，均值的0.95置信区