non-integer power for a model ?

8 次查看（过去 30 天）
Mohamed
Mohamed 2013-12-31
回答： Image Analyst 2014-1-2
I wanna to make a program that generate a CRONE regulation system, and it use fractional power for models and I can't find how to do it
  6 个评论
显示 4更早的评论
Image Analyst
Image Analyst 2014-1-1
Well g is gone, so that's good. I still don't know what is inside tf() since you didn't give the code for that function. And you still didn't answer my questions. I'm still not sure how to help. If you run your code, but fixed to use .^ instead of ^ you get this:
P=[1 3 6];
m = 1.4; %non integer power
P .^ m% (x²+3x+6)^m;
ans = 1 4.65553672174608 12.2860350664753
There are no unknowns. You have an m, you have a P, you raised P to the m power, end of story. What does this x in the comment have to do with anything???
Mohamed
Mohamed 2014-1-2
well P is a polynomial P= x^2+3x+6 and I want to raise it to the m power with m=1.4; if m is an integer (m=3) P raised to the 3 power is P*P*P in matlab conv(P,conv(P,P)) , but I don't know how to raise a polynomila to a non integer power i'm sorry for my english and thnx for your help

请先登录，再进行评论。

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

回答（2 个）

Image Analyst
Image Analyst 2013-12-31
So you want to have a model like y = a * x .^ n. Do you know the fractional power, n, or is it an unknown? Either way, one way is to take the log of the equation which turns it into a linear equation - that's pretty common and works reasonably well, though I'm not sure if it gives you the absolute minimum sum of residuals as compared to some other non-linear method.
log(y) = log(a) + n * log(x)
yNew = aNew + n * xNew
Now it's a linear equation in the unknowns aNew and n and you can just use least squares.
Image Analyst
Image Analyst 2014-1-2
You're doing it right. Just assign some x over some range you care about. Then construct P and finally raise it to a power, for every P (every x value) that you have in your array:
x = linspace(3, 25, 50); % Some arbitrary x values that we want.
P = x .^ 2 + 3 * x + 6;
m = 1.4; % non integer power
output = P .^ m;
plot(x, output, 'bo-');
grid on;
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
Is that what you want?

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by