You are now following this question
- You will see updates in your followed content feed.
- You may receive emails, depending on your communication preferences.
How to solve this non-linear equation?
2 views (last 30 days)
Show older comments
Can you please provide a code for this equation? k is in the range of 3<=k<=16, and yj[k] is a known value at given time,k. S[k] is a sequence of binary matrix that have different matrix when k is changing, I am not sure how to code this please help me!
Answers (1)
Walter Roberson
on 15 May 2020
Which of those variables are vectors or matrices? would normally indicate a transpose but if is a scalar then it does not make any sense to transpose it. I can only make sense out of it if is a vector being transposed.
Which is the variable being solved for?
13 Comments
Eric Chua
on 15 May 2020
Hi, thanks for replying my question. From the equation, S[k] is a vector , or I should say it is a 7x1 row matrix of binary sequence where each of the elements in the vector contain only 1 or 0.
yj[k] is a value that i have calculated out. For example, yj[3]= 141 (number of molecules observed at time k=3. So my question here is to estimate a set of C values, so there will be fourteen values of C at each time interval. (k=3 until k=16). I hope my explanation is clear, I am not sure how to write the code because it involve a lot of matrix and i need your help please.
Walter Roberson
on 15 May 2020
What you describe sounds like S[k] would be a scalar, with S being a vector? Or is S a 2D array and S[k] would be a row of it?
If there is one equation and 14 variables, it is greatly under-specified unless there are constraints on the C values.
I notice you have -- does that indicate complex conjugate? Are the C expected to be complex valued?
Eric Chua
on 16 May 2020
Hi, thanks for the follow up. S[k] is a vector. For example at time k=3, S[3] = [1 1 1 1 1 1 1]^T , at time =4, S[k] = [0 0 1 1 1 1 1]^T. The transpose here just for writing paper purpose, so S[k] is a 7x1 vector. yj[k] is a known value at specific time. For example at time =3, yj[3] = 141, at time k=4, yj[4] = 41. Cj is the 7x1 vector that I need to estimate.
In here, j means the receiver number. For example, y1[3]=141 would be the number of molecules observed at receiver 1 at time k=3. So given the value of S[k] and yj[3], I would like to estimate my Cj.
Walter Roberson
on 16 May 2020
Are there constraints such as C needing to be positive or integer valued? Without constraints you can solve exactly for one of the C values in terms of arbitrary choices for the other C values.
Eric Chua
on 16 May 2020
Thanks for follow up.
Cj are positive semidefinite. However maximum likelihood estimator could estimate negative value for some elements of Cj. Sub optimal solution is to set to zero all the negative entries of the estimated C. (This sentence is the explanation from the paper).
I think Cj will come out as a 7x1 vector.
Walter Roberson
on 21 May 2020
"The ML estimate of the CIR for the D-MIMO channel at receiver j is obtained by solving a system of non-linear equations given below [14]:"
Notice that they are solving systems of such equations, which means more than one such equation, each of which has that form. You cannot solve just one at a time.
Walter Roberson
on 21 May 2020
Edited: Walter Roberson
on 21 May 2020
I am having difficulty figuring out what the proper sizes are for the variables. I see that the paper discusses the notation in more detail; I will try again later to understand it when my head is a bit clearer.
But the subtraction of implies that the thing being subtracted from is also a vector, and that you would get a vector of results, and you would add a number of those together, getting a vector of results. But you want that summation, which would be a vector, to equal 0, which is a mismatch in sizes (vector compared to scalar.)
See Also
Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!An Error Occurred
Unable to complete the action because of changes made to the page. Reload the page to see its updated state.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list
How to Get Best Site Performance
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Asia Pacific
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)