24小时热门版块排行榜    

查看: 1166  |  回复: 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的回帖

shuitong110

银虫 (小有名气)

引用回帖:
6楼: Originally posted by alorsmail at 2014-04-25 16:12:52
如果要得到一个曲线的话,你需要加一个循环,
示例如下,结果见图片
==========================
L=0:0.02*10^-6:140*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 ;
...

谢谢您的回复,
MM = zeros(size(L)) ;
for i = 10:length(L)
    Lp = L(1:i) ;
这三行是什么意思?
9楼2014-04-25 19:41:16
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见