surf_integral.zip

  • segment
    了解作者
  • matlab
    开发工具
  • 2KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 10 积分
    下载积分
  • 1
    下载次数
  • 2020-08-10 13:32
    上传日期
MATLAB语言并未直接提供曲面积分的现成函数,因此,此处给出计算曲面积分的函数。
surf_integral.zip
  • surf_integral.m
    1.9KB
内容介绍
function I = surf_integral(f,vars,t,a,b) %surf_integral %第一类曲面积分 % I = surf_integral(f, z, [x,y], [x_m,x_M], [y_m,y_M]) % I = surf_integral(f, [x,y,z], [u,v], [u_m,u_M], [v_m,v_M]) % Examples: % 计算int_int(x^2*y+z*y^2)dS, 积分曲面如下: % x=ucosv, y=usinv, z=v, 0<=u<=a, 0<=v<=2*pi % MATLAB求解语句 % syms u v; syms a positive; % x=u*cos(v); y=u*sin(v); z=v; f=x^2*y+z*y^2; % I = surf_integral(f,[x,y,z],[u,v],[0,a],[0,2*pi]) % %第二类曲面积分 % I = surf_integral([P,Q,R], z, [x,y], [x_m,x_M], [y_m,y_M]) % I = surf_integral([P,Q,R], [x,y,z], [u,v], [u_m,u_M], [v_m,v_M]) % 注意:I = int_int_S(P*dydz+Q*dxdz+R*dxdy) % Examples: % 计算曲面积分int_int(x*y+z)dxdy, 积分曲面如下: % (x/a)^2+(y/b)^2+(z/c)^2=1的上半部,且积分沿椭球面的上面。 % 引入参数方程:x=a*sin(u)*cos(v),y=b*sin(u)*sin(v),z=c*cos(u); % 且0<=u<=pi/2, 0<=v<=2*pi 。 % MATLAB求解语句 % syms u v; syms a b c positive; % x=a*sin(u)*cos(v); y=b*sin(u)*sin(v); z=c*cos(u); % I = surf_integral([0, 0, x*y+z],[x,y,z],[u,v],[0,pi/2],[0,2*pi]) if length(f)==1 if length(vars)~=1 E = simplify(sum(diff(vars,t(1)).^2)); F = sum(diff(vars,t(1)).*diff(vars,t(2))); G = simplify(sum(diff(vars,t(2)).^2)); else E = simplify(1+diff(vars,t(1))^2); F = diff(vars,t(1))*diff(vars,t(2)); G = simplify(1+diff(vars,t(2))^2); end I = int(int(simplify(f*sqrt(E*G-F^2)),t(1),a(1),a(2)),t(2),b(1),b(2)); else if length(vars)~=1 A = diff(vars(2),t(1))*diff(vars(3),t(2)) - diff(vars(3),t(1))*diff(vars(2),t(2)); B = diff(vars(3),t(1))*diff(vars(1),t(2)) - diff(vars(1),t(1))*diff(vars(3),t(2)); C = diff(vars(1),t(1))*diff(vars(2),t(2)) - diff(vars(2),t(1))*diff(vars(1),t(2)); else A = - diff(vars,t(1)); B = - diff(vars,t(2)); C = 1; end f = f(:); abc = [A, B, C]; I = int(int(simplify(abc*f),t(1),a(1),a(2)),t(2),b(1),b(2)); end
评论
    相关推荐
    • MATLAB:MATLAB教程和MATLAB代码
      MATLAB:MATLAB教程和MATLAB代码
    • matlabcnhelp.rar
      matlab中文帮助很难找的,快速下载
    • MobilePolice.rar
      移动警察,车牌识别,车牌定位系统源代码,已经运用在移动车载稽查系统中。
    • SVM(matlab).rar
      支持向量机(SVM)实现的分类算法源码[matlab]
    • svm.zip
      用MATLAB编写的svm源程序,可以实现支持向量机,用于特征分类或提取
    • Classification-MatLab-Toolbox.rar
      模式识别matlab工具箱,包括SVM,ICA,PCA,NN等等模式识别算法,很有参考价值
    • VC++人脸定位实例.rar
      一个经典的人脸识别算法实例,提供人脸五官定位具体算法及两种实现流程.
    • QPSK_Simulink.rar
      QPSK的Matlab/Simulink的调制解调仿真系统,给出接收信号眼图及系统仿真误码率,包含载波恢复,匹配滤波,定时恢复等重要模块,帮助理解QPSK的系统
    • LPRBPDemo2009KV.rar
      车牌识别,神经网络算法,识别率高达95%,识别时间低于80ms。
    • MODULATION.RAR
      这个源程序代码包提供了通信系统中BPSK,QPSK,OQPSK,MSK,MSK2,GMSK,QAM,QAM16等调制解调方式 用matlab的实现,以及它们在AWGN和Rayleigh信道下的通信系统实现及误码率性能