matlab实验五1308060384聂川

%1.求解系统方程的解 disp('第一题'); disp('求解下面系统方程的解') disp('3a+6b+4c=1') disp('a+5b=2') disp('7b+7c=3') disp('解法一：左除（\）') A=[3,6,4;1,5,0;0,7,7]; B=[1,2,3]'; x=A\B; disp('a、b、c 分别为：'); disp(x) disp('解法二：求逆（inv）') y=inv(A)*B; disp('a、b、c 分别为：') disp(y) %2.求积分 disp('第二题')
第五题结果：
实 验 总 结
通过这次实验我掌握并理解 Maltab 在数值计算的基本用法，还理解了 matlab 在工程领域解题的一般过程。通过不断的测试每个函数和语句的用 法，掌握 Matlab 符合运算的基本方法等。过程很认真，结果，我自己也很 满意。
指 导 教 师 意 见 签名： 年 月 日
环 境
1．线性系统方程：分别使用左除（\）和求逆（inv）求解下面系统方 程的解： 实 验 内 容 较。 要求：1）使用 quad 求积分时，要求使用两种方式要求分别 ①用函数文件撰写被积函数，用 quad 求积分 ②在 quad 函数中使用匿名函数求积分
3a 6b 4c 1 a 5b 2 7b 7c 3
2. 使用 quad 和 trapz 求解 xe x / 3 dx 的数值积分，并与其解析解相比
0
5
提示： 1）这里 trapz 是求解离散积分，因此要先模拟采样，即需先分别对 x 和 y 数据进行采样， 产生相应的数据点， 利用 trapz 对离散采样数据求积分） 2）解析解需利用符号运算进行求解，参考函数 int，并利用 subs 函数 计算求解区间的函数值。 3.求下面函数在[0.5,4]区间内的过零点。（用 fzero 函数）
年利率为 10%。试制定出资金的使用计划，以使 4 年效益之和为最大。 （提示：这是一个有约束条件的最优化问题的求解，需利用线性规划建立函 数方程，假设第一年至第四年使用的资金分别为x1,x2,x3,x4，则可建立以下 方程：）
max : f ( x1 , x2 , x3 , x4 ) x1 x2 x3 x4 x1 400 x (400 x ) 1.1 2 1 x3 ((400 x1 ) 1.1 x2 ) 1.1 x4 (((400 x1 ) 1.1 x2 ) 1.1) x3 ) 1.1
x0=[1;1;1;1]; [x,fval]=fmincon('mycon',x0,A,b)
其中所用的函数：
function f=mycon(x) f=sqrt(x(1))+sqrt(x(2))+sqrt(x(3))+sqrt(x(4)); end
第一题结果：
实 验 结 果
第二题：
第三题：
第四题结果：
提示：首先用函数 ezplot 绘制曲线 f(x)在[0.5,4]上的图像，命令如下 ezplot('x^3-2*x^2*sin(x)+5*x*cos(x)+1/x',[0.5,4]) ；其次，根据 图像得到函数上接近零点的坐标，利用 fzero 函数求其零点。 4. 对边长为 3m 的正方形铁板，在 4 个角剪去相等的正方形以制成方形 无盖水槽，问何种剪法可使水槽的容积最大？（提示：设剪去的正方形边长 为 x，则水槽的体积应为(3-x)2x，利用 fminbnd 求解，考虑实际情况，注意 fminbnd 的求解区间的设置） 5. 设有 400 万元资金，要求 4 年内使用完，若在一年内使用资金 x 万 元，则可得效益 万元（效益不能再使用），当年不用的资金可存入银行，
%第三题，求零点 ezplot('x^3-2*x^2*sin(x)+5*x*cos(x)+1/x',[0.5,4]) [xx,yy]=ginput(2); [x1,y1]=fzero(@fzerol,xx(1)); 其中的 fzerol 函数为： function y=fzerol(x) y=x^3-2*x^2*sin(x)+5*x*cos(x)+1/x end %4.求边长为多少时，水槽的容量最大 x=0:0.1:3; [x,fval]=fminbnd(@myfun1,0.5,1.5) Myfun1 函数为： function y=myfun1(x) y=x.*((3-x).^2); end %5.求效益的最大值 A=[1,0,0,0;1.1,1,0,0;1.21,1.1,1,0;1.331,1.21,1.1,1]; b=[400;400;484;532.4];
实 验 要 求
本次实验涉及参考教材第四章、第五章的学习内容，需要理解相关函数的 含义及调用方法。实验上机过程中会对相关函数进行举例讲解，同学们可在 实验上机时， 边学习边完成实验任务， 也可提前预习相关内容， 完成实验题。 记得多用帮助。
实 验
计算机一台（带有 MATLAB7.0 以上的软件环境）。
贵州大学实验报告
学院：计算机科学与信息学院 姓名 实验时间 实验项目名称 实 验 目 的 聂川 专业： 学号 指导教师 信息安全 班级：131 实验组 成绩
1308060384
实验五 MATLAB 数数值计算的基本用法。 2．理解 matlab 在工程领域解题的一般过程。 3. 掌握 Matlab 符合运算的基本方法。
实 disp('求解 xe^(-x/3)dx 的数值积分，x∈[0,5];并与其解析解相比较') 验 [S,n]=quad(@fesin,0,5); 程
disp('解：S=') disp(S); function f=fesin(x)
序 其中所用函数为：
设 f=x.*exp(-x/3); 计 end