Построение графика с логарифмической осью - 3 значения y для каждого значения x

Question

Я пытаюсь построить график значения значения переменной x по оси x; Однако, есть три значения для каждой точки, то есть, собственная частота 1, частота собственных колебаний 2 и собственная частота 3. Это код, и он не работает:

M = [3 0 0; 0 2 0; 0 0 0.5] 
% disp('Mass matrix [M]'); 
% disp(M); 

for i=1:1000:60e06 
    K = [i+7e06 i -6e06; i i+3e06 -3e06; -6e06 -3e06 10e06]; 
    % disp(K) 
    [V,L]=eig(K,M);  % eig is a standard Matlab function 
    values=diag(L);  % diag gets the values from the leading diagonal 
    [values,I]=sort(values); % sort into ascending order of frequency 
    V(:,I)=V;     % now sort the mode shape vectors 
    freq=sqrt(values)/(2*pi); % this divides all elements by 2 pi 

    % disp('Eigenvalues [s^(-2)]'); 
    % disp(values');  % the quote mark after values prints the column vector as a row 

    % disp('Frequencies [Hz]'); 
    % disp(freq'); 

    % disp('Mode shape vectors in the columns of matrix [V]'); 
    % disp(V); 
     loglog(i, freq(1:1), i, freq(2:1), i, freq(3:1) 
end 

Извиняюсь за ошибки, Я начинающий.

Answer 1

Вы пропускаете закрывающую скобку в конце вашей линии

loglog(i, freq(1:1), i, freq(2:1), i, freq(3:1) 

Кроме того, я не вижу причин для вашего 2:1, 3:1, ...?!

Заменить его

loglog(i, freq(1), i, freq(2), i, freq(3)) 

, и я не вижу ничего плохого с вашим черчения.

В целом я предлагаю хранить значения freq в некотором массиве и рассчитать его отдельно. Это ускорит сценарий до много, т.е.

%preallocate freq, the number of iterations in the loop is important here 
freq = zeros(3,length(1:1e3:60e6)) 

for i = ... 
    %your loop goes here 
    %remove freq = ... 
    %remove loglog(... 
    freq(:,i) = sqrt(values)/(2*pi); % this divides all elements by 2 pi 
end 

loglog(1:1e3:60e6, freq(1,:)) %plot first curve 
loglog(1:1e3:60e6, freq(2,:)) %plot second curve 
loglog(1:1e3:60e6, freq(3,:)) %plot third curve 

Answer 2

Ваша проблема заключается в том, что у вас есть loglog внутри цикла.

Перед цикла, я поставил:

i_vals = 1:1000:60e06; 
freq = zeros(3, length(i_vals)); 

изменить верхнюю часть петли на:

for n=1:length(i_vals) 
    i = i_vals(n); 
    K = [i+7e06 i -6e06; i i+3e06 -3e06; -6e06 -3e06 10e06]; 
    ... 

Также изменить freq назначение:

freq(:, n)=sqrt(values)/(2*pi); 

Удаление loglog и поместите это после/вне вашей петли.

loglog(i_vals, freq') 

Делая все, что дает мне это: