12第12章-其他建模方法_第1页
12第12章-其他建模方法_第2页
12第12章-其他建模方法_第3页
12第12章-其他建模方法_第4页
12第12章-其他建模方法_第5页
已阅读5页,还剩25页未读 继续免费阅读

下载本文档

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

文档简介

(12.32),求得的本地直角坐标和直接利用MATLAB工具箱求得的本地直角坐标差异很大,当然由于算法不同,有些差异是正常的,但这里的误差实在是太大了,我们这里就不给出计算结果了。计算的MATLAB程序如下:clc,cleara=6378137;e=0.081819190842622;B=[40.0614,40.0613,40.0487,39.9310,39.8429,39.6983,39.6140,39.6482]';%纬度L=[78.6817,78.6814,78.6396,77.9549,76.6504,76.1164,75.9389,75.7390]';%经度H=[1135,1139,1153,1167,1228,1342,1513,1737]';%高程N=a./sqrt(1-e^2*(sind(B)).^2);X=(N+H).*cosd(B).*cosd(L);%计算地心直角坐标系的X坐标Y=(N+H).*cosd(B).*sind(L);%计算地心直角坐标系的Y坐标Z=(N*(1-e^2)+H).*sind(B);%计算地心直角坐标系的Z坐标X0=X(1);Y0=Y(1);Z0=Z(1);%本地直角坐标系的坐标原点xyz=[-sind(B(1))*cosd(L(1)),-sind(B(1))*sind(L(1)),cosd(B(1))-sind(L(1)),cosd(L(1)),0-cosd(B(1))*cosd(L(1)),-cosd(B(1))*sind(L(1)),-sind(B(1))]*[X-X0,Y-Y0,Z-Z0]';xyz2=lla2flat([BLH],[B(1)L(1)],0,-H(1));%利用工具箱计算,输入顺序:纬度,经度,高程check=[xyz;xyz2']%两种计算结果比较习题1212.1利用MonteCarlo方法,模拟掷骰子各面出现的概率。12.2利用MonteCarlo方法,求积分,并与数值解的结果进行比较。12.3使用MonteCarlo方法,求椭球面所围立体的体积。12.4某家庭从现在着手从每月工资中拿出一部分资金存入银行,用于投资子女的教育,并计划20年后开始从投资账户中每月支出1000元,直到10年后子女大学毕业用完全部资金。要实现这个投资目标,20年内共要筹措多少资金?每月要向银行存入多少钱?假定投资的月利率为0.4%。12.5某大型企业1999年至2004年的产品销售额如REF_Ref630235\h表12.9所列,试建立GM(1,1)预测模型,并预测2005年的产品销售额。表12.SEQ表12.\*ARABIC9产品

温馨提示

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

评论

0/150

提交评论