三重积分计算--切片法
摘要：
以上半球为例，切片法将三重积分，转换成对z的定积分的以及对xy的二重积分 阅读全文
posted @ 2017-06-20 17:03 遗忘海岸 阅读(3825) 评论(0) 推荐(0) 编辑
2017年6月20日 #
2017年6月6日 #
matlab 向量场线积分
摘要：
syms t x y z F x=cos(t); y=sin(t); z=2*sin(t)^2-1; F=[x^2*y , (1/3)*x^3,x*y ] ; %场函数 V=[diff(x,t),diff(y,t),diff(z,t)]; r= int(dot(F,V),t,0,pi*2) %点积后 阅读全文
posted @ 2017-06-06 07:16 遗忘海岸 阅读(562) 评论(0) 推荐(0) 编辑
2017年6月3日 #
平面切圆柱面的椭圆绘制
摘要：
t=-sqrt(12):0.0001:sqrt(12); x=t; y=sqrt(12-t.^2); z=2*t + 3*y; plot3(x,y,z); hold on y1=-sqrt(12-t.^2); z=2*t + 3*y1; plot3(x,y1,z); grid on axis equ 阅读全文
posted @ 2017-06-03 07:23 遗忘海岸 阅读(938) 评论(0) 推荐(0) 编辑
2017年5月26日 #
抛物面倾斜体积积分
摘要：
假设抛物面里方满水，现在倾斜30度后，一部分水到出，那里余下的水平面正好过边缘点(0,1,1/2) 现在建立坐标，抛物面是y=1/2(x^2+y^2)，而水面是倾斜的。 将水到空的倾斜角度时，碗口处的一点的切线斜率对应的角度， 将碗逐渐倾斜，那么水面椭圆的y轴长度逐渐缩短，即割线在碗口变成切线时水倒 阅读全文
posted @ 2017-05-26 13:29 遗忘海岸 阅读(936) 评论(0) 推荐(0) 编辑
2017年5月25日 #
2017年5月24日 #
2017年5月5日 #
计算误差函数的积分--erf(x)
摘要：
注意exp(-x^2-y^2) 可以拆成exp(-x^2) * exp(-y^2) 对 dx积分时，另外一个可以当常数提出到积分号外，另外 -b到b的积分看做常数提取到积分号外 阅读全文
posted @ 2017-05-05 15:24 遗忘海岸 阅读(8817) 评论(0) 推荐(0) 编辑
2017年4月27日 #
三棱椎的体积
摘要：
另外可以以平行yz平面在x轴上切直角三角形类计算体积，结果一样都是 abc对应长方体体积的1/6 , 将一个长方体盒子按各对角线剪刀成三棱椎，刚好是4个，但是体积不是1/4 abc,因为如果拿实心的长方体剪裁中间还有一个四面体（四个面是全等三角形） 正方体如此剪裁后中间的是正四面体，应次可以有这个办 阅读全文
posted @ 2017-04-27 08:50 遗忘海岸 阅读(379) 评论(0) 推荐(0) 编辑
2017年4月26日 #
2017年4月25日 #
matlab 设定坐标比例
摘要：
figure(1) u=-0.1:0.005:0.1; v=-0.1:0.005:0.1; [x,y]=meshgrid(u,v); z=sin(x-y)./abs(x)+abs(y); surf(x,y,z); hold on % contourf(x,y,z) % ezplot('x+y=3') 阅读全文
posted @ 2017-04-25 08:53 遗忘海岸 阅读(772) 评论(0) 推荐(0) 编辑
2017年4月19日 #
2017年4月18日 #
Lagrange 乘子法求最优解
摘要：
clc clear syms x y z r1 r2 w f=x^2+y^2+z^2+w^2; g1=2*x-y+z-w-1; g2=x+y-z+w-1; h=f-r1*g1 -r2*g2; hx=diff(h,x); hy=diff(h,y); hz=diff(h,z); hw=diff(h,w) 阅读全文
posted @ 2017-04-18 16:24 遗忘海岸 阅读(252) 评论(0) 推荐(0) 编辑
2017年4月14日 #
二元函数临界点的局部极值
摘要：
close figure(2) syms x y f=2.*x.^4+y.^4-2.*x.^2-2.*y.^2+3; fx=diff(f,x); fy=diff(f,y); fxx=diff(f,x,2); fyy=diff(fy,y); fxy=diff(fx,y); fx fy fxx fyy 阅读全文
posted @ 2017-04-14 13:38 遗忘海岸 阅读(418) 评论(0) 推荐(0) 编辑
2017年4月13日 #