0

収集したアクセラレータデータを、記録されたLVDT変位と比較して検証しようとしています。これを行うために、加速を得るためにLVDTデータを2回区別しようとしています。ただし、コードを実行するとエラーが発生します

Undefined function or variable 'DiffDiffLVDT'.

これを調査すると、Matlabが2番目のforループを処理していないため、変数DiffDiffLVDTが作成されないことがわかりました。

2番目のforループをスキップするのはなぜですか？

enter code here
clear all

clc

test2 = csvread('HZ1.csv'); %raw data from excel
time = test2(:,1);
%% Filtering 
f=250;%  sampling frequency
f_cutoff = 2; % cutoff frequency
f_cutoff2 = 5;

fnorm =f_cutoff/(f/2); % normalized cut off freq, you can change it to any value depending on your requirements
fnorm2 =f_cutoff2/(f/2);

[b1,a1] = butter(4,fnorm,'low'); % Low pass Butterworth filter of order 4
[b2,a2] = butter(4,fnorm2,'low'); % Low pass Butterworth filter of order 4

filtaccZ = filtfilt(b2,a2,test2(:,6));
filtLVDTZ = filtfilt(b1,a1,test2(:,7));

%% Remove Offset

Accz = filtaccZ -mean(filtaccZ);
LVDTz = filtLVDTZ - mean(filtLVDTZ);

%% Unit Conversion

LVDTm = LVDTz/1000; % mm to m

%% LVDT Displacement to Acc

% Displacement to Velocity
for d = 2:1:(size(LVDTm)-1)

    x = LVDTm(d+1);

    y = LVDTm(d-1);

    z = x-y; %differnce in y

    a = time (d+1);

    b = time (d-1);

    c = a-b; %differnce in x

    DiffLVDT(d)= (z/c); % Displacement to velocity

end

velocity = DiffLVDT;

% Velocity to Acceleration
for e=1:1:(size(velocity)-1)
    x2 = velocity(e+1); 

    y2 = velocity(e-1);

    z2 = x2-y2; %differnce in y

    a2 = time (e+1);

    b2 = time (e-1);

    c2 = a2-b2; %differnce in x

    DiffDiffLVDT(e)= (z2/c2) %velocity to acc.

end

Acc= DiffDiffLVDT
%% Plotting
close all 
figure
hold on
plot(time(1:5000),Acc(1:5000),'b')
plot(time(1:5000),Accz(1:5000),'r')

grid on;box on
legend('DiffDiffLVDTFilter','Accz')

    enter code here
4

1 に答える 1

4

なぜなら

1:1:(size(velocity)-1)

あなたが望むことをしません。

velocityは1xN配列であるsize(velocity)ため、を返します[1 N]。コロン演算子は配列の最初の値のみを考慮し、ループする配列は空になります。

1:1:[1 N]-1 == 1:1:1-1 == 1:1:0 == Empty matrix

このような問題があるため、size(var)ループを作成するときは、常に次のいずれかを使用する必要があります。

size(var, n) % Returns the size of the nth dimension of var
length(var) % Returns the size of the largest dimension of var
numel(var) % Returns the number of elements in var
于 2013-03-26T22:32:47.057 に答える