你正在与曲线,而不是功能。
你的
翼型
是(平面)曲线,或样条函数工具箱措辞2-valued, 1-variate样条函数。这个函数是由两个1-valued 1-variate函数插值对x (t(我),(我),(t (i), y(我))。与
fnval(机翼,s1)
你要求的(x, y)坐标的值
t = 0.1
。
你可以看到这个情节两个组件
airfoil_x = fncmb(机翼(1,0));
airfoil_y = fncmb(机翼,[0,1]);
图(3)
clf
h1 =次要情节(2,1,1);
fnplt (airfoil_y)
持有在
情节(t y“o”)
包含(“t”),ylabel (“y”)
h2 =次要情节(2,1,2);
fnplt (airfoil_x)
持有在
情节(t x,“o”)
包含(“t”),ylabel (“x”)
你要的是发现在这些情节
errorbar (h1, s1, s2 (2), s2 (2) - h1.YLim (1), 0,“k *’)
errorbar (h2, s1, s2 (1), s2 (1) - h2.YLim (1), 0,“k *’)
和你原来的情节
图(2)
情节(s2 (1), (2),“gsq”,“MarkerFaceColor”,‘g’)
如果你想点(s)在机翼在x = 0.1你可以这样做:
- 找到0 (s)坐标的函数-所需的x值
- 评估你的函数的值t
x0 = 0.1;
t0 =意味着(fnzeros (fncmb (airfoil_x,“- - -”x0) 1);
y0 = fnval (airfoil_y、t0);
airfoil_0 = fnval(机翼、t0);
情节(h2, [0, 1], [x0, x0],“g -”)
errorbar (h2, t0, x0, x0, 5个,5个,5个,5个,“gsq”,“MarkerFaceColor”,‘g’)
h2。YLim = [0, 1];
h1。YLim = h1.YLim;
errorbar (h1、t0, y0, 5个,5个,(0,0),“gsq”,“MarkerFaceColor”,‘g’)