2011 级物理专业计算物理上机考试题.docVIP

下载本文档

3
0
约 7页
2017-08-13 发布于安徽
举报
版权申诉

2011 级物理专业计算物理上机考试题.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

2011 级物理专业计算物理上机考试题 th=RandomInteger[{1,12}] Date[ ] nm=ToString[Input[ 请输入姓名 :]] Export[nm .dat,{th,Date[]}] yy=Import[nm .dat] 解初值问题 1 eq={y[t]+0.01y[t]+25?^2y[t]==Sin[? t]Exp[-t/5],y[0]==10,y[0]==-2}; 答案： y (t)=0.0202819 E^(-0.1 t)+9.97972 E^(-1.125 t) Cos[15.6676 t]+0.589063 E^(-1.125 t) Sin[15.6676 t] 2 数据拟合问题自变量采样为 ts={0.,0.02,0.04,0.06,0.08,0.1,0.12,0.14,0.16,0.18,0.2,0.22,0.24,0.26,0.28,0.3,0.32,0.34,0.36,0.38,0.4,0.42,0.44,0.46,0.48,0.5,0.52,0.54,0.56,0.58,0.6,0.62,0.64,0.66,0.68,0.7,0.72,0.74,0.76,0.78,0.8,0.82,0.84,0.86,0.88,0.9,0.92,0.94,0.96,0.98,1.,1.02,1.04,1.06,1.08,1.1,1.12,1.14,1.16,1.18,1.2,1.22,1.24,1.26,1.28,1.3,1.32,1.34,1.36,1.38,1.4,1.42,1.44,1.46,1.48,1.5,1.52,1.54,1.56,1.58,1.6,1.62,1.64,1.66,1.68,1.7,1.72,1.74,1.76,1.78,1.8,1.82,1.84,1.86,1.88,1.9,1.92,1.94,1.96,1.98,2.}; 函数值采样为 ys={10.,9.48032,8.07801,5.96613,3.37805,0.582433,-2.14399,-4.54332,-6.40028,-7.56101,-7.94503,-7.54978,-6.44739,-4.77448,-2.71629,-0.486786,1.69309,3.61687,5.11167,6.05329,6.37595,6.07598,5.2095,3.88436,2.24758,0.469541,-1.27336,-2.81584,-4.01909,-4.78286,-5.05339,-4.8265,-4.14578,-3.09638,-1.79496,-0.377188,1.01611,2.25265,3.22096,3.84016,4.06633,3.89513,3.36042,2.52937,1.49456,0.364007,-0.749866,-1.74118,-2.52045,-3.02242,-3.21124,-3.08263,-2.66286,-2.00494,-1.18231,-0.280978,0.60933,1.40387,2.03081,2.4375,2.59471,2.49837,2.16884,1.64793,0.993922,0.275259,-0.436425,-1.07332,-1.57776,-1.90729,-2.03807,-1.96633,-1.70785,-1.2956,-0.775799,-0.202941,0.365805,0.876176,1.28189,1.54871,1.65721,1.60396,1.40118,1.07485,0.661641,0.204925,-0.249674,-0.658737,-0.985118,-1.20122,-1.29121,-1.25203,-1.09312,-0.834951,-0.506616,-0.14263,0.220595,0.548325,0.810746,0.985613,1.06004}; data=Transpose[{ts,ys}]; 建议目标函数选为yg[t_]=a1 ?^(-a2 t)+a3 ?^(-a4 t)Cos[a5 t+a6]; 请做非线性数据拟合; 给出参数 {a1, a2, a3, a4, a5, a6} 的值; 给出目标函数; 画出函数曲线; 答案为; ; ; {a1?0.0202822,a2?0.100003,a3?9.99709,a4?1.125,a5?15.6676,a6?-0.0589578} 3 为求边值问题3 求边值问题的数值解