Write a loop that can generate all data in one M-file with one variable change
1 次查看(过去 30 天)
显示 更早的评论
dt= a variable that is going to change.
% loop
q= 0
for i = 0:dt:10
q=q+1
for a= 1:3
for b = 1:3
if a~=b
then
.....calculate x, y ,
..... calculate h ( h depends on x, y)
.... hh(q)= h
.....tt(q) = i+dt;
I want to write the loop for the whole thing where it accounts dt is changing , and save data for each run then plot hh vs tt for any given dt change , dt can be [.1, .01, .001, .0001, .00001]. Know the loop to calculate h cant change. Anybody knows how to write loop as i describe above and still can plot hh vs tt @ dt = .1 , dt = .01 , etc in the same graph . Thanks
0 个评论
采纳的回答
Kevin Holst
2012-2-7
To answer your comments:
After looking back at the code, I made a big assumption that's not the best to make I think. I assumed that you're only interested in positive hh values. Try this code instead:
for j = -1:-1:-5
q = 0;
dt = 10^j;
for i = 0:dt:10
q=q+1
for a = 1:3
for b = 1:3
if a~=b
then
.....calculate x, y ,
.....calculate h ( h depends on x, y)
.... hh.(['dt' num2str(abs(j))])(q) = h
.....tt.(['dt' num2str(abs(j))])(q) = i+dt;
This will create two structures, hh and tt, with fields dt1, dt2, etc. To access the data you can plot it as:
plot(tt.dt1,hh.dt1,tt.dt2,hh.dt2)
You mentioned that it's only plotting the first column of data, is it actually only plotting the first column, or are the two sets of data overlapping exactly? This function appears to be a parametric study of dt, studying how dt size affects the plot. I would expect the plots to be on top of one another.
3 个评论
Kevin Holst
2012-2-7
I'm not sure I understand what you're asking there. When dt=0.1 data is collected at every 1 step, dt=0.01 data is collected at every 10th step? Is that correct?
更多回答(1 个)
Kevin Holst
2012-2-6
You could just add another dimension to your hh and tt calculations like this:
for j = -1:-1:-5
q = 0;
dt = 10^j;
for i = 0:dt:10
q=q+1
for a = 1:3
for b = 1:3
if a~=b
then
.....calculate x, y ,
..... calculate h ( h depends on x, y)
.... hh(q,abs(j))= h
.....tt(q,abs(j)) = i+dt;
Your hh and tt matrices would have a lot of zeros that you wouldn't want to plot, so your plot functions would need to be something like:
plot(tt(tt(:,1)>0),hh(tt(:,1)>0))
I don't have access to matlab right now so I can't say for 100% that that plot call will work as expected, but hopefully that'll point you in the right direction.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!