Ir para o conteúdo
  • Revista PROGRAMAR: Já está disponível a edição #58 da revista programar. Faz já o download aqui!

ssofia

Duvida de implementaçao ODE equaçao de van der pol

Mensagens Recomendadas

ssofia    0
ssofia

Ola a todos

Eu estou a tentar implementar uns scripts em matlab...e esta sempre a darM erro....e nao estou a perceber o que é, quer dizer é na equaçao mas nao estou a conseguir resolver....

Ou seja estou a implementar a a eq. Van der Pol para exportar 1 grafico com o comportamento da equaçoes - 4 equacoes ....

o A que representa a amplitude é um array com diferentes valores de amplitude medidos em 6 tempos para cada Oscilador.... o wt representa a frequencia natural que neste caso é 24h....

Espero que alguém me possa ajudar

Obrigado

___________________________________________________________________________________________

tspan=[0 100];
y0=[1 0 1 1];

global A1 A2 A3 wt

A1 = [0.000413401 , 0.0020457, 0.018875653, 0.008827253, 0.001098447, 0.00034614];
A2 = [0.005770291, 0.014026409, 0.01266695, 0.0032231, 0.000410485, 0.000287991];
A3 = [0.007822632, 0.005662967, 0.000258098, 0.000917877, 0.004825493, 0.007161058];

wt = 24;

[T,Y] = ode15s(@Month_3,[0 100],[1 0 1 1]);

plot(t,y(:,1),'-',t,y(:,2),'--',t,y(:,3),'.',t,y(:,4),'..')

title('Comportamento da Amplitude ao longo tempo');
xlabel('t');
ylabel('Amplitude');
legend('X','y(2)','Y','Z')

____________________________________________________________________________________________________

function dy = Month_3(t,y,A1,A2,A3,wt)


dy = zeros(4,1);  

dy(1) = y(2);

dy(2) = 0.03*(1 - y(1)^2)*y(2) - y(1) - A1 * sin(wt); % Comportamento Amplitude X

dy(3) = 0.03*(1 - y(3)^2)*y(2) - y(3) - A2 * sin(wt); % Comportamento Amplitude Y

dy(4) = 0.03*(1 - y(4)^2)*y(2) - y(4) - A3 * sin(wt); % Comportamento Amplitude Z

Editado por Rui Carlos
Formatação do código.

Partilhar esta mensagem


Link para a mensagem
Partilhar noutros sites

Crie uma conta ou ligue-se para comentar

Só membros podem comentar

Criar nova conta

Registe para ter uma conta na nossa comunidade. É fácil!

Registar nova conta

Entra

Já tem conta? Inicie sessão aqui.

Entrar Agora


×

Aviso Sobre Cookies

Ao usar este site você aceita os nossos Termos de Uso e Política de Privacidade. Este site usa cookies para disponibilizar funcionalidades personalizadas. Para mais informações visite esta página.