how to plot two graphs in one plot

6 ビュー (過去 30 日間)

HARIKRISHNA B YADULADODDI 2021 年 4 月 23 日

0
リンク

Direct link to this question

https://jp.mathworks.com/matlabcentral/answers/811530-how-to-plot-two-graphs-in-one-plot

コメント済み: Adam Danz 2021 年 4 月 27 日

採用された回答: Adam Danz

I have created the following plot but i want to reverse the y direction(i,e top left corner should start with 50 in y direction) but if i reverse the y direction the vectors will point in downward direction and i dont want the vectors to point in downward direction if i reverse the y direction kindly help me . .

5 件のコメント
表示非表示 4 件の古いコメント

Adam Danz 2021 年 4 月 24 日

50112.txt

HARIKRISHNA B YADULADODDI 's duplicate question moved here

The Y-axis ticks go like 0,5,10- - - - - 50

I want the plot to remain as it is, but flip the ticks so that they go like 50 45 40 - - - - - - 5,0.

Can this be done? please find the code and file for reference.

                                clearall;closeall;clc
                               
                                formatlong
                               
                                %--------Reading velocity data-------%
                               
                                fid=fopen('50112.txt');
                               
                                fulltext=textscan(fid,'%s','delimiter','\n');
                               
                                numberofgridpoints=length(fulltext{1});
                               
                                fclose(fid);
                               
                                n=numberofgridpoints/200;
                               
                                fid=fopen('50112.txt');
                               
                                datamatrix = fscanf (fid,'%f%f%f%f',[4,inf]);
                               
                                %-------stroing velocity-----------%
                               
                                fori=1:200
                               
                                c=i-1;
                               
                                forj=1:n
                               
                                d=(c*n)+j;
                               
                                x(j)=datamatrix(1,j);
                               
                                y(j)=datamatrix(2,j);
                               
                                u(i,j)=datamatrix(3,d);
                               
                                v(i,j)=datamatrix(4,d);
                               
                                end
                               
                                end
                               
                                x1(1:13,1:20)=0;
                               
                                y1(1:13,1:20)=0;
                               
                                a=-19;
                               
                                b=10;
                               
                                forj=1:20
                               
                                fori=1:13
                               
                                x1(i,j)=a;
                               
                                end
                               
                                a=a+2;
                               
                                end
                               
                                x1=reshape(x1,1,260)
                               
                                fori=1:13
                               
                                forj=1:20
                               
                                y1(i,j)=b;
                               
                                end
                               
                                b=b+1.92307692307;
                               
                                end
                               
                                y1=reshape(y1,1,260)
                               
                                %-----------Ensemble mean------------%
                               
                                addu(1:n)=0;
                               
                                addv(1:n)=0;
                               
                                forj=1:n
                               
                                fori=1:200
                               
                                addu(1,j)=addu(1,j)+u(i,j);
                               
                                addv(1,j)=addv(1,j)+v(i,j);
                               
                                end
                               
                                avgu(1,j)=addu(1,j)/200;
                               
                                avgv(1,j)=addv(1,j)/200;
                               
                                end
                               
                                %---------------------plot--------%
                               
                                figure('name','ENSEMBLE AVARAGE FOR 50112')
                               
                                rectangle('Position',[-26.75 0 53.55 50]);
                               
                                xlabel('x (mm)')
                               
                                ylabel('y (mm)')
                               
                                axisequal
                               
                                axis ([-25 25 0 50])
                               
                                set(gca,'ytick',0:10:50)
                               
                                axistight
                               
                                holdon
                               
                                rectangle('Position',[-26.75 40 53.5 5],'FaceColor',[0.502 0.502 0.502])
                               
                                text (-2.5,42.5,'Piston','Color','black','FontSize',14)
                               
                                y=yline(1.2,'r','TDC')
                               
                                y.LineWidth = 1;
                               
                                y.LabelHorizontalAlignment='center'
                               
                                y.LabelVerticalAlignment='middle'
                               
                                holdon
                               
                                magnitude=sqrt(avgu.^2+avgv.^2);
                               
                                magnitude=reshape(magnitude,[13,20]);
                               
                                [X,Y]=meshgrid(-19:2:19,10:1.9230769:33.0769230);
                               
                                contourf(X,Y,magnitude,100,'linestyle','none')
                               
                                holdon
                               
                                h=quiver(x1,y1,avgu,avgv,'color','black');
                               
                                set(gca,'ydir',“反向”)
                               
                                axisequal
                               
                                c = colorbar;
                               
                                c.Label.String ='V [m/s]'
                               
                                c.Label.FontSize =14;
                               
                                c.FontWeight ='normal'
                               
                                holdoff

output: [ see figure in main question ]

サインインしてコメントする。

サインインしてこの質問に回答する。

採用された回答

Adam Danz 2021 年 4 月 24 日

0
リンク

この回答への直接リンク

https://jp.mathworks.com/matlabcentral/answers/811530-how-to-plot-two-graphs-in-one-plot#answer_683855

編集済み:Adam Danz 2021 年 4 月 24 日

You're reversing the y-axis with this line in your code.

set(gca,'ydir',“反向”)

If you want the y-axis direction to be normal but leave the contour and quiver plots in their current orientation, that means the data are in the wrong order. Ideally you would alter the code so that the matrices and vectors are constructed in the correct order. Bu here's a way to flip the data as it is.

Remove the line that reverses the y-direction of the y-axis.
Flip the y-coordinates for the contour plot (see below)
Flip the y-coordinatesandthe y-vector-component of the quiver plots in polar coordinates (see below)

Step 2:

                              Y = flipud(Y);% <--- add this
                             
                              contourf(X,Y,magnitude,100,'linestyle','none')

Step 3:

                              fori=1:13
                             
                              forj=1:20
                             
                              y1(i,j)=b;
                             
                              end
                             
                              b=b+1.92307692307;
                             
                              end
                             
                              y1 = flipud(y1);% <--- add this
                             
                              y1=reshape(y1,1,260)
                             
                              。..
                             
                              。..
                             
                              。..
                             
                              [theta, radius] = cart2pol(avgu,avgv);% <--- add this
                             
                              [avgu, avgv] = pol2cart(-theta, radius);% <--- add this
                             
                              h=quiver(x1,y1,avgu,avgv,'color','black');

Result

6 件のコメント
表示非表示 5 件の古いコメント

HARIKRISHNA B YADULADODDI 2021 年 4 月 26 日

50112.txt

Thank you for your kindful information @Adam Danz 。我试图实施同样的事情code of mine but it not giving the correct plot. Please find below the code for reference and also the attached data file. Kindly help me for the same mentioned above problem . . Thank You @Adam Danz

                                   %------------CODE---------------------------------%
                                  
                                   clearall;closeall;clc
                                  
                                   formatlong;
                                  
                                   %--------Reading velocity data-------%
                                  
                                   fid=fopen('50112.txt');
                                  
                                   fulltext=textscan(fid,'%s','delimiter','\n');
                                  
                                   numberofgridpoints=length(fulltext{1});
                                  
                                   fclose(fid);
                                  
                                   n=numberofgridpoints/200;
                                  
                                   fid=fopen('50112.txt');
                                  
                                   datamatrix = fscanf (fid,'%f%f%f%f',[4,inf]);
                                  
                                   %-------stroing velocity-----------%
                                  
                                   fori=1:200
                                  
                                   c=i-1;
                                  
                                   forj=1:n
                                  
                                   d=(c*n)+j;
                                  
                                   x(j)=datamatrix(1,j);
                                  
                                   y(j)=datamatrix(2,j);
                                  
                                   u(i,j)=datamatrix(3,d);
                                  
                                   v(i,j)=datamatrix(4,d);
                                  
                                   end
                                  
                                   end
                                  
                                   %-----------Ensemble mean------------%
                                  
                                   addu(1:n)=0;
                                  
                                   addv(1:n)=0;
                                  
                                   forj=1:n
                                  
                                   fori=1:200
                                  
                                   addu(1,j)=addu(1,j)+u(i,j);
                                  
                                   addv(1,j)=addv(1,j)+v(i,j);
                                  
                                   end
                                  
                                   avgu(1,j)=addu(1,j)/200;
                                  
                                   avgv(1,j)=addv(1,j)/200;
                                  
                                   end
                                  
                                   %---------------------plot--------%
                                  
                                   figure('name','ENSEMBLE AVARAGE FOR 50112')
                                  
                                   rectangle('Position',[0 0 53.5 50]);
                                  
                                   xlabel('x (mm)')
                                  
                                   ylabel('y (mm)')
                                  
                                   axisequal
                                  
                                   axis ([-25 25 0 50])
                                  
                                   set(gca,'ytick',0:10:50)
                                  
                                   axistight
                                  
                                   holdon
                                  
                                   rectangle('Position',[0 40 53.5 5],'FaceColor',[0.502 0.502 0.502])
                                  
                                   text (23,42.5,'Piston','Color','black','FontSize',14)
                                  
                                   yv=yline(1.2,'r','TDC')
                                  
                                   yv.LineWidth = 1;
                                  
                                   yv.LabelHorizontalAlignment='center'
                                  
                                   yv.LabelVerticalAlignment='middle'
                                  
                                   holdon
                                  
                                   magnitude=sqrt(avgu.^2+avgv.^2);
                                  
                                   magnitude=reshape(magnitude,[13,20]);
                                  
                                   [X,Y]=meshgrid(11.461148:1.580848000000000:41.497260000000018, 9.975784000000001:1.564829000000000:28.753731999999996);
                                  
                                   contourf(X,Y,magnitude,100,'linestyle','none')
                                  
                                   holdon
                                  
                                   quiver(x,y,avgu,avgv,'color','black');
                                  
                                   set(gca,'ydir',“反向”);
                                  
                                   axisequal
                                  
                                   axison
                                  
                                   c = colorbar;
                                  
                                   c.Label.String ='V [m/s]'
                                  
                                   c.Label.FontSize =14;
                                  
                                   c.FontWeight ='normal'
                                  
                                   holdoff
                                  
                                   %--------------END OF ENSEMBLE AVAERAGE---------%
                                  
                                   Output:

HARIKRISHNA B YADULADODDI 2021 年 4 月 27 日

編集済み:HARIKRISHNA B YADULADODDI 2021 年 4 月 27 日

After your saying i edited the first mentioned code in question and i was able to get the required output but when i tried to do the same thing in above mentioned code in comment section i was not able to get it.. Kindly help me in getting the output. . @Adam Danz 。.... I have attached edited but i am not getting the plot after follwowing your instructions.

                                   clearall;closeall;clc
                                  
                                   formatlong;
                                  
                                   %--------Reading velocity data-------%
                                  
                                   fid=fopen('50112.dat');
                                  
                                   fulltext=textscan(fid,'%s','delimiter','\n');
                                  
                                   numberofgridpoints=length(fulltext{1});
                                  
                                   fclose(fid);
                                  
                                   n=numberofgridpoints/200;
                                  
                                   fid=fopen('50112.dat');
                                  
                                   datamatrix = fscanf (fid,'%f%f%f%f',[4,inf]);
                                  
                                   %-------stroing velocity-----------%
                                  
                                   fori=1:200
                                  
                                   c=i-1;
                                  
                                   forj=1:n
                                  
                                   d=(c*n)+j;
                                  
                                   x(j)=datamatrix(1,j);
                                  
                                   y(j)=datamatrix(2,j);
                                  
                                   u(i,j)=datamatrix(3,d);
                                  
                                   v(i,j)=datamatrix(4,d);
                                  
                                   end
                                  
                                   end
                                  
                                   y=flipud(y)
                                  
                                   y=reshape(y,1,260)
                                  
                                   %-----------Ensemble mean------------%
                                  
                                   addu(1:n)=0;
                                  
                                   addv(1:n)=0;
                                  
                                   forj=1:n
                                  
                                   fori=1:200
                                  
                                   addu(1,j)=addu(1,j)+u(i,j);
                                  
                                   addv(1,j)=addv(1,j)+v(i,j);
                                  
                                   end
                                  
                                   avgu(1,j)=addu(1,j)/200;
                                  
                                   avgv(1,j)=addv(1,j)/200;
                                  
                                   end
                                  
                                   %---------------------plot--------%
                                  
                                   figure('name','ENSEMBLE AVARAGE FOR 50112')
                                  
                                   rectangle('Position',[0 0 53.5 50]);
                                  
                                   xlabel('x (mm)')
                                  
                                   ylabel('y (mm)')
                                  
                                   axisequal
                                  
                                   axis ([-25 25 0 50])
                                  
                                   set(gca,'ytick',0:10:50)
                                  
                                   axistight
                                  
                                   holdon
                                  
                                   rectangle('Position',[0 40 53.5 5],'FaceColor',[0.502 0.502 0.502])
                                  
                                   text (23,42.5,'Piston','Color','black','FontSize',14)
                                  
                                   yv=yline(1.2,'r','TDC')
                                  
                                   yv.LineWidth = 1;
                                  
                                   yv.LabelHorizontalAlignment='center'
                                  
                                   yv.LabelVerticalAlignment='middle'
                                  
                                   holdon
                                  
                                   magnitude=sqrt(avgu.^2+avgv.^2);
                                  
                                   magnitude=reshape(magnitude,[13,20]);
                                  
                                   [X,Y]=meshgrid(11.461148:1.580848000000000:41.497260000000018, 9.975784000000001:1.564829000000000:28.753731999999996);
                                  
                                   Y = flipud (Y)
                                  
                                   contourf(X,Y,magnitude,100,'linestyle','none')
                                  
                                   holdon
                                  
                                   [theta, radius] = cart2pol(avgu,avgv);% <--- add this
                                  
                                   [avgu, avgv] = pol2cart(-theta, radius);% <--- add this
                                  
                                   quiver(x,y,avgu,avgv,'color','black');
                                  
                                   axisequal
                                  
                                   axison
                                  
                                   c = colorbar;
                                  
                                   c.Label.String ='V [m/s]'
                                  
                                   c.Label.FontSize =14;
                                  
                                   c.FontWeight ='normal'
                                  
                                   holdoff
                                  
                                   %--------------END OF ENSEMBLE AVAERAGE---------%

Output:

Adam Danz 2021 年 4 月 27 日

First , in future posts, please format your code using the [>] button (see help ). I've formatted your code above for you.

Second , when I run your code I get the error (below). Is '50112.dat' exactly the same as 50112.txt ?

                                   Indexin position 2 is invalid. Array indices must be positive integers or logical values.
                                  
                                   Errorin temp (line 17)
                                  
                                   u(i,j)=datamatrix(3,d);

Third , you've made lots of changes from the original code you shared 。My answer applies to the original code. You didn't add the lines in the correct place.

For example, look at step 2 in my answer, do you see the " <--- add this " ?

                                   Y = flipud(Y);% <--- add this
                                  
                                   contourf(X,Y,magnitude,100,'linestyle','none')

That means you're supposed to add it just before the contourf() line....

In step 3, there are three " <--- add this " . Make sure you're adding them in the right place. Try again.

                                   y1 = flipud(y1);% <--- add this
                                  
                                   y1=reshape(y1,1,260)
                                  
                                   。..
                                  
                                   。..
                                  
                                   。..
                                  
                                   [theta, radius] = cart2pol(avgu,avgv);% <--- add this
                                  
                                   [avgu, avgv] = pol2cart(-theta, radius);% <--- add this
                                  
                                   h=quiver(x1,y1,avgu,avgv,'color','black');

Adam Danz 2021 年 4 月 27 日

50112.txt

Great!!

Here's what I did with the attached data. Look for the 5 leftward arrows " % <------ "

%--------Reading velocity data-------%

fid=fopen('50112.txt');

fulltext=textscan(fid,'%s','delimiter','\n');

numberofgridpoints=length(fulltext{1});

fclose(fid);

n=numberofgridpoints/200;

fid=fopen('50112.txt');

datamatrix = fscanf (fid,'%f%f%f%f',[4,inf]);

%-------stroing velocity-----------%

fori=1:200

c=i-1;

forj=1:n

d=(c*n)+j;

x(j)=datamatrix(1,j);

y(j)=datamatrix(2,j);

u(i,j)=datamatrix(3,d);

v(i,j)=datamatrix(4,d);

end

x1(1:13,1:20)=0;

y1(1:13,1:20)=0;

a=-19;

b=10;

forj=1:20

fori=1:13

x1(i,j)=a;

end

a=a+2;

end

x1=reshape(x1,1,260);

fori=1:13

forj=1:20

y1(i,j)=b;

end

b=b+1.92307692307;

end

y1 = flipud(y1);% <------ ADD THIS

y1=reshape(y1,1,260);

%-----------Ensemble mean------------%

addu(1:n)=0;

addv(1:n)=0;

forj=1:n

fori=1:200

addu(1,j)=addu(1,j)+u(i,j);

addv(1,j)=addv(1,j)+v(i,j);

end

avgu(1,j)=addu(1,j)/200;

avgv(1,j)=addv(1,j)/200;

end

%---------------------plot--------%

figure('name','ENSEMBLE AVARAGE FOR 50112')

rectangle('Position',[-26.75 0 53.55 50]);

xlabel('x (mm)')

ylabel('y (mm)')

axisequal

axis ([-25 25 0 50])

set(gca,'ytick',0:10:50)

axistight

holdon

rectangle('Position',[-26.75 40 53.5 5],'FaceColor',[0.502 0.502 0.502])

text (-2.5,42.5,'Piston','Color','black','FontSize',14)

y=yline(1.2,'r','TDC');

y.LineWidth = 1;

y.LabelHorizontalAlignment='center';

y.LabelVerticalAlignment='middle';

holdon

magnitude=sqrt(avgu.^2+avgv.^2);

magnitude=reshape(magnitude,[13,20]);

[X,Y]=meshgrid(-19:2:19,10:1.9230769:33.0769230);

Y = flipud(Y);% <------ ADD THIS

contourf(X,Y,magnitude,100,'linestyle','none')

holdon

[theta, radius] = cart2pol(avgu,avgv);% <------ ADD THIS

[avgu, avgv] = pol2cart(-theta, radius);% <------ ADD THIS

h=quiver(x1,y1,avgu,avgv,'color','black');

% set(gca,'ydir','reverse') % <------ REMOVE THIS

axisequal

c = colorbar;

c.Label.String ='V [m/s]';

c.Label.FontSize =14;

c.FontWeight ='normal';

holdoff

サインインしてコメントする。

その他の回答 (0 件)

サインインしてこの質問に回答する。

カテゴリ

MATLAB Graphics 2-D and 3-D Plots

Find more on2-D and 3-D PlotsinHelp CenterandFile Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by

how to plot two graphs in one plot

5 件のコメント
表示非表示 4 件の古いコメント

採用された回答

6 件のコメント
表示非表示 5 件の古いコメント

その他の回答 (0 件)

参考

カテゴリ

タグ

Community Treasure Hunt

how to plot two graphs in one plot

5 件のコメント 表示非表示 4 件の古いコメント

採用された回答

6 件のコメント 表示非表示 5 件の古いコメント

その他の回答 (0 件)

参考

カテゴリ

タグ

Community Treasure Hunt

5 件のコメント
表示非表示 4 件の古いコメント

6 件のコメント
表示非表示 5 件の古いコメント