目前在做一个HP滤波分析，教材上对HP滤波中lambda的选择，月度数据通常用14400.但是看一些国外的文献中采用的是129600.想请问如何在月度的时间序列中选择最优的lambda?求程序。谢谢！
clear
data=xlsread('D:\china.xlsx');

lambda = 100:10:100000;
N = length(lambda);
TS = zeros(length(data),N);
CS = TS;
for i = 1:N
    [TS(:,i),CS(:,i)] = hpfilter(data(:,2),lambda(i));
end
lambda_s = (std(diff(CS))./std(diff(TS))).^2;
[lambda_opt,maxi] = max(lambda./lambda_s);
[T,C] = hpfilter(data,lambda(maxi));



这样做下来，lambda_opt=87.0130，问题在哪里？求解答。谢谢！