有关于Matlab的两道计算题1.某个函数（类似于z=2(x-y^2)^3）在x^2+(y/4)^2=1这个平面曲线上的积分.2.f=x^2+y^2+z^2在函数z=2(x-y^2)^3的空间曲面上的积分.其中x跟y的积分区间均为[-6,6].

来源：学生作业帮助网 编辑：作业帮 时间：2024/10/03 18:44:06
有关于Matlab的两道计算题1.某个函数（类似于z=2(x-y^2)^3）在x^2+(y/4)^2=1这个平面曲线上的积分.2.f=x^2+y^2+z^2在函数z=2(x-y^2)^3的空间曲面上的积分.其中x跟y的积分区间均为[-6,6].
xX[oY+Avұ23Bi!Z +vrc$$& !&jiЧ_>m&> SUUUtkgy_;3glRvqmnVPy8gW}Ϋ{%ws^u/Y~"8{Wm_u~vwveOˤ41ח V}߀>Xi'-2{o$S0?nWvw;]9@!ӉzRdzO$~Otd.'x܌ $C١LLSpL $ u(ρ,̲LbdC8M$MtM1u..a1UTT-YQUEY$bbXX@)sHD*JX2𢊄c銧>Bw}աN6Ʈ Ughd;7Ff2{+\ݻnnYI W[y cEyϽKdcl/L^,/' r6;o[g Ʈ sS7ꢽ_uo?B}$%R{̙y},0<+ ^uJZ%ʙ(Ԗ dg )칻ុܟFUb.̜LitO'9KL$K%3H.,ơgZCNt911YSWpkчE0Yepg:QTDcRBtTUƂ, 1E9ARX@@oCղT"|50??y(=e Qx #EH帣Ѐzk^0gx3Ϻ]0oyjwW+4MuZkE2]'/n^ĜMx,+k>j|/t"@2=I^# M)nRQ~qܕu¨ rgdTػE`Q˵q<e{{^8|7E WI ,$<~ =ax X3 e¹EޙXV_AFjxȠ@K)mA j cX ɮݛw-R4ѐ1{o8μwk KaK vXa ЇS3a=.'ˉl#L)^vu-PRvcocwX"2$^xib-K9*ɢ I8%C0yAºi:LT!Xn3}ї#Lp@4G:dUƺ :qiu8TaAwUW~Bgza8y1mhIϴJ~6? )> 12:  +C*>*ٛuʷtVxe {qcG{dZT1go^= }KJ! p =~TΓEEudw +PIAx3~ai# 8tuE-ܠh \}>6F1SƩnt}_seC::P6ǜ7IfvxM#8~O+@O_[/:&r1ӏY8UA: Kg9Si j7y0%ִ$Q/cA9B"WD 6싾)f7MaVD~9KXrTk(ЧHXo \@MJׁMhgQL7=n 6?h-[df"۲꿋J&feF6To2rt&ܠ.=ؼCBü&#`nÚϭ05-@rSц-(f2Ҝڶ6FF/Oi KXC5:,ryNfuCR0$HdaM2X~2aVa2`Y؃gyMc8"W,~:

有关于Matlab的两道计算题1.某个函数（类似于z=2(x-y^2)^3）在x^2+(y/4)^2=1这个平面曲线上的积分.2.f=x^2+y^2+z^2在函数z=2(x-y^2)^3的空间曲面上的积分.其中x跟y的积分区间均为[-6,6].
有关于Matlab的两道计算题
1.某个函数（类似于z=2(x-y^2)^3）在x^2+(y/4)^2=1这个平面曲线上的积分.2.f=x^2+y^2+z^2在函数z=2(x-y^2)^3的空间曲面上的积分.其中x跟y的积分区间均为[-6,6].

有关于Matlab的两道计算题1.某个函数（类似于z=2(x-y^2)^3）在x^2+(y/4)^2=1这个平面曲线上的积分.2.f=x^2+y^2+z^2在函数z=2(x-y^2)^3的空间曲面上的积分.其中x跟y的积分区间均为[-6,6].
第1题
属于第一类曲线积分,计算原理可以参考附件ppt的92-93页（内容源自薛定宇《高等数学问题的MATLAB求解》,下同）.
使用MATLAB符号积分函数int计算,但无法求出解析解,可通过vpa得到所需精度的数值结果.
代码如下：
syms t
x = cos(t);
y = 4*sin(t);
z = 2*(x-y^2)^3;
F = z * sqrt(diff(x,t)^2+diff(y,t)^2);
I = int(F,t,0,2*pi);
vpa(I)
得到结果为
ans = 
-25865.95635790842439629780459932
也可以直接使用 数值方法求解,把上述代码最后两行改成
quadl(inline(vectorize(F)),0,2*pi)
或
quadl(matlabFunction(F),0,2*pi)
得到结果为
ans =
 -2.5866e+004
 
第2题
属于第一类曲面积分,计算原理可以参考附件ppt的98-99页.
用MATLAB计算可以使用符号积分函数int,但无法求出解析解,然后使用vpa得到所需精度的数值结果.
代码如下：
syms x y
z = 2*(x-y^2)^3;
f = x^2+y^2+z^2;
F = f * sqrt(1+diff(z,x)^2+diff(z,y)^2);
I = int(int(F,x,-6,6),y,-6,6)
vpa(I)
得到结果：
ans = 
8829165547218478.3258650115855033
 
需要说明的是,使用int+vpa的方法耗时较长,而且不同MATLAB版本上的情况可能存在差别.我在2008b上计算两重int耗时超过半小时,得到最后结果共耗时秒（和具体硬件、软件环境相关,仅供参考）.
 
如果采用数值积分,由于在x=-6,y=±6附近函数的值较大（下图为把曲面积分转换成普通二重积分的被积函数曲面图形,使用 ezmesh(F,[-6 6]) 绘制）,如果用默认的误差限（1e-6）,计算时会出现大量的“Maximum function count exceeded; singularity likely”警告信息,计算精度无法得到保证.

使用大一些的误差限能避免该问题：
>> dblquad(inline(vectorize(F)),-6,6,-6,6,0.1,@quadl)
ans =
  8.8292e+015
或者,使用匿名函数的形式比起inline函数速度更快一些：
dblquad(matlabFunction(F),-6,6,-6,6,0.1,@quadl)
 
第3题（a）
属于常微分方程的初值问题,可用ode系列函数求解.
微分方程看不太清楚,等号右边第1项是不是y的二阶导数?
dy = inline('[y(2); y(3); 2*y(3)+6*x*y(1)]','x','y');
y0 = [1;1;1];
[T1, Y1] = ode45(dy,[0 10], y0);
[T2, Y2] = ode45(dy,[0 -10], y0);
T = [T2(end:-1:2); T1];
Y = [Y2(end:-1:2, :); Y1];
plot(T,Y)
得到结果看上去是发散的（不确定方程是否有问题）：

 
第3题（b）
属于边值问题（BVP）,可用bvp4c求
dy = inline('[y(2); y(3); 2*y(3)+6*x*y(1)]','x','y');
bc = inline('[ya(1)-2; yb(1); yb(2)]','ya','yb');
solinit = bvpinit(linspace(0,5,5),[0 0 0]);
sol = bvp4c(dy,bc,solinit);
x = linspace(0,5);
y = deval(sol,x);
subplot 211
plot(x,y(1,:));
ylabel('y')
subplot 212
plot(x,y(2,:));
xlabel('x')
ylabel('y''')

 

有关于Matlab的两道计算题1.某个函数（类似于z=2(x-y^2)^3）在x^2+(y/4)^2=1这个平面曲线上的积分.2.f=x^2+y^2+z^2在函数z=2(x-y^2)^3的空间曲面上的积分.其中x跟y的积分区间均为[-6,6]. 有没有关于某个电视的读后感 关于杠杆的两道计算题, 有没有关于密度的计算题? 南京有没有什么关于某个名胜或某个地方名称由来的民间传说 关于静电场中的导体的两道计算题 MATLAB中有表示某元素属于某个集合的函数吗?例如：表示,某数A属于集合B{...} 在matlab中,在一个矩阵中,如要想导出一个有某个数的行或列,求编程 MATLAB怎么求图像中某个灰度值的坐标? 有没有关于化学方程式的计算题不是化学式,越多越好! 有没有关于压强、密度的科学计算题 谁有关于电压锅的物理计算题 有没有关于一元二次方程的计算题和应用题, 谁有关于动能及功率的计算题 关于数学题的计算题 关于对数的计算题, 关于分式的计算题 关于密度的计算题