张量功率谱CAMB参数调试

最新推荐文章于 2021-03-17 02:20:24 发布

Mr-Cat伍可猫

最新推荐文章于 2021-03-17 02:20:24 发布

阅读量689

点赞数

分类专栏： CAMB参数调动文章标签：调试 power camb

本文链接：https://blog.csdn.net/Mr_Cat123/article/details/78550058

版权

CAMB参数调动专栏收录该内容

6 篇文章 0 订阅

订阅专栏

（对应的原始tm文章在CMB\CosmoMC\camb\power spectra下，名为：primordial_power_spectrum_analysis.tm）
这里写图片描述

上图是在导入数据后再绘图的，matlab代码如下：

clc
figure
plot(L,sqrt(TT))
hold on
plot(L1,sqrt(TT1))
hold on 
plot(L2,sqrt(TT2))
hold on 
plot(L3,sqrt(TT3))
hold on 
plot(L4,sqrt(TT4))
axis([0,1000,-inf,inf])
legend('nt=0','nt=0.01','nt=0.05','nt=0.1','nt=0.5')
title('Tensor power spectrum CTT,nrun=0')

第三部分
此部分在谱指数nt=0，改变第一个跑动nrun的情况下查看角功率谱如何改变，由第一部分数学分析已经可以知道，此项跟标量功率谱第一个跑动完全类似，因此是谱是随着nrun增加而增加的。
以下分别为TT,EE,TE谱。
1)TT

clc
figure
loglog(L,(TT))
hold on
loglog(L1,(TT1))
hold on 
loglog(L2,(TT2))
hold on 
loglog(L3,(TT3))
hold on 
loglog(L4,(TT4))
axis([0,2200,-inf,inf])
legend('nr=0','nr=0.01','nr=0.05','nr=0.1','nr=0.5')
title('Tensor power spectrum CTT,nt=0,log coordinate')

这里写图片描述
2)EE

clc
figure
loglog(L,(EE))
hold on
loglog(L1,(EE1))
hold on 
loglog(L2,(EE2))
hold on 
loglog(L3,(EE3))
hold on 
loglog(L4,(EE4))
axis([0,1000,-inf,inf])
legend('nr=0','nr=0.01','nr=0.05','nr=0.1','nr=0.5')
title('Tensor power spectrum CEE,nt=0,log coordinate')

这里写图片描述
3）TE

clc
figure
semilogx(L,(TE))
hold on
semilogx(L1,(TE1))
hold on 
semilogx(L2,(TE2))
hold on 
semilogx(L3,(TE3))
hold on 
semilogx(L4,(TE4))
axis([0,2200,-40,40])
legend('nr=0','nr=0.01','nr=0.05','nr=0.1','nr=0.5')
title('Tensor power spectrum CTE,nt=0,log coordinate')