我正在尝试使用 GNU 科学库 (GSL) 绘制从 A 到 B 的平滑路径。我使用的 API 返回少量(在本例中为 8 个)不规则间隔的点(红色) )，您可以在下图中看到:
紫色点代表我希望看到从 GSL 返回的点。
首先，这种2D B-Spline形状可以用GSL得到吗？我不太了解 B 样条曲线，更不用说 2D B 样条曲线了。我能够获得显示的 B 样条曲线示例 here运行并创建一个流畅的 .ps 文件没有问题，但该示例使用带有以下代码的统一断点:
/* use uniform breakpoints on [0, 15] */
gsl_bspline_knots_uniform(0.0, 15.0, bw);
在我的例子中，鉴于我提供的数据不稳定且间隔不均匀，我是否必须使用非均匀结？我尝试使用 gsl_bspline_knots()，以便在以下测试代码中使用非统一断点，但我真的不确定这是否是正确的方向。
#define NCOEFFS 8 // not sure what this number should be - number of data points?
#define NBREAK (NCOEFFS - 2)
const size_t nbreak = NBREAK;
int main (void) {
// (example code)...
gsl_vector *non_uniform = gsl_vector_alloc(nbreak);
// create some random breakpoint values
for (i=0; i<nbreak; i++) {
double val = gsl_ran_gaussian(r, 2.0);
printf("val: %f\n", val);
gsl_vector_set(non_uniform, i, val);
}
gsl_bspline_knots(non_uniform, bw);
// (more example code)...
}
此外，我将如何翻译上面在 2D x/y 坐标空间中绘制 B 样条曲线的示例？如果 GNU 科学库不适合这个，有人可以推荐更合适的 C/C++ 库吗？
任何方向上的帮助或指示将不胜感激。
最佳答案
首先:一维基样条
给定一组 NBREAK断点 (t_1, ..., t_{NBREAK}) , 有 NCOEFFS=NBREAK+2三次 b 样条组件 B_j(t) .这些函数及其一阶和二阶导数始终是连续的，即使在断点处也是如此。因此，线性组合 f(t) = \sum m_j B_j(t) 给出的任何拟合也将共享这些属性(类似于自然三次样条)。 b样条分量的数量，NCOEFFS ，不需要等于数据点的数量，NDATA .如果NCOEFFS < NDATA ，您应该使用最小二乘最小化来获得拟合(GSL 文档有一个很好的最小二乘计算示例来获得 b 样条拟合 here )。 NCOEFFS < NDATA当数据包含噪声时是不错的选择，而您的情况似乎并非如此
系数个数不等于断点个数的原因，NCOEFFS=NBREAK+2 , 与处理基样条时没有指定边界条件有关。鉴于人们通常更熟悉自然三次样条，因此值得评论的是 natural cubic splines施加边界条件 d^2f(x)/dx^2=0 .这就是为什么使用三次多项式基础的自然三次样条的任何表示都会有 NCOEFFS=NBREAK . Here is a link很好地解释了由三次多项式的系数给出的自由度的计数，这些系数表示自然 b 样条和施加连续性所需的方程数 f(t) , df(t)/dt和 d^2f(t)/dt^2 ).
最后:使用 b 样条拟合参数曲线。
你有一组“数据”点(x_1, y_1)....(x_{NDATA},y_{NDATA})并且您想构建参数拟合 P(t)=( f_1(t), f_2(t) ) .如果 NCOEFFS<NDATA，则 B 样条拟合不会遍历所有数据点。 (如果您仔细选择断点和 NCOEFFS=N_DATA ，您可以要求这样做)。在我的研究中，我只使用一维非参数拟合 (y=f(x))，但我相信这个参数案例是相似的。我会尝试以下操作
第 1 步:创建“数据”点集 (t, x) = {(1, x_1), (2, x_2)...(NDATA, x_{NDATA})}并使用 gsl 1D b 样条曲线来拟合它们。这种合身会给你f_1(t) = sum_{i=1}^{NCOEFFS} mx_j B_j(t)与 t \in [1,NDATA] .
第 2 步:现在构建“数据”点集 (t, y) = {(1, y_1), (2, y_2)...(NDATA, y_{NDATA})}并使用 b 样条曲线来拟合它们。这将为您提供f_2(t) = sum_{i=1}^{NCOEFFS} my_j B_j(t)与 t \in [1,NDATA]
现在绘制P(t)=( f_1(t), f_2(t) ), t \in [1,NDATA] .基本上，我在 2 个一维非参数拟合(这是 GSL 提供的)中映射了一个二维参数曲线问题。
最后一点是第一步和第二步的断点选择(以及基元数量 NCOEFFS )，只要覆盖范围 t\in[1, NDATA] 即可。和 NCOEFFS <= NDATA ，断点的选择是任意的。我相信，如果您选择断点为 {1, 3, ..., NDATA-2, NDATA }拟合将通过数据点(注意我跳过了内部点 t=2 和 t=NDATA-1 这样 NBREAK=NDATA-2 和 NCOEFFS=NDATA )。就是这样NAG库选择断点以获得插值拟合(意思是:通过数据点的拟合)。
关于c++ - 使用 GNU 科学库 (GSL) 使用不均匀分布的点绘制 2D B 样条路径，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/24194909/
我正在学习如何使用Nokogiri，根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div
我有一个Ruby程序，它使用rubyzip压缩XML文件的目录树。gem。我的问题是文件开始变得很重，我想提高压缩级别，因为压缩时间不是问题。我在rubyzipdocumentation中找不到一种为创建的ZIP文件指定压缩级别的方法。有人知道如何更改此设置吗？是否有另一个允许指定压缩级别的Ruby库？ 最佳答案 这是我通过查看rubyzip内部创建的代码。level=Zlib::BEST_COMPRESSIONZip::ZipOutputStream.open(zip_file)do|zip|Dir.glob("**/*")d
类classAprivatedeffooputs:fooendpublicdefbarputs:barendprivatedefzimputs:zimendprotecteddefdibputs:dibendendA的实例a=A.new测试a.foorescueputs:faila.barrescueputs:faila.zimrescueputs:faila.dibrescueputs:faila.gazrescueputs:fail测试输出failbarfailfailfail.发送测试[:foo,:bar,:zim,:dib,:gaz].each{|m|a.send(m)resc
很好奇，就使用rubyonrails自动化单元测试而言，你们正在做什么？您是否创建了一个脚本来在cron中运行rake作业并将结果邮寄给您？git中的预提交Hook？只是手动调用？我完全理解测试，但想知道在错误发生之前捕获错误的最佳实践是什么。让我们理所当然地认为测试本身是完美无缺的，并且可以正常工作。下一步是什么以确保他们在正确的时间将可能有害的结果传达给您？ 最佳答案 不确定您到底想听什么，但是有几个级别的自动代码库控制:在处理某项功能时，您可以使用类似autotest的内容获得关于哪些有效，哪些无效的即时反馈。要确保您的提
假设我做了一个模块如下:m=Module.newdoclassCendend三个问题:除了对m的引用之外，还有什么方法可以访问C和m中的其他内容？我可以在创建匿名模块后为其命名吗(就像我输入“module...”一样)？如何在使用完匿名模块后将其删除，使其定义的常量不再存在？ 最佳答案 三个答案:是的，使用ObjectSpace.此代码使c引用你的类(class)C不引用m:c=nilObjectSpace.each_object{|obj|c=objif(Class===objandobj.name=~/::C$/)}当然这取决于
我正在尝试使用ruby和Savon来使用网络服务。测试服务为http://www.webservicex.net/WS/WSDetails.aspx?WSID=9&CATID=2require'rubygems'require'savon'client=Savon::Client.new"http://www.webservicex.net/stockquote.asmx?WSDL"client.get_quotedo|soap|soap.body={:symbol=>"AAPL"}end返回SOAP异常。检查soap信封，在我看来soap请求没有正确的命名空间。任何人都可以建议我
关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点？我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等)，但我确实想创建一个输出文件。
我在我的项目目录中完成了compasscreate.和compassinitrails。几个问题:我已将我的.sass文件放在public/stylesheets中。这是放置它们的正确位置吗？当我运行compasswatch时，它不会自动编译这些.sass文件。我必须手动指定文件:compasswatchpublic/stylesheets/myfile.sass等。如何让它自动运行？文件ie.css、print.css和screen.css已放在stylesheets/compiled。如何在编译后不让它们重新出现的情况下删除它们？我自己编译的.sass文件编译成compiled/t
我想将html转换为纯文本。不过，我不想只删除标签，我想智能地保留尽可能多的格式。为插入换行符标签，检测段落并格式化它们等。输入非常简单，通常是格式良好的html(不是整个文档，只是一堆内容，通常没有anchor或图像)。我可以将几个正则表达式放在一起，让我达到80%，但我认为可能有一些现有的解决方案更智能。 最佳答案 首先，不要尝试为此使用正则表达式。很有可能你会想出一个脆弱/脆弱的解决方案，它会随着HTML的变化而崩溃，或者很难管理和维护。您可以使用Nokogiri快速解析HTML并提取文本:require'nokogiri'h
在railstutorial中，作者为什么选择使用这个(代码list10.25):http://ruby.railstutorial.org/chapters/updating-showing-and-deleting-usersnamespace:dbdodesc"Filldatabasewithsampledata"task:populate=>:environmentdoRake::Task['db:reset'].invokeUser.create!(:name=>"ExampleUser",:email=>"example@railstutorial.org",:passwo