From linear to logarithmic distribution

Question

Mohamed Abdalmoaty 2015-6-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/223174-from-linear-to-logarithmic-distribution

评论： Walter Roberson 2015-6-10

Given a linear frequency scale denoted by the vector f = 1:fmax; fmax is known, say 1000.

How can I construct, out of this vector, a logarithmic selection of the entries such that the ratio between two consecutive odd numbers is 1.1?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2015-6-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/223174-from-linear-to-logarithmic-distribution#answer_182172

在 MATLAB Online 中打开

You cannot. The only solution to x*11/10 = x+2 is x = 20, and that is an even number.

If the task is to construct a vector of length fmax such that the ratio

x(2*k+3)/x(2*k+1)

is 11/10 for k in nonnegative integers, then as logarithmic progression is x(n+1) = C * x(n) for some constant C, it follows that x(n+2) = C*(C*x(n)) = C^2 * x(n). You want that value, C^2 to be 1.1, so your C is sqrt(1.1). Then you get

initial_value * sqrt(1.1).^(0:fmax-1)

2 个评论
显示无隐藏无

Mohamed Abdalmoaty 2015-6-9

编辑：Walter Roberson 2015-6-10

在 MATLAB Online 中打开

fmax = 1000; Ratio = 1.1;
Lines_linear = 1:fmax/2;
currentIdx = 1; selected = 1;
Lines = 1;
while currentIdx < fmax/2
 [~,nextIdx]=min(abs(Lines_linear(currentIdx)*Ratio-...
        Lines(currentIdx+1:length(Lines_linear))));
  currentIdx = nextIdx+currentIdx;
  % a logarithmic step!
  if Lines_linear(currentIdx) > sqrt(Ratio)*Lines(length(Lines)),
    selected = [selected;currentIdx];
    Lines = [ Lines; Lines_linear(currentIdx)];
  end
end

Walter Roberson 2015-6-10

I have no idea what that code is intended to do. It is uncommented and the relationship to the Question you posted is not obvious.

You start with Lines=1 which is length 1. In the first statement in your while you try to index Lines from 2 to length(Lines_linear) which would be 2:500. But Lines only has 1 entry, so you cannot access elements 2 to 500 of it.

请先登录，再进行评论。