matlab在积分中的应用报告第五组.docVIP

MATLAB在积分中的应用报告第五组 ———————————————————————————————— 作者： ———————————————————————————————— 日期： 2021—2021学年第1学期 MATLAB在积分中的应用报告 课程名称： 工程应用数学B 案例名称： MATLAB在积分中的应用 专业班级： 11级自动化 起止时间： 第12-15周 组 别： 第五组 指导教师： 成 绩： 一. 小组成员〔具体分工〕 姓名 学号 具体分工 实验步骤，实验程序 实验结果，实验总结 提炼题目，问题分析 二. 案例描述 1、电路有关： RC串联电路，开关S闭合前电路处于零初始状态。在T=0时刻，开关S闭合，电流接入直流电压源，求电阻消耗的电能为多少： 2、数学有关： ，其中。 3、物理学有关： 有一外半径R1为10cm，内半径R2为7cm的金属球壳，在球壳中放一半径R3为5cm的同心金属球壳。假设使球壳和球均带有的正电荷，问球心的电势为多少？ 三. 提炼问题，问题分析〔将实际问题数学化〕 RC电路的零状态响应，开关S闭合前Uc〔0-〕=0，t=0时刻，开关S闭合，电路接入直流电压源Us。根据KVL有，，。 ，。 简单的三重积分问题，首先可以进展笔算先单后重，得到Ｉ。接着可以用MATLAB进展验证。 　 为了计算球心的电势，必须先计算出各点的电场强度，由于在所讨论的范围内，电场具有求对称性，因此可以用高斯定理计算各点的电场强度。 从球内开场：取的球面为高斯面，由导体的静电平衡条件，球内的电场强度为。 　　　　　　当做一个高斯面，在此高斯面内的电荷仅是半径为的球上的电荷，由高斯定理。 　　　　　对于的球面上的各点，。 　　　　　当时，。 再由电势的定义 = 四. 解决方案〔提供多种备选〕 1.电路有关问题解决方案： ①根据提炼的数学问题，进展笔算求解一阶非齐次方程，得到通过电阻的电流，再通过积分求解电阻消耗的电能W。 ，。 ==。在此电路中，无论电路中电容C与电阻R的数值是多少，在充电过程中，电源提供的电能只有一半转变成电场能量储存于电容中，另一半那么为电阻所消耗，即电容的效率只有50%。 ②应用MATLAB求解一阶非齐次方程，得到通过电阻的电流，再用MATLAB求解积分电阻消耗的电能W。 2.数学有关问题解决方案： ①简单的三重积分问题，首先可以进展笔算先单后重，得到I。 ②应用MATLAB编程，求解三重积分问题。可将三重积分问题首先根据题目要求转化为三个定积分，最后再求定积分的值。 3.物理有关问题解决方案： ①直接运用物理公式抽象出数学公式，进展数学运算。 解决思路与过程见提炼问题，问题分析。 ②使用MATLAB计算，仍然是简单积分问题。编写符号积分函数int，问题最终还是转化为定积分的计算。 方案实施及结果〔实验程序等形式〕 1.电路有关程序代码： syms t Us R C W=int(((Us/R)*exp(-t/(R*C)))^2*R,t,0,inf) 根据R、C的取值范围〔均为正数〕可得W= 2.数学有关问题解决方案： syms x y z F2=int(int(int(z*sin(x+y+z),z,0,(pi/2-x-y)),y,0,(pi/2-x)),x,0,pi/2) VF2=vpa(F2) VF2 = pi^2/8 – 1 x=-pi:0.1:pi; y=-pi:0.1:pi; z=-pi:0.1:pi; [X,Y]=meshgrid(x,y); Z=X+Y-pi/2; surf(X,Y,Z);colormap(cool) 平面区域 x+y+z=pi/2 3.物理学有关问题程序代码： syms a r q a=8.85*10^-12; q=10^-8; Vo=int(q/(4*pi*a*r*r),r,0.05,0.07); Vo = 2417851639229258349412352/4705663114835422578125 Vo2=vpa(Vo) Vo2 = 513.81741111184944956426766213 syms a r q a=8.85*10^-12; q=10^-8; Vo=int(2*q/(4*pi*a*r*r),r,0.1,inf); Vo = 15474250491067253760/8604641124270487 ] Vo2=vpa(Vo) Vo2 = 1798.36