קובץ:Ispline order3.svg
תוכן הדף אינו נתמך בשפות אחרות.
מראה
מתוך ויקיפדיה, האנציקלופדיה החופשית
גודל התצוגה המקדימה הזאת מסוג PNG של קובץ ה־SVG הזה: 720 × 540 פיקסלים. רזולוציות אחרות: 320 × 240 פיקסלים | 640 × 480 פיקסלים | 1,024 × 768 פיקסלים | 1,280 × 960 פיקסלים | 2,560 × 1,920 פיקסלים.
לקובץ המקורי (קובץ SVG, הגודל המקורי: 720 × 540 פיקסלים, גודל הקובץ: 43 ק"ב)
זהו קובץ שמקורו במיזם ויקישיתוף. תיאורו בדף תיאור הקובץ המקורי (בעברית) מוצג למטה. |
תקציר
תיאורIspline order3.svg |
English: Plot showing I-splines of order 3. |
תאריך יצירה | |
מקור | נוצר על־ידי מעלה היצירה |
יוצר | Skbkekas |
SVGהתפתחות InfoField | Matplotlib עם נוצרה ה גרפיקה וקטורית |
קוד מקור InfoField | Python codefrom __future__ import division
import numpy as np
import matplotlib.pyplot as plt
def mspline(x, k, i, T):
"""Evaluate the kth order m-spline at x.
Arguments:
----------
x : the point at which the spline is evaluated
k : the order of the spline (k>=1)
i : the interval on which the spline is defined (i<len(T)-k)
T : the knots
Returns:
--------
The value of the spline at x
Note: Use mspline_knots(...) to create T.
"""
if x<T[i]: return 0.
if x>=T[i+k]: return 0.
if k==1: return 1/(T[i+1]-T[i])
d1 = x-T[i]
d2 = T[i+k]-x
v = d1*mspline(x, k-1, i, T) + d2*mspline(x, k-1, i+1, T)
v *= k
v /= (k-1)*(T[i+k]-T[i])
return v
def ispline(x, k, i, T):
"""Evaluate the kth order i-spline at x.
Arguments:
----------
x : the point at which the spline is evaluated
k : the order of the spline
i : the interval on which the spline is defined (i<len(T)-k-1)
T : the knots
Returns:
--------
The value of the spline at x
Note: Use ispline_knots(...) to create T.
"""
j = np.flatnonzero(T<=x)[-1]
if j<i: return 0.
if j-k+1>i: return 1.
v = 0.
for m in range(i,j+1):
v += mspline(x, k+1, m, T)*(T[m+k+1]-T[m])/(k+1)
return v
def mspline_knots(xmin, xmax, k, m):
"""Create an order list of points in the inteval (xmin,xmax)
suitable for use as knots in m-splines and i-splines.
Arguments:
----------
xmin : the lower bound of the interval
xmax : the upper bount of the interval
k : the order of the m-spline for which the knots will be used
m : the number of intervals into which (xmin,xmax) will be partitioned
"""
T = [xmin]*(k-1)
h = (xmax-xmin)/m
T.extend(np.arange(xmin, xmax, h))
T.extend([xmax]*k)
return T
def ispline_knots(xmin, xmax, k, m):
"""Create an order list of points in the inteval (xmin,xmax)
suitable for use as knots in m-splines and i-splines.
Arguments:
----------
xmin : the lower bound of the interval
xmax : the upper bount of the interval
k : the order of the m-spline for which the knots will be used
m : the number of intervals into which (xmin,xmax) will be partitioned
"""
return mspline_knots(xmin, xmax, k+1, m)
def test_ispline():
plt.clf()
plt.axes([0.1,0.1,0.7,0.8])
k = 3
T = ispline_knots(0, 1, k, 5)
A = []
for i in range(1,len(T)-k-1):
F = []
for x in np.arange(0,1.01,0.01):
F.append([x,ispline(x, k, i, T)])
F = np.array(F)
a = plt.plot(F[:,0], F[:,1], '-')
A.append(a)
if i==0: print F
plt.hold(True)
b = plt.figlegend(A, ("i=1", "i=2", "i=3", "i=4", "i=5", "i=6", "i=7"),\
'center right')
b.draw_frame(False)
plt.xlabel(r"$x$")
plt.ylabel(r"$M_i(x
|
רישיון
אני, בעל זכויות היוצרים על עבודה זו, מפרסם בזאת את העבודה תחת הרישיון הבא:
הקובץ הזה מתפרסם לפי תנאי רישיון קריאייטיב קומונז ייחוס-שיתוף זהה 3.0 לא מותאם.
- הנכם רשאים:
- לשתף – להעתיק, להפיץ ולהעביר את העבודה
- לערבב בין עבודות – להתאים את העבודה
- תחת התנאים הבאים:
- ייחוס – יש לתת ייחוס הולם, לתת קישור לרישיון, ולציין אם נעשו שינויים. אפשר לעשות את זה בכל צורה סבירה, אבל לא בשום צורה שמשתמע ממנה שמעניק הרישיון תומך בך או בשימוש שלך.
- שיתוף זהה – אם תיצרו רמיקס, תשנו, או תבנו על החומר, חובה עליכם להפיץ את התרומות שלך לפי תנאי רישיון זהה או תואם למקור.
פריטים שמוצגים בקובץ הזה
מוצג
6 ביוני 2009
היסטוריית הקובץ
ניתן ללחוץ על תאריך/שעה כדי לראות את הקובץ כפי שנראה באותו זמן.
תאריך/שעה | תמונה ממוזערת | ממדים | משתמש | הערה | |
---|---|---|---|---|---|
נוכחית | 01:53, 7 ביוני 2009 | 540 × 720 (43 ק"ב) | Skbkekas | {{Information |Description={{en|1=Plot showing I-splines of order 3.}} |Source=Own work by uploader |Author=Skbkekas |Date=2009-06-06 |Permission= |other_versions= }} <!--{{ImageUpload|full}}--> |
שימוש בקובץ
אין בוויקיפדיה דפים המשתמשים בקובץ זה.
שימוש גלובלי בקובץ
אתרי הוויקי השונים הבאים משתמשים בקובץ זה:
- שימוש באתר en.wikipedia.org