How i do my RMS function on Embedded Matlab Function

2012 8 28
1 个回答
回答已采纳
2 次查看（30 天）

0 个投票

Hello.
I need to create my RMS block to use with simulink. This is the first part of a protection relay project. So i'll test many algoritms to do this. I already did this on m-file project but im want to do this a dinamic mode.
I've read about of coder.extrinsic.
Example:
% Mann & Morrison algoritm
f0=60;
w=2*pi*f0;
b=(1/w)^2;
c=1/sqrt(2);
fs=32*f0; % sample frenq. 32 sample per cycle.
dt=1/fs;
der=(I(n) - I(n-2))/(2*dt); %I is a imput (sine wave block)
Ief(n-1)=c*(sqrt(I(n-1)^2 + b*(der)^2)); %Ief is a RMS value at sample n
n=n+1;
I dont know how i work with arrays on the function block. I need a buffer?
Thanks

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

 采纳的回答

Ryan G
Ryan G 2012-8-29

0 个投票

If you don't need to generate code from this Embedded MATLAB block and you are not overly concerned with simulation speed then you can call the entire function extrinsically as long as the output of the block is always the same size. You will need to initialize the output regardless by doing something like:
y = 0;
before the function call.
If the output is variable size you will need to buffer the output. Say you know the output will never exceed 100 elements, you can initialize the output like so:
y = zeros(1,100);
and fill in the block as needed with the actual values as they become available.

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Arthur
Arthur 2012-8-29
Thanks a lot Ryan. I've working hard to learn about the S-Functions.
I used a buffer of 32 samples, the same way that you said. Take a look:
function [sys,x0] = RMS_Square(t,x,u,flag)
n=32;
if flag == 0,
sys = [0 32 1 1 0 0]; % length 32
x0 = zeros(1,32);
elseif abs(flag) == 2,
sys = zeros(32,1);
sys = x;
sys = circshift(sys,[1,0]); % rotate de array to get the new sample
sys(1) = u;
elseif abs(flag) == 3 && max(t)>0.016, %0,016 is a cycle t/60hz (wait until the buffer is complete)
sys=sqrt(sum(x.^2)/n); %rms
else
sys = [];
end
Ryan G
Ryan G 2012-8-30
I'm assuming this is still giving you issues I don't see x0 anywhere but the first branch of the if-statement. x0 will need to be initialized in all branches of the if statement because otherwise, say flag never equaled 0, x0 would never be initialized.
You may have an issue with circshift if it's not supported by code generation, however, you can set it with
coder.extrinsic('circshift')
if needed. Same may go for abs, but I bet that is supported.
Mike Hosea
Mike Hosea 2012-9-2
Support for circshift was added in R2008b.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 General Applications 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by