Determining several unknown variables in an equation using measured data.

Question

John Carroll 2021-6-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/857220-determining-several-unknown-variables-in-an-equation-using-measured-data

评论： John Carroll 2021-6-18

I am trying to calculate several unknown variable in an equation. I have measured data that is to fit an equation and I want to calculate the unknown variables. The data is measured S-Parameter data (y axis) and frequency (x axis). This data is converted to impedance and graphed with frequency on the x axis and the impedance on the y axis. This data should fit the following equation

Zd(w) is the measured data. I know f[GHz], w is the converted frequency, I can calculated C, tano is also known and j in the imaginary number. Using this data and the known values I would like to use Matlab to calculate Rc, Rf and Ls.

I have read a number of examples for doing something similar though not the same. I have an average understanding of Matlab but I am no expert and many of the answers I've read assumed a certain level of programming knowledge that I do not have.

I was hoping someone could help get me started in the right direction. I've looked at things like the cftool and optimizer but I'm not knowlegable enough to know how to use them or if that is the correct tool for this case. I've also looked at things like the nonlinear curve fitting tool with no success. Any help would be greatly apprieated.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2021-6-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/857220-determining-several-unknown-variables-in-an-equation-using-measured-data#answer_725740

在 MATLAB Online 中打开

Create a function

function goodness = fitparameters(RcRfLs, C, f, tandelta, w, Zd)
  Rc = RcRfLs(1); Rf = RcRfLs(2); Ls = RcRfLs(3);  
  Zd_projected = Rc + Rf and so on
  goodness = sum( (Zd_projected - Zd).^2 );
end

Now in the function that sets up all of the input values you can create

%assign to C, f, tandelta, w, Zd before this part
obj = @(RcRfLs) fitparameters(RcRfLs, C, f, tandelta, w, Zd);

Now you can minimize obj using your favorite minimizer, such as fmincon() or fminsearch()

13 个评论
显示 11更早的评论隐藏 11更早的评论

John Carroll 2021-6-17

在 MATLAB Online 中打开

I think I'm following. So I start with the Zd, y12_0, freq_0, cap_0, w, fGHz and tandelta all of which I either measured or calculated before hand. Then I create the function

function goodness = fitparameters(RcRfLs, cap_0, fGHz, tandelta, w, Zd)
  Rc = RcRfLs(1); 
  Rf = RcRfLs(2); 
  Ls = RcRfLs(3);  
  Zd_projected = Rc+Rf*sqrt(fGHz)+1i*(w*Ls-(1/(w*cap_0.*(1-(1i*tandelta)))));
  goodness = sum(abs(Zd_projected - Zd).^2);
end

which should show me the difference between the measured data and the equation I'm trying to fit. The difference being the values of Rc, Rf and Ls.

Then I run

obj = @(RcRfLs) fitparameters(RcRfLs, cap_0, fGHz, tandelta, w, Zd);
cfs0 = [.25,.25,1e-11];
[best, fval] = fminsearch(obj, cfs0);

I set cfs0 = to my inital guess and it runs. However the results aren't really what I'm expecting. In my case Zd, y12_0, freq_0, cap_0, w, fGHz and tandelta are all 400x1 vectors. Should I set this up to evaluate the functions at a single point in the vector one at a time?

Thanks

Walter Roberson 2021-6-17

在 MATLAB Online 中打开

There are two different possibilities when you say that those are 400 x 1 vectors.

One possibility is that you are doing a multivariate fitting in which you are trying to find single Rc, Rf, Ls values that explain all of the variation. For example if you were trying to determine the specific heat capacity of something, then you would have a number of observations times in which each observation recorded time, temperature, pressure, whatever (acidity... whatever), and you combine all of the information to model the single factor(s) you are looking for that are constant across the experiment. Each time reading has a corresponding temperature, corresponding pressure, corresponding whatever.

For that situation, in which you have multiple variables each with a corresponding recording, then you just need to vectorize.

  Zd_projected = Rc + Rf .* sqrt(fGHz) + 1i .* (w .* Ls - (1./(w .* cap_0 .* (1 - (1i .* tandelta)))));

The information from all of the readings is combined by the sum() of the squares down to a single error.

The other possibility is that you are running 400 different experiments, each of which has a separate Rc, Rf, Ls. If that was what you were doing, you would program differently to evaluate at single points. However... remember that if you have N parameters, then you need at least N different sets of input in order to be able to fit the parameters, so if you were doing 400 different experiments, you would need 1200 readings.

400 is not divisible by 3, so we can deduce that you are not doing 133 1/3 different experiments, with 3 sets of readings each.

We cannot rule out the possibility that you are doing 100 different experiments with 4 different readings each.

How many different Rc, Rf, Ls sets are you expecting based on your 400 readings? The answer is the same as the number of different sets you would split your data into for separate calculations.

John Carroll 2021-6-17

I run one experiment at different voltages so cap_0 is 0 volts cap_5 would be 5 volts for this example I am just running the 0 volt. I collect just the S-parameters at each bias voltage. The S-parameters and frequency information are part of the file I am collecting. With that data I can calculate the Zd and y12_0 for each point so I would have a 400 element vector with the z and y information. Then I calculate the tandelta at each point again giving me a 400 element vector. W is just a conversion of the frequency. And finally I calculate the cap_0 from the y12_0 at each point. So it is really just one experiment for each bias voltage. In this case I'm am just trying look at the single experiment of the 0 volt experiment. I expect to get a single Rc, Rf and Ls but I wasn't sure if it made more sense to get a vector of each and find the average or just allow the optimizer to give me single values for each. The problem is that whe I run it on some test data I get results approx Rc = 300, Rf = -60 and Ls = 7e-13 and I'm expecting closer to Rc=0.5, Rf=0.025 and Ls=30e-11. Roughly. But none of them can be negative. I'm asssuming I have something wrong but I feel like I'm close.

Thanks