24小时热门版块排行榜    

查看: 1439  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

shuitong110

银虫 (小有名气)

[求助] 程序有误!!MTLAB 对函数积分(变量L),积分后的函数M再进行sin函数计算得到函数MM 已有2人参与

程序有误!!MTLAB 对函数积分(变量L),积分后的函数M再进行sin函数计算得到函数MM,最后绘制关于L,MM的图。 绘制的代码如下,能出图,但是和书上的结果差别很大,出来的图形应该是近似于正弦函数的图。我运算出来的怎么是一条直线,求大侠解答!!!!
clear;
clear;
L=0:0.02*10^-6:70*10^-6;

d=0.1*10^-6;
R=100*10^-6;
P=0.3;
k=2*pi./(1.55*10^-6);
B=11*10^6;
n1=3.45;
n2=1;
w=0.5*10^-6;
y1=sqrt(abs(k.^2*n1.^2-B.^2));
y2=sqrt(abs(B.^2-k.^2*n2.^2));

A=2*P*y1.^2*y2.^2;
C=B*k*k*(n1.^2-n2.^2)*(2+y2*w);
D=exp(-y2*d-y2*R+y2.*sqrt(R*R-L.*L));
M=A.*D./C;

S=trapz(L,M);
MM=sin(S);
plot(L,MM);

程序有误!!MTLAB 对函数积分(变量L),积分后的函数M再进行sin函数计算得到函数MM
1111.gif
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alorsmail

铁杆木虫 (著名写手)

引用回帖:
14楼: Originally posted by shuitong110 at 2014-05-09 08:18:05
??? Undefined function or method 'integral' for input arguments of type 'function_handle'.

Error in ==> xiaimuchong2 at 17
    M(i)=sin( integral(fun,0, v) ) ;

仿不出来。...

可能你的MATLAB版本比较旧,没有这个函数
我现在用的是2012B,
把其中相关的两个函数复制出来,(我没有旧版的MATLAB,无法测试是否足够)
你放到同一个目录下试试看,

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : integral.zip
  • 2014-05-09 12:44:46, 6.72 K
15楼2014-05-09 12:47:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

alorsmail

铁杆木虫 (著名写手)

最好能提供积分的原方程,这样方便大家对比
2楼2014-04-25 12:30:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuitong110

银虫 (小有名气)

引用回帖:
2楼: Originally posted by alorsmail at 2014-04-25 12:30:48
最好能提供积分的原方程,这样方便大家对比

已经将函数上传至图片格式 ,希望给与指点。谢谢!
3楼2014-04-25 14:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuitong110

银虫 (小有名气)


fegg7502: 金币+1, 鼓励交流 2014-05-09 08:10:51
图示函数
程序有误!!MTLAB 对函数积分(变量L),积分后的函数M再进行sin函数计算得到函数MM-1
函数.jpg

4楼2014-04-25 14:52:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见