发布绘图流线(按升序输入)

1次观看(过去30天)

显示旧的注释

艾迪生柯林斯 2021年8月13日

编辑: 艾迪生柯林斯2021年8月13日

答:接受艾迪生柯林斯

final_data_run1.mat

你好,

我在绘制流线时遇到了问题。当我试图绘制我的流线时，我得到一个错误，因为它使用interp1和我的一个输入流线不是“升序”。我已经附上了我的代码和数据文件。流线应该像箭袋图中的矢量一样从左到右。

箭袋图是这样的:

错误信息如下:

我需要对streamlined()的输入做什么才能使其正常工作?

                          clc;明确;接近所有；
                         
                          %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                         
                          %设置情节的东西:
                         
                          集(0,“DefaultLineLineWidth”, 1.5)
                         
                          %集(0“DefaultLineColor”(1,1,1))
                         
                          集(0,“DefaultLineMarkerSize”15)
                         
                          集(0,“DefaultAxesFontSize”, 20)
                         
                          集(0,“DefaultFigureColor”, 1, 1, 1)
                         
                          集(0,“DefaultTextFontSize”, 18)
                         
                          %设置(0,“DefaultTextInterpreter”,“乳胶”)
                         
                          集(0,“DefaultTextFontName”，“times roman那样”）
                         
                          集(0,“DefaultAxesFontName”，“times roman那样”）
                         
                          %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                         
                          runname =“M.25_0Y_A”；
                         
                          加载数据
                         
                          负载([“final_data_run1.mat”])
                         
                          Uavg(Uavg==0) = NaN;将空数据设置为NaN
                         
                          Vavg(Vavg==0) = NaN;将空数据设置为NaN
                         
                          Yi = sort(Y(1，:))';
                         
                          创建xq和yq的组合
                         
                          如果文件中的x和y数据不在网格中，则插入数据
                         
                          Uq = rot90(Uavg);
                         
                          Vq = rot90(Vavg);
                         
                          %%探针壁(需要最终旋转这些)
                         
                          Ytip = 4.18-1;%左探测壁的左下顶点
                         
                          Xtip = -7.703;%左探测壁的左下顶点
                         
                          宽度= 12.7;%毫米
                         
                          高度= 1;%毫米
                         
                          差距= 3;% mm -探针壁之间的距离(壁的左端)
                         
                          %%新的坐标(原点将在探针尖端)
                         
                          Xcenter = xtip;探针尖端的x中心点
                         
                          Ycenter = ytip-height;探针尖端y型中心点
                         
                          Xtip = Xtip -xcenter;%左探测壁的左下顶点
                         
                          Ytip = Ytip -ycenter;%左探测壁的左下顶点
                         
                          Xq = Xq+ycenter;
                         
                          Yq = Yq-xcenter;
                         
                          %%绘制向量、探测壁和流线
                         
                          图=图(渲染器的，“画家”，“位置”， [1300 10 947 900])
                         
                          箭袋(Yq -Xq、Vq uq, 1);%， 'LineWidth'，2) %绘制向量
                         
                          持有在；
                         
                          yline(0);参照线(0);持有在绘图探测中心线
                         
                          L =标题({“test_line”}，“翻译”，“没有”）;
                         
                          标题([平均速度和流线(runname,“)”])
                         
                          % %绘制探测壁
                         
                          矩形(“位置”, xtip ytip,宽度,高度,“FaceColor”，“k”）
                         
                          矩形(“位置”, xtip ytip-gap,宽度,高度,“FaceColor”，“k”）
                         
                          包含(“x(毫米)）
                         
                          ylabel (“y”(毫米)）
                         
                          S = 0.10;
                         
                          启动= -2.5:s:2.5;
                         
                          startx = (min(min(Yq))+.25)*ones(size(starty));
                         
                          %%图形流线
                         
                          % hgridinterp = streamlined (Yq,Xq,Vq,Uq,startx,starty);抓住;
                         
                          hgridinterp =流线(Yq，-Xq,Vq，-Uq,startx,starty);持有在；
                         
                          集(hgridinterp,“颜色”，“红色”，“线宽”3);