MATLAB-finding roots and summations

  • Thread starter Thread starter Dsab123
  • Start date Start date
  • Tags Tags
    Roots
AI Thread Summary
The discussion revolves around solving the polynomial equation p(z) = z^4 + 4.25*z^2 + 1 using MATLAB to find its roots. The user successfully computes the roots using the 'solve' function, obtaining four complex roots. However, they struggle with the second part of the assignment, which involves defining a complex-valued function Z(t) as a summation of the roots multiplied by the exponential of each root times t. Clarification is sought on how to implement the summation in MATLAB, particularly regarding the use of a for loop and the meaning of 'complex constants ck.' The conversation highlights the user's challenges with MATLAB's symbolic computations and summation techniques.
Dsab123
Messages
2
Reaction score
0

Homework Statement



(all this is included in the attachment, btw)

p(z)=z^4+4.25*z^2+1

1. use MATLAB to find the four roots Ck of
p(z)= 0

2. For all complex constants ck (different than Ck) define the complex valued function
Z(t)= 4\Sigmak=1 ck*exp(Ck*t)

^^that is a sigma sign with variable k going from one to four


Homework Equations



p(z)=z^4+4.25*z^2+1
a=4.25

The Attempt at a Solution



Number 1 I can do fine, my code is as follows:

>> syms x
>> y=(x^4+4.25*x^2+1)
>> >> solve (y)
ans =
(-2)*i
2*i
-i/2
i/2

The problem is when I get to problem 2, because I have just about no idea of how to compute a sum in matlab. I've tried a for loop with k going from one to for, after naming each of the roots C1, C2, C3, C4.

Also, I'm having some trouble figuring out what it means by 'complex constants ck'..
sorry I'm not much help in understanding this haha.


3 and 4 i think I can figure out
Any help would be very much appreciated
 

Attachments

Physics news on Phys.org
I don't do much with symbolic in Matlab, but something like this should work.

>> syms x
>> y=(x^4+4.25*x^2+1)
>> roots = solve (y)
Z=0;
syms t
for j=1:length(roots)
Z = Z + roots(j)*exp(roots(j)*t)
end
 

Similar threads

Back
Top