How can I get this function to stop graphing a thousand straight lines and instead graph the data points it's supposed to?

1 次查看（过去 30 天）
tk27182
tk27182 2014-12-9
评论： Chad Greene 2014-12-9
When I call the following function, it should give me two nice bell curves. However, it's giving me a thousand straight lines on a single graph. How do I fix it so that I get the data points?
Here is the section of code:
function xprime = dcoupledeta(t,x)
c1 = 1;
c2 = 1;
E = .001;
m = 5;
k = 1;
s=size(t);
eta1 = zeros(1000,1);
eta2 = zeros(1000,1);
P = zeros(1,1000);
P2 = zeros(1,1000);
for j=1:s
for i=1:s
eta1(i) = .0001*i;
eta2(i) = .001*i;
xprime(1)= x(3);
xprime(2)= x(4);
xprime(3)= -eta1(i)*(((x(3)^2)/2) + (c1/(2*m))*x(1)^2 -(E/m))*x(3) - (c1/m)*x(1) + k*(x(2)-x(1));
xprime(4)= -(eta2(i)/m)*x(4) - (c2/m)*x(2) + k*(x(1)-x(2));
xprime = xprime(:);
end
P(j) = (1/100)*eta1(i)*x(4,j)^2;
P2(j) = (1/100)*eta2(i)*x(4,j)^2;
end
figure(j);
hold on;
plot(eta1,P,eta2,P2);
xlabel('\eta')
ylabel('Power')
title('Power vs. \eta')
end

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

回答（1 个）

Chad Greene
Chad Greene 2014-12-9
What are your x and t values?
  2 个评论
Chad Greene
Chad Greene 2014-12-9
We need some numbers to work with, and the ode45 line throws an error without tspan. From your question, I have some suspicions about what the problem might be, but I need to be able to replicate the problem to be sure. What values of x and t are you putting into dcoupledeta(t,x)?

请先登录，再进行评论。

类别

在 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