Please help me solve this four equations with methode newton-raphson urgent (i have MATLAB R2010a)
Show older comments
hi ,i need file fonction and the code matlab to solve those 4 equations [cos(3.θ1)+cos(3.θ2)+cos(3.θ3)+cos(3.θ4) ]=0 ; [cos(5.θ1)+cos(5.θ2)+cos(5.θ3)+cos(5.θ4) ]=0 ; [cos(7.θ1)+cos(7.θ2)+cos(7.θ3)+cos(7.θ4) ]=0 ; [cos(9.θ1)+cos(9.θ2)+cos(9.θ3)+cos(9.θ4) ]=0 ; with first condition θ1=0.314 ; θ2=0.628 ; θ3=0.945 ; θ4=1.225 ; in interval [0 pi/2] ; the angel's with radiane not degrees thnx
4 Comments
John D'Errico
on 7 May 2016
Edited: John D'Errico
on 7 May 2016
If you need code, then you should write code. We are not a code writing service. You won't learn to use MATLAB until you bother to use it.
Start with the symbolic toolbox, so solve is the tool for this problem. If it fails to give you something useful, then use vpasolve.
You might also try the function fsolve, from the optimization toolbox.
Of course, there are many poorly written Newton-Raphson tools on the file exchange. You might try one of them on your problem. In those cases, free code is probably worth just what you paid for it.
By the way, you need to learn basic MATLAB syntax too. . is not a multiplication operator. * is what is used for multiplication. So again, you need to start learning to use MATLAB.
yahia mehda
on 7 May 2016
Jan
on 7 May 2016
@yahia mehda: Please read this thread: Answers: Why your question is not urgent . You will see, that trying to push the voluntary helpers reduce their motivation to post an answer.
The strange notation with "θ1" is simply confusing. Posting the values directly would be a much better idea when you want others to assist you.
The forum usually ignore people, how do not show any own effort. Together with the "urgent" and the annoying "plllz" the chances to get an answer are really small. So better edit the question and show whatyou have tried so far, even if it does not work.
yahia mehda
on 7 May 2016
Answers (1)
Maneet Goyal
on 7 May 2016
Your problem requires the output to be in a certain range:
[0 2*pi]
I suggest using fmincon. However, the your equations can be put in a function as below:
function out = proby(theta)
out = zeros(4,1);
out(1) = cos(3*theta(1))+...);
out(2) = cos(5*theta(1))+...);
out(3) = cos(7*theta(1))+...);
out(4) = cos(9*theta(1))+...);
end
The fsolve command can be easily applied over such a function. Since you are new, try the following to learn about using fsolve.
help fsolve
3 Comments
yahia mehda
on 7 May 2016
yahia mehda
on 7 May 2016
Edited: Walter Roberson
on 7 May 2016
Maneet Goyal
on 12 May 2016
fsolve(@theta proby(theta), [initial_theta(1);initial_theta(2),...])
Try this!
Categories
Find more on Numeric Solvers in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!