Continue to Site

Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronics Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

Help me solve a symbolic polynomial with complex coefficient

Status
Not open for further replies.

aibrahim

Junior Member level 2
Junior Member level 2
Joined
Aug 14, 2006
Messages
23
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,429
Dear All,

I want to solve a symbolic polynomial with complex cofficient and find its roots in terms of symboles

example:

(((s-I*w)^2+g1)*((s-I*w)^2+g2)+(a1+a2*(s-I*w))+(a3+a4*(s-I*w))*((s-I*w)^2+g1))

Where I=√-1 and g1,g2,w,a1,a2,a3,a4 are symboles

and i want to solve this equation in s

i try maple functions but i don't reach roots of the polynomial in terms of this symboles if any can help me through suggesting another function in maple or matlab i will be greatful.

Regards,
 

Re: symbolic polynomial

Dear friend

You can use MATLAB to solve it:

syms g1 g2 w a1 a2 a3 a4 ;
R = solve('(((s-i*w)^2+g1)*((s-i*w)^2+g2)+(a1+a2*(s-i*w))+(a3+a4*(s-i*w))*((s-i*w)^2+g1))','s')



Regards
 

Re: symbolic polynomial

Thanks alot,

I try it but the solution is no longer simplified and i try simplify and pretty function in matlab but still solution not simplified can you suggest me away to get readable form for the solution.

Regards,
 

Re: symbolic polynomial

Dear friend

You have 7 unknown variables in your equation, so the roots can not have very simple terms.

In this program:
syms g1 g2 w a1 a2 a3 a4 ;
R = solve('(((s-i*w)^2+g1)*((s-i*w)^2+g2)+(a1+a2*(s-i*w))+(a3+a4*(s-i*w))*((s-i*w)^2+g1))','s');


R : ( size= 4x1 and Bytes=2141680).

Using expand function:

K = expand(R);

K : ( size= 4x1 and Bytes=1397784).

I think to simplify more, you should substitute some of unknown variables with known values. For example:

N=subs(K,{a1,a2,a3,a4},{'1', '2', '3', '4'});
L=simplify(N);
H=expand(L);


H : ( size= 4x1 and Bytes=96568).

Best
 

Status
Not open for further replies.

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Back
Top