这是后续我之前遇到的一个问题 https://stackoverflow.com/questions/45688607.
我想给信号一个偏移量,然后在其中添加一些延迟并计算 RMSE,但是在求差时我遇到以下问题:
我想请教以下问题:
- 我该如何解决上述问题?
- 有人可以用简单的话解释一下什么吗
iddata
确实如此 - 因为我研究了包括 MATLAB 在内的不同门户,但仍然无法获得一个好的概念。
- 如何存储类型的数据
iddata
在我的代码的最后部分进行减法的单元格中?
有问题的代码:
drv(1)=load('123.mat');
t = drv(1).x;
ref = drv(1).y;
angle = drv(1).z;
Fs = 1000;
t1 =t';
ref1= ref';
d_data = iddata(ref1, t1, 1/Fs);
%% Add offset:
x = 1;
afterOffset1= {};
for i = 100:10:130
T = getTrend(d_data);
% <detrend data if needed>
T.InputOffset = i;
T.OutputOffset = i;
afterOffset = retrend(d_data,T);
afterOffset1{x,1}= afterOffset;
x= x+1 ;
end
%% Add delay:
y=20;
afterDelay1= {};
for i = 1:1:4
% delaySamples = i; % Must be a non-negative value
% afterDelay = iddata([NaN(delaySamples,1); d_data.OutputData],...
% [d_data.InputData; NaN(delaySamples,1)], 1/Fs);
afterOffset1{i}.Tstart = y;
afterDelay1{i,1}= afterOffset1{i};
y= y+10;
end
%% Plot:
n = size(afterDelay1,1);
figure();
for i=1:1:n
subplot(2,2,i);
plot(d_data);
hold all
plot(afterDelay1{i});
end
sig_diff = angle(1)-afterDelay1;
square_error(i,:) = (sig_diff(i)).^2;
mse(i,:)= mean(square_error(i));
rmse(i,:) = sqrt(mse(i));
sig_diff = d_data_1 - afterDelay; % <<<<<<<<<<<<<<<<<<<<<< Problem is here
% square_error = (sig_diff).^2;
% mse= mean(square_error);
% rmse = sqrt(mse);
end
您很可能想要OutputData
属性来自于iddata
对象是输出或y
您的问题的信号:
sig_diff = angle(1)-afterDelay1.OutputData;
另请注意,这将为您提供column向量,但稍后您的代码会假设它是行向量。执行上述计算后,您可能需要转置此数据,然后再继续:
sig_diff = angle(1)-afterDelay1.OutputData;
sig_diff = sig_diff.';
一般来说,iddata
是创建表示输入和输出时域或频域数据的对象的函数。请注意,当您创建iddata
对象,输入矩阵可能有多个源,因此每个column规定了一个来源。对于输出也可以这样说,其中每个column规定了输出。因此,在使用此函数之前转置数据非常重要,以确保每个信号位于单独的列中,或者仅使用单个列来表示一个输入/输出。
对象内部具有多种属性,包括采样时间或采样频率、函数所采用的有效域和范围以及最终访问输入和输出数据。OutputData
是这些字段之一。我建议您查看有关您可以访问的所有属性的文档iddata
. OutputData
这里有明确的定义:https://www.mathworks.com/help/ident/ref/iddata.html https://www.mathworks.com/help/ident/ref/iddata.html
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)