CMB标量功率谱第二个谱指数跑动项n(2)跑动带来的影响

由前一篇文章的分析知道，在n(1)=const或者取零（此处取0）的情况下，增加n(2)会出现谱前面下降后面上升的情况，因此本文从图形上来分析是否出现这种情况。
导入数据，然后执行以下代码：

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('nrr=0','nrr=0.01','nrr=0.05','nrr=0.1','nrr=0.5')
title('scalar power spectrum CTT,nrun=0,log coordinate')

不难发现TT谱满足所描述的那样，在前半部分nrunrun越大，谱越往下移动，后半部分nrunrun越大，谱越小；但总体来看后期变化幅度是很小的，说明小尺度上近似不变的。在谱上升和下降的交叉口即是k=ks=0.05的时候
同样给出EE谱和TE谱如下：
1）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,2200,-inf,inf])
legend('nrr=0','nrr=0.01','nrr=0.05','nrr=0.1','nrr=0.5')
title('scalar power spectrum CEE,nrun=0,log coordinate')

2)TE谱

%{
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,-inf,inf])
legend('nrr=0','nrr=0.01','nrr=0.05','nrr=0.1','nrr=0.5')
title('scalar power spectrum CTE,nrun=0,log coordinate')