How to Solve Complex PDEs and Calculate Wiener Filter Using Matlab PDE Toolbox?

  • Context: MATLAB 
  • Thread starter Thread starter RandomGuy88
  • Start date Start date
  • Tags Tags
    Gui Matlab Pde
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 1K views
RandomGuy88
Messages
404
Reaction score
6
I am attempting to solve the following PDE using the GUI for Matlab's PDE toolbox.

[tex] <br /> \newcommand{\pd}[3]{ \frac{ \partial^{#3}{#1} }{ \partial {#2}^{#3} } }<br /> <br /> \pd{\Psi}{y}{} <br /> + \pd{\Psi}{x}{2} + \pd{\Psi}{y}{2}=0 <br /> [/tex]

Is this possible? I have been able to use the PDE toolbox for other simpler PDEs, for example Laplace's Equation with the same boundary conditions I am using for the above equation. But I can't seem to get it to work once I add first partial of Psi w.r.t y

Does anyone know how I can do this?

Thanks.
 
Physics news on Phys.org
Hi all!

Problem:
I am currently trying to calculate a Wiener filter for a stochastic system. The model is an ARX with determined parameters.

Where I am:
I have access to the transfer functions of the ARX model.
I need to calculate the optimal causal filter:
H(s)=1/Fi(s)*[Sx(s)/Fi(-s)]L

I know that:
Fi(s)Fi(-s)=Sx(s)+Sn(s)
where Sx(s) is the power spectral density of the output signal and Sn(s) is the power spectral density of the aditive noise. To find Sx and Sn I found the square root of the absolute value of the transfer functions of the model and the noise filter respectively, in the jω domain.

I have the whitening filter (1/Fi(s)), I determined Fi(s) by taking out all poles and zeros on the right plane of Fi(s)Fi(-s).

Question:
Is the bode diagram of the whitening filter supposed to be the symetric, relative to magnitude, of the bode diagram of the noise filter of the model?

Now I have to determine:
[Sx(s)/Fi(-s)]L
I have Sx(s) and Fi(-s), but my question is how do I determine the transfer function of the non-causal part only? I know I could use partial fraction expansion by hand but my Sx(s) and Sn(s) have 12th order polynomials so I will certainly not go that way.

Please help.
Thank you.

Gonçalo