MATLAB计算引擎应用.ppt

下载文档 降价啦

29
0
约1.28万字
约 90页
2018-02-18 发布于浙江
举报
版权申诉
保障服务

MATLAB计算引擎应用.ppt

1、本文档共90页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

MATLAB计算引擎应用

012 C 013 integer ep, T, D, V 014 real*8 time(11), dist(11), finitV 015 integer engPutVariable, engEvalString, engClose 016 integer temp, status 017 data time / 0.0,0.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0 / 018 write(6,*) 开始MATLAB计算引擎应用程序...... 019 C 020 C 打开计算引擎 021 C 022 ep = engOpen(\0) 023 　　　　　 if (ep .eq. 0) then 024 write(6,*) 无法打开计算引擎! 025 stop 026 　 endif 027 C 创建数据 028 T = mxCreateDoubleMatrix(1, 11, 0) 029 call mxCopyReal8ToPtr(time, mxGetPr(T), 11) 030 finitV = 20.0 031 V = mxCreateScalarDouble(finitV) 032 C 033 C 将变量T,V放置到MATLAB工作空间 034 C 035 status = engPutVariable(ep, T, T) 036 if (status .ne. 0) then 037 write(6,*) 无法设置变量 --- T! 038 stop 039 endif 040 status = engPutVariable(ep, V, V) 041 if (status .ne. 0) then 042 write(6,*) 无法设置变量 --- V! 043 stop 044 endif 045 C 046 C 执行MATLAB指令 047 C 048 if (engEvalString(ep, D = V.*T + 0.5.*(-9.8).*T.^2;) .ne. 0) 049 * then 050 　　write(6,*) 无法正确执行指令! 051 　　　stop 052 endif 053 C 绘制结果 054 if (engEvalString(ep, plot(T,D);) .ne. 0) then 055 write(6,*) 无法执行指令! 056 stop 057 endif 058 if (engEvalString(ep, title(位置 vs. 时间)) .ne. 0) then 059 write(6,*) 无法执行指令! 060 stop 061 endif 062 if (engEvalString(ep, xlabel(时间 (秒))) .ne. 0) then 063 write(6,*) 无法执行指令! 064 stop 065 endif 066 if (engEvalString(ep, ylabel(位置 (米))) .ne. 0)then 067 　　　write(6,*) ‘无法执行指令!’ 068 　　　stop 069 endif 070 write(6,*) 按回车键将继续...... 071 pause 072 C 获取计算结果 073 D = engGetVariable(ep, D) 074 call mxCopyPtrToReal8(mxGetPr(D), dist, 11) 075 write(6,*) MATLAB 的计算结果如下: 076 write(6,*) time(s) distance(m) 077 do 10 i=1,11 078 print 20, time(i), dist(i) 079 20 format( , G10.3, G10.3) 080 10 continue 081 write(6,*) 按回车键将关闭计算引擎...... 082 pause 083 C 关闭图形窗体 084 if (engEvalString(ep, close;) .ne. 0) then 085