实验程序 直立台阶.doc_第1页
实验程序 直立台阶.doc_第2页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

%Gravityball.m%功能：直立台阶模型的重力异常%入口参数：x,y-平面坐标 rd-球体剩余密度% h地面到台阶上部距离 H-地面到台阶底部距离%出口参数：g-重力异常 Vxz, Vzz-二阶导数 Vzzz-三阶导数 (剖面）% clear; %清除内存变量x=-15:1:15; %y=x; 相当于测线上的测点y=-10:1:10; %相当于测线X,Y=meshgrid(x,y); %-矩阵的方法进行计算(方法一的语句)rd=0.8; %球体剩余密度 Remain Dentisity (rd: g/cm3)h=3; %地面到台阶上部距离 (h:m)H=4; %地面到台阶底部距离 (H:m)G=6.67; %万有引力常数 (G:m3/(kg*s2) %用矩阵的方法计算（方法一）%g1=G*rd*(pi*(H-h)+2*H*atan(X/H)-2*h*atan(X/h)+X.*log(X.2+H.2)./(X.2+h.2); %计算重力异常,用矩阵的方法进行计算(方法一的语句)g=g1*10-2; %(g: g.u.)(方法一的语句)%-以下均用方法一来计算其他分量,也可改写成循环(方法二)来计算-V1xz=2*G*rd* log(X.2+H.2)./(X.2+h.2); %计算重力位X方向二阶导数Vxz=V1xz*10; %(Vxz: E )1E=10-9 1/s2V1zz=2*G*rd*( atan(X/h)- atan(X/H); %计算重力位Z方向二阶导数Vzz=V1zz*10;V1zzz=2*G*rd*X.*(1./( X.2+h.2)-1./(X.2+H.2); %计算重力位三阶导数Vzzz=V1zzz*104; %(Vzzz:pMKS) 1pMKS=10-12 1/m*s2%-画等值线图-figure(1);%C,H=contour(g,20); colormap(hsv) ;clabel(C,H); title(Gravity contour) %等值线有标注C,H=contour(g,20); title(Gravity contour) %等值线没标注%-画立体图-figure(2);subplot(2,2,1);surf(Vxz);title(Vxz);. subplot(2,2,2);surf(Vzz);title(Vzz);. subplot(2,2,3);surf(Vzzz);title(Vzzz);.%-画剖面图-figure(3);subplot(2,1,1);plot(x, g(1,:)*20,b*,x, g(11,:),r-);legend(g(1),g(11);title(*)subplot(2,1,2);plot(x,Vxz(
人人文库> 全部分类> 生活休闲 > 科普知识

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论