查看: 98|回复: 3|关注: 0

[已答复] 关于级数及其函数图像

[复制链接]

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 4 天前 | 显示全部楼层 |阅读模式
L=1:1:300;
mu=0.1;etaBob=0.045;alpha=0.2;Pe=3.2*10.^(-7);pdark=1.7*10.^(-6);
eta=(mu*10.^((-alpha.*L)/10).*etaBob)./((mu*10.^((-alpha.*L)/10).*etaBob)+2*Pe);edetector=0.033;
Y0=pdark;
e=exp(1);
Qu=Y0.*e.^(-mu);
for n=1:1:100
    etan=1-(1-eta).^n;
    Yn=etan+pdark;
    en=(edetector.*etan+((1./2).*pdark))./Yn
    Qu=Qu+Yn*e.^(-mu)*((mu.^n)./factorial(n));
    Eu=(e.^(-mu)*((mu.^n)./factorial(n))*en)./Qu;
end
Q1=(pdark+eta)*0.045*(mu./(1+mu).^2);
e1=((0.5*pdark)+(etaBob*eta))./(pdark+eta);
a1=log2(Eu);
a2=log2(1-Eu);
b1=log2(e1);
b2=log2(1-e1);
R=-Qu.*1.15.*(-Eu.*a1-((1-Eu).*a2))+Q1.*(1-(-e1.*b1-((1-e1).*b2)));
plot(L,R)
第一种:当n取值到100的时候能够画出图像,但将n取值到200就无法画出图像,只有坐标。
第二种:用sym n就无法画出图像,显示“错误使用 plot
数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。
出错 Untitled2 (line 12)
plot(x,R)”
回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

7

主题

1479

帖子

316

最佳答案
  • 关注者: 75
发表于 4 天前 | 显示全部楼层
仅供参考
  1. L=1:1:300;
  2. mu=0.1;etaBob=0.045;alpha=0.2;Pe=3.2*10.^(-7);pdark=1.7*10.^(-6);
  3. eta=(mu*10.^((-alpha.*L)/10).*etaBob)./((mu*10.^((-alpha.*L)/10).*etaBob)+2*Pe);edetector=0.033;
  4. Y0=pdark;
  5. e=exp(1);
  6. Qu=Y0.*e.^(-mu);
  7. for n=1:1:200
  8.     etan=1-(1-eta).^n;
  9.     Yn=etan+pdark;
  10.     en=(edetector.*etan+((1./2).*pdark))./Yn;
  11.     Qu=Qu+Yn*e.^(-mu)*((mu.^n)./factorial(sym(n)));
  12.     Eu=(e.^(-mu)*((mu.^n)./factorial(sym(n)))*en)./Qu;
  13. end
  14. Q1=(pdark+eta)*0.045*(mu./(1+mu).^2);
  15. e1=((0.5*pdark)+(etaBob*eta))./(pdark+eta);
  16. a1=log2(Eu);
  17. a2=log2(1-Eu);
  18. b1=log2(e1);
  19. b2=log2(1-e1);
  20. R=-Qu.*1.15.*(-Eu.*a1-((1-Eu).*a2))+Q1.*(1-(-e1.*b1-((1-e1).*b2)));
  21. plot(L,R)
复制代码
回复此楼 已获打赏: 0 积分

举报

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
 楼主| 发表于 4 天前 | 显示全部楼层

感谢!请问如果n取到无穷可以做吗,我考虑到n取无穷做不出来,所以将n的取值固定
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

7

主题

1479

帖子

316

最佳答案
  • 关注者: 75
发表于 4 天前 | 显示全部楼层
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表