How to plot values from a for loop that don't start at x=0?
1 回表示 (過去 30 日間)
古いコメントを表示
I am trying to plot values from multiple for loops for different values of x that are greater than 0. I have one for loop from x=100:215 and x=215:230. When I try to plot them in the same plot the two lines are not connected as I would like them to be and they both have a line connecting to the origin. Is there anyway to plot them without the lines connecting to the origin and the two lines connected at x=215
1 件のコメント
Jan
2016 年 3 月 5 日
Without seeing your code, it is hard to guess, why the lines have a connection to the origin.
回答 (1 件)
Jan
2016 年 3 月 5 日
I do not see the reason for any problem. Perhaps this helps:
x1 = 100:215;
x2 = 215:230;
y1 = rand(size(x1));
y2 = rand(size(x2)) + 2;
plot([x1, x2], [y1, y2]);
3 件のコメント
Jan
2016 年 3 月 5 日
Please, m13, give us a chance to understand your problem and post the code, which reproduces the effect. We cannot guess what you have written.
Shannon Hemp
2016 年 3 月 13 日
Here is the code that I have. I want to plot time vs thrust.
dt_1=0; dt_2=0;dt_3=0;dt_4=0;
for t=0:100; Pc2_RD_1(t+1) = Pc_RD*(1+(.015*dt_1)^2); Pe_RD_1(t+1) = Pc2_RD_1(t+1)*(1+((gamma_RD-1)/2)*Me_RD^2)^(gamma_RD/(1-gamma_RD)); Isp_RD_1(t+1) = (lamda_RD*c_RD/g0)*(gamma_RD*sqrt((2/(gamma_RD-1))*((2/(gamma_RD+1))^((gamma_RD+1)/(gamma_RD-1)))*(1-((Pe_RD_1(t+1)/Pc2_RD_1(t+1))^((gamma_RD-1)/gamma_RD))))+(e_RD/Pc2_RD_1(t+1))*(Pe_RD_1(t+1)-Pa)); mprop_RD1(t+1) = At_RD*Pc2_RD_1(t+1)/c_RD; thrust_RD1(t+1) = Isp_RD_1(t+1)*g0*mprop_RD1(t+1)/1000; dt_1=dt_1+.2; time(t+1)=t+.2; end
for t2=100:215; Pc2_RD_2(t2) = Pc_RD*(.95+(.008*dt_2)^2); Pe_RD_2(t2) = Pc2_RD_2(t2)*(1+((gamma_RD-1)/2)*Me_RD^2)^(gamma_RD/(1-gamma_RD)); Isp_RD_2(t2) = (lamda_RD*c_RD/g0)*(gamma_RD*sqrt((2/(gamma_RD-1))*((2/(gamma_RD+1))^((gamma_RD+1)/(gamma_RD-1)))*(1-((Pe_RD_2(t2)/Pc2_RD_2(t2))^((gamma_RD-1)/gamma_RD))))+(e_RD/Pc2_RD_2(t2))*(Pe_RD_2(t2)-Pa)); mprop_RD2(t2) = At_RD*Pc2_RD_2(t2)/c_RD; thrust_RD2(t2) = Isp_RD_2(t2)*g0*mprop_RD2(t2)/1000; dt_2=dt_2+.2; time2(t2)=t2+.2; end
for t3=215:230; Pc2_RD_3(t3+1) = Pc_RD*(.98-(.04*dt_3)); Pe_RD_3(t3+1) = Pc2_RD_3(t3+1)*(1+((gamma_RD-1)/2)*Me_RD^2)^(gamma_RD/(1-gamma_RD)); Isp_RD_3(t3+1) = (lamda_RD*c_RD/g0)*(gamma_RD*sqrt((2/(gamma_RD-1))*((2/(gamma_RD+1))^((gamma_RD+1)/(gamma_RD-1)))*(1-((Pe_RD_3(t3+1)/Pc2_RD_3(t3+1))^((gamma_RD-1)/gamma_RD))))+(e_RD/Pc2_RD_3(t3+1))*(Pe_RD_3(t3+1)-Pa)); mprop_RD3(t3+1) = At_RD*Pc2_RD_3(t3+1)/c_RD; thrust_RD3(t3+1) = Isp_RD_3(t3+1)*g0*mprop_RD3(t3+1)/1000; dt_3=dt_3+.2; time3(t3+1)=t3+.2; end
for t4=230:238; Pc2_RD_4(t4+1) = Pc_RD*(.75-(.05*dt_4)); Pe_RD_4(t4+1) = Pc2_RD_4(t4+1)*(1+((gamma_RD-1)/2)*Me_RD^2)^(gamma_RD/(1-gamma_RD)); Isp_RD_4(t4+1) = (lamda_RD*c_RD/g0)*(gamma_RD*sqrt((2/(gamma_RD-1))*((2/(gamma_RD+1))^((gamma_RD+1)/(gamma_RD-1)))*(1-((Pe_RD_4(t4+1)/Pc2_RD_4(t4+1))^((gamma_RD-1)/gamma_RD))))+(e_RD/Pc2_RD_4(t4+1))*(Pe_RD_4(t4+1)-Pa)); mprop_RD4(t4+1) = At_RD*Pc2_RD_4(t4+1)/c_RD; thrust_RD4(t4+1) = Isp_RD_4(t4+1)*g0*mprop_RD4(t4+1)/1000; dt_4=dt_4+.2; time4(t4+1)=t4+.2; end
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!