How to solve the equation x:
x = sum of k from -10 to 10 (1/4*k * ((cos(2*k*(t/T))-1)^2 + (sin(2*k*(t/T))))^2)
where t and T are unknown and plot the answer.
MATLAB: How to sum unknowns and plot
sum
Related Solutions
What is A?
Try this:
function x = inf_sum(t,A)k = 1:200;[k,t] = meshgrid(k,t);x = ((-1).^(k-1).*((2*A)./(k*pi)).*sin(k.*pi.*t));x = sum(x,2);Then from the command-line you do this:
A = 4;t = 0:.1:10;plot(t,inf_sum(t,A))I have a feeling, however, that your A is a Fourier Coefficient, thought you seem to leave no confirmation of this. If that is true you will also need to calculate A in the function.
No, you cannot "fit a curve" to this. At least not anything remotely useful, nor do you need to do so.
That a solution is generated using rootof is not bad actually. Just a statement that no analytical solution was found in the form of roots.
You can use vpa to resolve the rootof forms into numbers.
As an addendum, since this problem seems to be hanging around with no resolution.
Again, this is not the purpose of curve fitting, nor can it easily be used for that purpose. Lets see if we can understand SOME of the problems. To do that, I'll need to look in more depth at your problem.
syms a b k tE(1) = tan(k)==cos(0.43)*cos(b)/(cos(0.43)*sin(b)*sin(a)+sin(0.43)*cos(a));E(2) = tan(t)*sin(k)==cos(0.43)*cos(b)/(cos(0.43)*sin(b)*cos(a)-sin(0.43)*sin(a));pretty(vpa(E(1),5)) 0.90897 cos(b)tan(k) == -------------------------------------- 0.41687 cos(a) + 0.90897 sin(a) sin(b)pretty(vpa(E(2),5)) 0.90897 cos(b)sin(k) tan(t) == - -------------------------------------- 0.41687 sin(a) - 0.90897 cos(a) sin(b)So there WILL be infinitely many solutions, if any solution exists. But if we restrict the domain of a and b to lie within [0,2*pi], can we do anything?
Even that is difficult, since there may STILL be multiple solutions.
Lets pick some arbitrary value of k and t, just to see what is happening. Remember that k and t should also be expected to live in ranges. So lets say that
k lives in the interval [-pi,pi]t lives in the interval [-pi/2,pi/2]Create a couple of functions.
fun1 = @(a,b) cos(0.43)*cos(b)./(cos(0.43)*sin(b).*sin(a)+sin(0.43)*cos(a));fun2 = @(a,b) cos(0.43)*cos(b)./(cos(0.43)*sin(b).*cos(a)-sin(0.43)*sin(a));And a gridded lattice for a and b.
[A,B] = meshgrid(linspace(0,2*pi,100));f1AB = fun1(A,B);f2AB = fun2(A,B);Now, lets do some plots that will be HIGHLY informative.
contour(A,B,f1AB,[1 1]*tan(k),'b')hold oncontour(A,B,f2AB,[1 1]*sin(k)*tan(t),'r')grid onxlabel 'A'ylabel 'B'
So in blue, we see the solution locus for equation 1. All points on those curves in clue are the solutions of this problem for that equation.
Likewise, all points on those curves in red form the solution locus for equation 2. Any intersection of those curves is a solution to BOTH equations, and therefore an entirely valid solution to your problem.
This is only approximate, because contour was used to solve the problem, on a grid of limited size. But there are clearly something like 12 distinct solutions, where those curves cross.
Now, you MIGHT be able to decide that some of those solutions might be excluded, trying to find a fundamental solution of some sort. However, it is entirely possible that some other choices of t and k might have a different number of solutions. (By the way, I am quite confidant that I can find sets of k and t that have a different number of solutions.) Anyway, which of those solutions is somehow more fundamental than the others?
One thing I can tell you with certainty: CAN you use curve fitting to resolve this problem? NO WAY. PERIOD.
Worse, does there even exist a function g(k,t) that yields the "solution"? Again, this is a multi-valued function. It looks like there will be something like a dozen branches to that function, a dozen possible solutions. So no, you won't find some simple function, positively NOT using curve fitting.
If you were willing to use a table lookup with interpolation, you might consider spending the time to solve this problem for hundreds of values of k and t, isolating the specific pair for a and b that make you happy. Then create a table over a complete grid for k and t in the ranges I showed above. You might then try to use a tool like interp2 to interpolate. Even that approach can be fraught with problems, since you would need to be incredibly careful to find a consistent choice among the dozen (or more or less) solutions that may exist.
Other ideas might revolve around series approximations. But remember that since tan(k) and tan(t) can take on ANY values from -inf to inf, series approximations of any viability will be difficult as hell to find.
Don't think any of this will be easy to do.
Best Answer