RMS e MÉDIO EFICAZ.

[fabim]

DIZ A LENDA DE OBTENÇÃO DE VALOR RMS, CUJA ESTE É O VALOR EFICAZ DE i OU v EM t.

Sendo que a lenda diz que o valor RMS =

SQRT{ 1/T * (sum(N1²--Nn²)) }.

Sendo o meu T = 1.
E meu Nn=5860;

SQRT{ 1/1 * (sum(N1²--N5860²)) }/N5860.

Olhando ainda como exemplo em.

T0segundo T1/2segundo T1segundo.
20V 15V 10V

V RMS para N0 a N3.

SQRT{ 1/1 * (sum(N1²--N3²)) }/N3.
SQRT{ 725 }/N3.
26,925824035672520156253552457702 / N3.
VRMS=8,9752746785575067187511841525667

ONDE ESTOU ENGOLINDO CAROÇO ?

Porque.
Se ficar na média. 20+15+10/3 = médio 15 para 1 segundo.

Agora.
Vejam bem um retificador half bridge enxergando a saida de um TC.
levando em consideração que V- = V+ em forma e duração.
Uma amostragem de 5860 por segundo.
Então no final da acumulação de amostras.
ACC *= 2;

V média/eficaz = ACC / 5860;

CORRETO ?

Consegui me expressar ?

Porque revejam o calculo lá encima.
Suponha que uma bateria em tensão constante de 10V por 1 segundo e foram feitas 10 amostras.

SQRT{ 1/1 * (sum(N1²--N10²)) }/N10.
SQRT{ 1000 }/N10.
RESU = 3,1622776601683793319988935444327

como pode isso ?

abraços

Fabio

[Djalma Toledo Rodrigues]

.........._________________________
----- \/ (x1² + x2² .... xn-1² + xn²) / n '

No seu exemplo: SQRT{ 725 /N3} e não SQRT{ 725 }/N3.

SQRT{ 725 /N3} = ~ 15,546

Abraxi Fabim e boa RMS

[fabim]

aff... POIS É DJ.
tava fazendo mercadoria.

eita diaxo.

Agora deu certo..

E outra coisa.

Se eu sair somando tudo e dividir pelo numero de amostras em um segundo, o valor é praticamente o mesmo só quebra o 0,XX.

fiz 168 exemplos de valores no delphi e gerei um TXT pra comparar as diferenças.

Onde soma de 5860amostras / 5860 amostras deu 15, o calculo em RMS deu 15,235...
Ou seja, toném aí.

Abração DJ, muito obrigado.

O pirex também vo ti contá viu.

O melhor mesmo é fazer o SQRT com o pic em 12 mips, demorou só 2.3 segundos.kkkk

Fabim

[Djalma Toledo Rodrigues]

Isso depende da comprexidade do Sinal.

Para DC 10 V e 3 amostragens

E(rms) = Sqr { (100 + 100 + 100) / 3 } = 10

Média Ariti = 10

Para Sinal complx

E(rms) = {( 100 + 0 + 100) / 3 } = ~ 8.165

Média Ariti (10 + 0 + 10) / 3 = ~ 6.667 dif. significativa
.

[fabim]

Isso depende da comprexidade do Sinal.

Para DC 10 V e 3 amostragem

E(rms) = Sqr { (100 + 100 + 100) / 3 } = 10

Média Ariti = 10

Para Sinal complx

E(rms) = {( 100 + 0 + 100) / 3 ) = ~ 8.165

Média Ariti (10 + 0 + 10) / 3 = ~ 6.667 dif. significativa
.

poisé, estava havendo um erro ilógico que eu não achava onde estava.

o osciloscopio me mostra V RMS 0.106mV.

O meu sistema esta me mostrando 56mV.

mechendo aqui agora sem querer eu selecionei a V média do osciloscopio. Esta tensão média é exatamente a diferença de erro do meu sistema para que o osciloscopio esta mostrando...

[proex]

Oh MANÉZÃO (Fabim), eu te passei a formula correta e te disse que nao da pra fazer isso com PIC.

Deixa de ser Zé Mané, pobre, miserável, faz 2 anos que estou te falando que PIC só serve pra medir tensão DC. Compra logo uma interface JTAG e aprende a mexer com ARM.

PIC nao vai te levar a lugar algum, vc só vai ganhar dinheiro com projeto quando abandonar PIC.

[fabim]

Poisé, esse ano depois de 9, eu vou pegar 30 dias de férias.
Minha cabeça esta com alguns dutos de raciocinio entupidos..

Fazendo 1000 coisas em paralelo, ja esta me prejudicando.

Vlw DJ, é pirex..
Eu até que gostaria, mais como eu ja te expliquei...................

Fabim

[fabim]

bom, o tchelo me deo uma idéia ontem a noite.

Alguém aí ja viu ou sabe de alguma função em C, que é a famosa fast root maxell ?
Se eu fizer 98 amostras por ciclo de 60hz, o valor da somatoria de 12 bits ao quadrado não estoura uma int32...
aí eu vou ficar fazendo o ping pong de ciclos, e vai dar tempo de ficar testanto a VP e a IP para a proteção em tempo real.
Achei uma rotina de root para int16, com o pic rodando a 12 mips, e eu jogando o valor 65535, ele demorou 56uS para retornar a raiz inteira.

Código: Selecionar todos

 
 unsigned int fast_root16(unsigned int x){
unsigned int a,b;
b = x;
a = x = 0x3f;
x = b/x;
a = x = (x+a)>>1;
x = b/x;
a = x = (x+a)>>1;
x = b/x;
x = (x+a)>>1;
return(x);
}


Fiquei ontem o dia inteiro fazendo calculos e tentando ganhar tempo minha pressão foi 18-12... Pior foi quando eu tentei imaginar uma forma de usar um calculo parecido com esse aí pra int32...

Alguém conhece ou ja viu algum root pra int32 ?

Abraços

Fabim

[fabim]

http://www.codecodex.com/wiki/index.php?title=Calculate_an_integer_square_root

acabei de achar este.
só estou dando uma escovada.

Abraços

Fabim

[Djalma Toledo Rodrigues]

O Sinal Analógico ao ser Atenuado não perde razoavelmente suas demais caracteristicas

( Subentada-se um bom Atenuador, isto é, Atenuador Compensado )

para que não ocorra Overflow no Somatório das amostragens ao quadrado, passe o Sinal Analógico por um atenuador

divida-o por 2, por 4, ou por 8 ...
Pode inclusive dividir também após o ADC mas, aqui ocorrerá perdas a menos que o uC tenha Ponto Flutuante, como os da Renesas 32. Depois de calculada a respectiva E(rms) é só multiplicar pelo mesmo fator
( por isso, para facilitar, coloquei o Atenuador em Potência de 2 e não em decimal )


Para extrair Raiz Quadrada:

Veja o Algoritmo

Teorema dos quadrados perfeitos: "todos os quadrados perfeitos são formados pela somatória dos n primeiros números ímpares"

Ex.:

sqr (4) => 4 = 1 + 3 => Duas parcelas sqr = 2
sqr (9) => 9 = 1 + 3 + 5 => Três parcelas sqr = 3

A partir dai você pode calcular facil na mão quantas casas decimais voce quiser da raiz de qq número


Algoritmo:

N(ímpar) = 1
Prcl = 0 ; Quantidade de parcelas
SN = 0 ; Somatório dos N(impar)

Roda até que SN > Sqr (N) ; Maior para o caso em não seja quadrado perfeito
Incr N(ímpar)
Incr N(ímpar) ; Próximo número ímpar = N(ímpar) + 2 --- Claro!
SN = SN + N(ìmpar)
Inc Prcl
Roda
Sqr (N) = Prcl - 1

Algo assim ...
.

[proex]

O que o Fabim precisa na realidade é medir o True Rms, e não o Rms.

Já mostrei pra ele como fazer mas ele insiste em usar PIC pra fazer isso.

Assim não vai chegar a lugar nenhum, no maximo até o posto de saúde mais proximo ehehhehehe.

.

[MOR_AL]

Olá Pessoal!!! Tudo bem?

Escrevi alguma coisa sobre valor eficaz. Acho que auxilia na compreenção do que ocorre. Deixei uma pergunta no final, por falta de conhecimento mesmo. Caso alguém saiba a resposta, agradeço se for respondida.

http://rapidshare.com/files/260262034/Valor_Eficaz.pdf

Espero que seja útil.
Abraços.
MOR_AL

[fabim]

intão pirex.. Vou fazer com LPC mesmo.
Como eu te dise na 6°F, é a porcaria do tekronix que fica me mostrando aquele valor.
Eu tenho certeza absoluta que o meu valor esta correto, porêm o tek diz que não.
amanhã de manhãzinha ja vou falar com o so augusto pro tchelo trazer o murtimetro preu.

Esse aqui ó.
http://www.protektest.com/ProdPDF2006/2006_25_Pg23_6500.pdf

Pelomenos assim não vou estar comparando tensão rms de um osciloscopio com 500msps, e que só mostra a tensao RMS XXX "?"..

fuizme

[Djalma Toledo Rodrigues]

Mor_al
É obvio que é impossível calcular uma Integral de mais infinito a menos infinito.

Nesses casos se abre um "Janela de Tempo" e se expecifica como uma unidade de tempo. Por segundo, por minuto, por hora etc. caso do kW/h, km/h, rpm, etc.

Nada impede que calcule a RMS de uma Tensão Contínua. Faço N amostragens em uma unidade de Tempo.

Caso do Circuito do Fabim, se não houver capacitor de acoplamento.
.

[Djalma Toledo Rodrigues]

... trazer o murtimetro preu. Esse aqui ó.

http://www.protektest.com/ProdPDF2006/2006_25_Pg23_6500.pdf

fuizme

Bom hein ? 4 a 1200 Ohms no dBm e + RSR232.
.