ช่วงนี้กำลังกลุ้มใจมาก เพราะต้องหา Jacobian ของ Spline ปัญหาคือมันไม่สามารถหาเชิงวิเคราะห์ได้ จำเป็นต้องหาเชิงเลขบน MATLAB เพราะต้องการหาจุดตรึงของ Spline ที่เหมาะสมที่สุด
ทางเลือกหนึ่งคือใช้ฟังก์ชัน fndir ซึ่งมากับ Spline Toolbox อันนี้ไม่ค่อยอยากใช้นักเพราะว่าเครื่องที่บ้านไม่มี Toolbox อันนี้ (เวอร์ชันนักเรียน) ต้องไปใช้ที่ห้องทำงาน อีกทางเลือกหนึ่งดูเป็นเรื่องทั่วไปกว่า โดยการใช้ Automatic Numerical Differentiation ของคุณ John D’Errico ซึ่งยังไม่รู้ว่าดีแค่ไหน
ตัวผมเองไม่ได้สนใจศึกษาเรื่อง Numerical Computing เลย จนหลัง ๆ อ่านมากชักจะชอบไปเจอบทเรียน Numerical Methods in Mechanics ของ Lund University ก็ขอจดไว้หน่อย รหัสต้นข้างล่างนี้ลอกมาแปะไว้เลย ตัวฟังก์ชันชื่อ jacobian ก็ควรจะเปลี่ยนนะครับเพราะไปซ้ำกับ jacobian แบบ symbolic ของ MATLAB เค้า
function [J]=jacobian(func,x)% [...]
Filed under: MATLAB, Numerical Computing | Tagged: MATLAB, Numerical Computing | Leave a Comment »
