From octave-sources-request at bevo dot che dot wisc dot edu Tue May 25 07:48:11 1999 Subject: Spline interpolation and Savitzki-Golay smoothing From: Zdenek Remes To: octave-sources at bevo dot che dot wisc dot edu Date: Tue, 25 May 1999 14:48:04 +0200 (MET DST) This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime at docserver dot cac dot washington dot edu for more info. ---559023410-959030623-927636484=:8994 Content-Type: TEXT/PLAIN; charset=US-ASCII Hi, I send you too useful function. They work well, however they are not compatible with Matlab. Zdenek Remes ---559023410-959030623-927636484=:8994 Content-Type: TEXT/PLAIN; charset=US-ASCII; name="spline.m" Content-Transfer-Encoding: BASE64 Content-ID: Content-Description: IyMgbmF0dXJhbC1jdWJpYy1zcGxpbmUgaW50ZXJwb2xhdGlvbg0KIyMgdXNh Z2U6IHlzcGxpbmUgPSBzcGxpbmUoeCx5LHhzcGxpbmUpDQojIyBleGFtcGxl OiANCiMjCXggPSAwOjEwOyB5ID0gc2luKHgpOw0KIyMJeHNwbGluZSA9IDA6 MC4xOjEwOyB5c3BsaW5lID0gc3BsaW5lKHgseSx4c3BsaW5lKTsNCiMjCXBs b3QoeCx5LCIrIix4c3BsaW5lLHlzcGxpbmUpOw0KIyMgR2l2ZW4gdGhlIHZl Y3RvcnMgeCBhbmQgeSwgd2hpY2ggdGFidWxhdGUgYSBmdW5jdGlvbiwgd2l0 aCANCiMjIHgoMSkgPCB4KDIpIDwgeCgzKSA8Li4uIG9yIHgoMSkgPiB4KDIp ID4geCgzKSA+Li4uLCBhbmQgZ2l2ZW4NCiMjIHRoZSB2ZWN0b3IgeHNwbGlu ZSwgdGhpcyBmdW5jdGlvbiByZXR1cm5zIGEgbmF0dXJhbC1jdWJpYy1zcGxp bmUgDQojIyBpbnRlcnBvbGF0ZWQgdmVjdG9yIHlzcGxpbmUuIA0KIyMgYXV0 aG9yOiBaZGVuZWsgUmVtZXMsIE1heSAyMiwgMTk5OQ0KDQoNCmZ1bmN0aW9u IHluZXcgPSBzcGxpbmUoeCx5LHhuZXcpDQpbeCxpbmRleF09c29ydCh4KTsN Cnk9eShpbmRleCk7DQpuPWxlbmd0aCh5KTsNCnkyKDEpPTAuMDsNCnkyKG4p PTAuMDsNCnUoMSk9MC4wOw0KZm9yIGk9MjpuLTENCiAgc2lnPSh4KGkpLXgo aS0xKSkvKHgoaSsxKS14KGktMSkpOw0KICBwPXNpZyp5MihpLTEpKzIuMDsN CiAgeTIoaSk9KHNpZy0xLjApL3A7DQogIHUoaSk9KHkoaSsxKS15KGkpKS8o eChpKzEpLXgoaSkpLSh5KGkpLXkoaS0xKSkvKHgoaSkteChpLTEpKTsNCiAg dShpKT0oNi4wKnUoaSkvKHgoaSsxKS14KGktMSkpLXNpZyp1KGktMSkpL3A7 DQplbmRmb3I7DQprPW4tMTsNCndoaWxlIChrID49IDEpDQogIHkyKGspPXky KGspKnkyKGsrMSkrdShrKTsNCiAgay0tOw0KZW5kd2hpbGU7DQoNCmkxPTE7 IGluPWxlbmd0aCh4bmV3KTsNCg0KI2lmICh4bmV3KDEpIDwgeCgxKSkgDQoj ICBlcnJvcigic3BsaW5lOiBiYWQgeHNwbGluZSIpOw0KI2VuZGlmOw0KI2lm ICh4bmV3KGluKSA+IHgobikpIA0KIyAgZXJyb3IoInNwbGluZTogYmFkIHhz cGxpbmUiKTsNCiNlbmRpZjsNCiAgDQppZiAoeG5ldygxKSA9PSB4KDEpKSAN CiAgeW5ldygxKT15KDEpOw0KICBpMT0yOw0KZW5kaWY7DQppZiAoeG5ldyhp bikgPT0geChuKSkgDQogIHluZXcoaW4pPXkobik7DQogIGluPWluLTE7IA0K ZW5kaWY7DQogIA0KDQpmb3IgaT1pMTppbiAgDQogIGtoaT1uOw0KICBrbG89 MTsNCiAgd2hpbGUgKChraGkta2xvKSA+IDEpDQogICAgaz1mbG9vcigoa2hp K2tsbykvMik7DQogICAgaWYgKHgoaykgPiB4bmV3KGkpKSANCiAgICAgIGto aT1rOw0KICAgIGVsc2UNCiAgICAgIGtsbz1rOyAgDQogICAgZW5kaWY7DQog IGVuZHdoaWxlOw0KICBoPXgoa2hpKS14KGtsbyk7DQogIGE9KHgoa2hpKS14 bmV3KGkpKS9oOw0KICBiPSh4bmV3KGkpLXgoa2xvKSkvaDsNCiAgeW5ldyhp KT1hKnkoa2xvKStiKnkoa2hpKSsoKGFeMy1hKSp5MihrbG8pKyhiXjMtYikq eTIoa2hpKSkqKGgqaCkvNi4wOw0KZW5kZm9yOw0KZW5kZnVuY3Rpb247DQo= ---559023410-959030623-927636484=:8994 Content-Type: TEXT/PLAIN; charset=US-ASCII; name="savgol.m" Content-Transfer-Encoding: BASE64 Content-ID: Content-Description: IyMgU2F2aXR6a3ktR29sYXkgc21vb3RoaW5nIGZpbHRlcg0KIyMgdXNhZ2U6 IFt4c2F2Z29sLHlzYXZnb2xdPXNhdmdvbCh4LHksbmwsbnIsbSkNCiMjIGV4 YW1wbGU6IHg9MDowLjAxOjM7eTE9c2luKHguXjMpO3k9eTErKHJhbmQoMSwz MDEpLTAuNSkvMzsNCiMjCSAgICBbeHNhdmdvbCx5c2F2Z29sXT1zYXZnb2wo eCx5LDEwLDEwLDIpOw0KIyMJICAgIHBsb3QoeCx5LCIrIix4c2F2Z29sLHlz YXZnb2wseCx5MSkNCiMjIEdpdmVuIHZlY3RvcnMgeCwgeSBjb250YWluaW5n IGEgdGFidWxhdGVkIGRhdGEgeT1mKHgpIHdpdGggDQojIyBlcXVhbGx5IHNw YWNlZCB4J3MgdGhpcyBmdW5jdGlvbiBjYWxjdWxhdGVzIHNtb290aGVkIGRh dGENCiMjIHlzYXZnb2w9Zyh4c2F2Z29sKSBieSBTYXZpdHpreS1Hb2xheSBz bW9vdGhpbmcgZmlsdGVyLg0KIyMgbmwgaXMgdGhlIG51bWJlciBvZiBsZWZ0 d2FyZCAocGFzdCkgZGF0YSBwb2ludHMgdXNlZCwgd2hpbGUNCiMjIG5yIGlz IHRoZSBudW1iZXIgb2YgcmlnaHR3YXJkIChmdXR1cmUpIGRhdGEgcG9pbnRz LCBtYWtpbmcNCiMjIHRoZSB0b3RhbCBudW1iZXIgb2YgZGF0YSBwb2ludHMg dXNlZCBubCtucisxLiBtIGlzIHRoZSBvcmRlcg0KIyMgb2YgdGhlIHNtb290 aGluZyBwb2x5bm9taWFsLCBhbHNvIGVxdWFsIHRvIHRoZSBoaWdoZXN0DQoj IyBjb25zZXJ2ZWQgbW9tZW50OyB1c3VhbCB2YWx1ZXMgYXJlIG09MiBvciBt PTQuDQojIyBUaGUgaWRlYSBvZiBTYXZpdHpreS1Hb2xheSBmaWx0ZXJpbmcg aXMgdG8gc21vb3RoIHRoZQ0KIyMgdW5kZXJseWluZyBkYXRhIHk9Zih4KSB3 aXRoaW4gdGhlIG1vdmluZyB3aW5kb3cgbm90IGJ5IGENCiMjIGNvbnN0YW50 ICh3aG9zZSBlc3RpbWF0ZSBpcyB0aGUgYXZlcmFnZSksIGJ1dCBieSBhIHBv bHktDQojIyBub21pYWwgb2YgaGlnaGVyIG9yZGVyLiBUaHVzIGZvciBhIHBv aW50IHkoaSkgdGhlIGZ1bmN0aW9uDQojIyBzYXZnb2wgZml0cyBieSBhIGxl YXN0LXNxdWFyZXMgbWV0aG9kIGEgcG9seW5vbWlhbCB0bw0KIyMgcG9pbnRz IHkoaS1ubCksIC4uLiwgeShpK25yKSBpbiB0aGUgbW92aW5nIHdpbmRvdywg YW5kDQojIyB0aGVuIHNldCBnKGktbmwrMSkgdG8gdGhlIHZhbHVlIG9mIHRo YXQgcG9seW5vbWlhbCBhdA0KIyMgcG9zaXRpb24geChpKS4gDQojIyBaZGVu ZWsgUmVtZXMsIE1haSAyMiwgMTk5OQ0KCQ0KZnVuY3Rpb24gW3huZXcseW5l d109c2F2Z29sKHgseSxubCwgbnIsIE0pDQ0KICAgIGlmIG1heChkaWZmKHgs MikpPjEwMCplcHMgDQoJZXJyb3IoIlRoZSB4J3MgbXVzdCBiZSBlcXVhbGx5 IHNwYWNlZC4iKQ0KICAgIGVuZGlmDQogICAgZm9yIGk9LW5sOm5yDQoJZm9y IGo9MDpNDQoJICAgIEEoaStubCsxLGorMSk9aV5qOw0KCWVuZGZvcg0KICAg IGVuZGZvcg0KICAgIEFBPWludihBJypBKTsNCiAgICBmb3IgaT0tbmw6bnIN CgljYz0wOw0KCWZvciBtPTA6TQ0KCSAgICBjYz1jYytBQSgxLG0rMSkqaV5t Ow0KCWVuZGZvcg0KICAgIGMoaStubCsxKT1jYzsNCiAgICBlbmRmb3INCiAg ICANCiAgICBueD1sZW5ndGgoeCk7DQogICAgZm9yIGk9bmw6bngtbnItMQ0K CXl5PTA7DQoJZm9yIGo9LW5sOm5yDQoJICAgIHl5PXl5K2MoaitubCsxKSp5 KGkraisxKTsNCgllbmRmb3INCiAgICAgICAgeG5ldyhpLW5sKzEpPXgoaSsx KTsNCgl5bmV3KGktbmwrMSk9eXk7DQogICAgZW5kZm9yICAgIA0KZW5kZnVu Y3Rpb24= ---559023410-959030623-927636484=:8994--