JP2016034483A - Golf club fitting device, method and program - Google Patents

Abstract

The present invention provides a fitting device that objectively and accurately determines an optimal swingability index, which is a golf club swingability index suitable for a golfer. A fitting device is a fitting device that determines an optimal swingability index that is a golf club swingability index suitable for a golfer, and includes an acquisition unit, a calculation unit, and a determination unit. The acquisition unit acquires a measurement value obtained by measuring a swing motion of a test club by the golfer using a measuring device. The calculation unit calculates a swing index indicating a feature amount of the swing motion based on the measurement value. The determination unit determines the optimal swing ease index according to the size of the swing index. [Selection] Figure 17

Description

  The present invention relates to a fitting device, a method, and a program for determining a golf club swingability index (hereinafter referred to as an optimal swingability index) suitable for a golfer.

  2. Description of the Related Art Conventionally, various fitting methods have been proposed in which a golf club is made a test club and the movement thereof is measured by a measuring device, and a golf club suitable for the golfer is determined based on the measured value. Patent Document 1 points out that it is important to consider the weight of the golf club, the length of the golf club, and the like in order to select a golf club suitable for the golfer.

JP 2013-226375 A

  The weight of the golf club is one index representing the ease of swinging the golf club. By the way, at the time of fitting, the weight of the golf club suitable for the golfer (hereinafter referred to as the optimum club weight) may be determined. The optimum club weight is determined by the weight of the golf club currently used by the golfer, In many cases, it is empirically determined by a skilled club fitter based on the result of a golf club trial hit, an apparent swing tempo, and the like. However, such a fitting method that depends on experience and intuition is not objective, and there is a problem that individual differences occur in the fitting results. In addition, the present inventors considered that it is important to pay attention to a swing ease index other than the weight of the golf club in order to determine a golf club suitable for a golfer by fitting.

  An object of the present invention is to provide a fitting device, a method, and a program that objectively and accurately determine an optimal swingability index that is a golf club swingability index suitable for a golfer.

  A fitting device according to a first aspect of the present invention is a fitting device that determines an optimal swingability index that is a golf club swingability index suitable for a golfer, and includes an acquisition unit, a calculation unit, and a determination unit. Is provided. The acquisition unit acquires a measurement value obtained by measuring a swing motion of a test club by the golfer using a measuring device. The calculation unit calculates a swing index indicating a feature amount of the swing motion based on the measurement value. The determination unit determines the optimal swing ease index according to the size of the swing index.

  The fitting apparatus which concerns on the 2nd viewpoint of this invention is a fitting apparatus which concerns on a 1st viewpoint, Comprising: The said calculation part calculates the said kind of said swing parameter | index. The determining unit determines the optimum swing ease index according to the sizes of the plurality of types of swing indices.

  A fitting device according to a third aspect of the present invention is the fitting device according to the first or second aspect, wherein the swing index is an index representing energy or torque exerted by the golfer during the swing operation, or Indices correlated with these are included.

  A fitting device according to a fourth aspect of the present invention is the fitting device according to the third aspect, wherein the swing index includes arm energy of the golfer during the swing operation, and a shoulder circumference of the golfer during the swing operation. And at least one of torque and head speed.

  A fitting device according to a fifth aspect of the present invention is the fitting device according to the fourth aspect, wherein the calculation unit calculates a head speed during the swing operation. The determining unit determines the optimal swingability index according to the magnitude of the head speed in addition to the magnitude of at least one of the arm energy and the torque around the shoulder.

  A fitting device according to a sixth aspect of the present invention is the fitting device according to any one of the first to fifth aspects, wherein the determining unit determines that the optimal swingability index increases as the swing index increases or decreases. Is determined to be a large value.

A fitting device according to a seventh aspect of the present invention is the fitting device according to any one of the first to sixth aspects, wherein the size of the swing index and the optimum swing ease index for each type of the test club A storage unit is further provided for storing correspondence data for determining a correspondence relationship with the size of the. The determination unit determines the optimal swingability index by referring to the correspondence relationship data in the storage unit according to the type of the test club.

  A fitting device according to an eighth aspect of the present invention is the fitting device according to any one of the first to seventh aspects, wherein the ease of swinging index includes the weight of the golf club, the moment of inertia of the golf club. And at least one moment of inertia about the golfer's shoulder.

A fitting method according to a ninth aspect of the present invention is a fitting method for determining an optimal swingability index that is a golf club swingability index suitable for a golfer, and includes the following steps.
(1) A step of measuring a swing motion of the test club by the golfer with a measuring device.
(2) A step of calculating a swing index indicating a characteristic amount of the swing motion based on a measurement value obtained by the measuring device.
(3) A step of determining the optimum swing ease index according to the size of the swing index.

A fitting program according to a tenth aspect of the present invention is a fitting program for determining an optimum swinging ease index that is a golf club swinging ease index suitable for a golfer, and causes a computer to execute the following steps.
(1) A step of acquiring a measured value obtained by measuring a swing motion of the test club by the golfer with a measuring device.
(2) A step of calculating a swing index indicating a characteristic amount of the swing motion based on a measurement value obtained by the measuring device.
(3) A step of determining the optimum swing ease index according to the size of the swing index.

  There may be a certain relationship between the optimal swing ease index and the characteristic amount of the predetermined swing motion. For example, the present inventors have determined the optimum club weight, the optimum value of the moment of inertia of the golf club (hereinafter referred to as optimum club MI), the optimum value of the moment of inertia around the golfer's shoulder (hereinafter referred to as optimum swing MI), It has been discovered that there is a correlation between the energy or torque exerted by the golfer. Therefore, according to the present invention, the optimal swing ease index is determined according to the characteristic amount of the swing motion and the size of an index (swing index) having a certain relationship with the optimal swing ease index. Thereby, the optimal swing ease index can be objectively determined with high accuracy. The optimum club MI includes an optimum value of the moment of inertia around the grip end of the golf club (hereinafter, optimum grip end MI) and an optimum value of the moment of inertia around the center of gravity of the golf club.

  In particular, in the invention according to both the third aspect and the eighth aspect, depending on the size of an index representing energy or torque exerted by the golfer during the swing operation of the test club by the golfer, or an index correlated therewith. At least one of the optimum club weight, the optimum club MI, and the optimum swing MI is determined. Thereby, at least one of the optimum club weight, the optimum club MI, and the optimum swing MI can be objectively determined with high accuracy.

The figure which shows a golf swing system provided with the fitting apparatus which concerns on 1st Embodiment of this invention. The functional block diagram of a fitting system. The figure explaining xyz local coordinate system on the basis of the grip of a golf club. (A) The figure which shows an address state. (B) The figure which shows a top state. (C) The figure which shows an impact state. (D) The figure which shows a finish state. The flowchart which shows the flow of 1st conversion processing. The flowchart which shows the flow of the process which derives | leads-out the time of an address. The figure explaining a swing plane. The flowchart which shows the flow of a shoulder behavior derivation | leading-out process. The figure which illustrates a double pendulum model notionally. The flowchart which shows the flow of an parameter | index calculation process. Another diagram conceptually explaining the double pendulum model. The figure explaining the arm energy which is one of the swing parameter | indexes. The flowchart which shows the flow of the optimal total weight determination process. The figure which plotted the value of the arm energy and average shoulder torque at the time of swing operation | movement by 10 professional model users and 10 average model users. The figure which shows a pro model area | region and an average model area | region. The figure which shows the relationship between the weight of a golf club, and head speed. The figure which shows the head speed and optimal club weight at the time of the swing operation | movement by the professional model user calculated | required by experiment. The figure which shows the head speed and the optimal club weight at the time of swing operation | movement by the average model user calculated | required by experiment. The figure which shows the division area for allocating the optimal club weight to a professional model user. The figure which shows the division area for allocating the optimal club weight to an average model user. The functional block diagram of the fitting system which concerns on 2nd Embodiment of this invention. The flowchart which shows the flow of the optimal swing MI determination process. The figure which shows the head speed at the time of swing operation | movement by several golf players, the optimal swing MI, and the optimal grip end MI calculated | required by experiment. The figure which shows the relationship between a swing inertia moment and a grip end inertia moment, and head speed. The figure which shows the division area for assigning the optimal swing MI. The functional block diagram of the fitting system which concerns on 3rd Embodiment of this invention. The flowchart which shows the flow of the optimal grip end MI determination process. The figure which shows the division area for assigning the optimal grip end MI. The figure explaining the various angles which are the swing parameter | index which concerns on a modification.

  Hereinafter, a golf club fitting device, method, and program according to some embodiments of the present invention will be described with reference to the drawings.

<1. First Embodiment>
<1-1. Schematic configuration of fitting system>
1 and 2 show an overall configuration of a fitting system 100 including a fitting device 2 according to the present embodiment. The fitting device 2 is a device that analyzes the swing motion of the golf club 4 based on measurement data obtained by measuring the swing motion of the golf club 4 by the golfer 7. In the present embodiment, the fitting device 2 is used for the purpose of supporting the fitting of the golf club 4. The measurement of the swing motion is performed by the sensor unit 1 (measuring device) attached to the grip 42 of the golf club 4, and the fitting device 2 constitutes the fitting system 100 together with the sensor unit 1.

  Hereinafter, after describing the configuration of the sensor unit 1 and the fitting device 2, the flow of the analysis process of the swing operation will be described.

<1-1-1. Configuration of sensor unit>
As shown in FIGS. 1 and 3, the sensor unit 1 is attached to an end portion of the grip 42 of the golf club 4 opposite to the head 41 and measures the behavior of the grip 42. The golf club 4 is a general golf club and includes a shaft 40, a head 41 provided at one end of the shaft 40, and a grip 42 provided at the other end of the shaft 40. The sensor unit 1 is configured to be small and light so as not to hinder the swing operation. As shown in FIG. 2, an acceleration sensor 11, an angular velocity sensor 12, and a geomagnetic sensor 13 are mounted on the sensor unit 1 according to the present embodiment. The sensor unit 1 is also equipped with a communication device 10 for transmitting measurement data from these sensors 11 to 13 to the external fitting device 2. In the present embodiment, the communication device 10 is wireless so as not to hinder the swing operation, but may be connected to the fitting device 2 in a wired manner via a cable.

The acceleration sensor 11, the angular velocity sensor 12, and the geomagnetic sensor 13 respectively measure grip acceleration, grip angular velocity, and grip geomagnetism in the xyz station coordinate system with the grip 42 as a reference. More specifically, the acceleration sensor 11 measures grip accelerations a _x , a _y , and a _z in the x-axis, y-axis, and z-axis directions. The angular velocity sensor 12 measures grip angular velocities ω _x , ω _y , and ω _z around the x axis, the y axis, and the z axis. Geomagnetic sensor 13 measures the x-axis, y-axis and z-axis direction of the grip geomagnetism m _x, m _y, a m _z. These measurement data are acquired as time-series data of a predetermined sampling period Δt. The xyz local coordinate system is a three-axis orthogonal coordinate system defined as shown in FIG. That is, the z axis coincides with the direction in which the shaft 40 extends, and the direction from the head 41 toward the grip 42 is the z axis positive direction. The x-axis is oriented as much as possible in the toe-heel direction of the head 41, and the y-axis is oriented as much as possible in the normal direction of the face surface of the head 41.

  In the present embodiment, measurement data obtained by the acceleration sensor 11, the angular velocity sensor 12, and the geomagnetic sensor 13 are transmitted to the fitting device 2 in real time via the communication device 10. However, for example, the measurement data may be stored in a storage device in the sensor unit 1, and the measurement data may be taken out from the storage device after the swing operation is finished and transferred to the fitting device 2.

<1-1-2. Configuration of fitting device>
The configuration of the fitting device 2 will be described with reference to FIG. The fitting device 2 is manufactured by installing the fitting program 3 according to the present embodiment stored in a computer-readable recording medium 20 such as a CD-ROM or USB memory from the recording medium 20 to a general-purpose personal computer. Is done. The fitting program 3 is software for analyzing the swing motion based on the measurement data sent from the sensor unit 1 and outputting information that assists in selecting a golf club suitable for the golfer 7. The fitting program 3 causes the fitting apparatus 2 to execute an operation described later.

  The fitting device 2 includes a display unit 21, an input unit 22, a storage unit 23, a control unit 24, and a communication unit 25. These units 21 to 25 are connected via the bus line 26 and can communicate with each other. In the present embodiment, the display unit 21 is configured with a liquid crystal display or the like, and displays information to be described later to the user. In addition, a user here is a general term for the person who requires the result of fitting, such as golfer 7 himself or its instructor. The input unit 22 can be configured with a mouse, a keyboard, a touch panel, and the like, and accepts an operation from the user on the fitting device 2.

  The storage unit 23 is configured by a nonvolatile storage device such as a hard disk. The storage unit 23 stores the fitting program 3 and the measurement data sent from the sensor unit 1. The storage unit 23 stores correspondence data 28. The correspondence relationship data 28, which will be described in detail later, is defined for each golf club 4 of various models and is data indicating conditions for determining the optimum club weight. The communication unit 25 is a communication interface that enables communication between the fitting device 2 and an external device, and receives data from the sensor unit 1.

  The control unit 24 can be composed of a CPU, ROM, RAM, and the like. The control unit 24 reads and executes the fitting program 3 in the storage unit 23 to virtually acquire the acquisition unit 24A, the grip behavior deriving unit 24B, the shoulder behavior deriving unit 24C, the index calculating unit 24D, the determining unit 24E, and the display It operates as the control unit 24F. Details of the operations of the units 24A to 24F will be described later.

<1-2. Analysis of swing motion>
Next, a swing operation analysis process for fitting the golf club 4 by the fitting system 100 will be described. The analysis processing according to the present embodiment includes the following seven steps.
(1) Grip acceleration a _x in the xyz local coordinate system, a _y, a _z, grip angular velocity ω _x, ω _y, ω _z and grip geomagnetism m _x, m _y, measuring step of measuring the measurement data m _z ( 2) First, the measurement data in the xyz local coordinate system obtained in the measurement process is converted into grip accelerations a _X , a _Y , a _Z and grip angular velocities ω _X , ω _Y , ω _Z in the XYZ global coordinate system. Conversion step (In the first conversion step, grip speeds v _X , v _Y , and v _Z in the XYZ global coordinate system are also derived.)
(3) The behavior of the grip 42 (grip angular velocities ω _X , ω _Y , ω _Z and grip velocities v _X , v _Y , v _Z ) in the XYZ global coordinate system is determined from the grip 42 in the swing plane P (described later). (4) Shoulder behavior deriving step for deriving a pseudo shoulder behavior of the golfer 7 in the swing plane P based on the behavior of the grip 42 in the swing plane P (4) 5) Based on the behavior of the grip 42 and the pseudo shoulder behavior, a swing index (in the present embodiment, the golfer 7) serving as an index for determining the optimal swingability index (in the present embodiment, the optimal club weight). (6) Index calculation step for calculating average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h (6). Optimal total weight determination for determining the optimal club weight based on the swing index Determination step (7) Optimal shaft determination step for calculating the weight (hereinafter referred to as optimal shaft weight) and bending rigidity (hereinafter referred to as optimal bending stiffness) of the shaft 40 suitable for the golfer 7 Hereinafter, these steps will be described in order.

  The XYZ global coordinate system is a three-axis orthogonal coordinate system defined as shown in FIG. That is, the Z-axis is a direction from the vertically lower side to the upper side, the X-axis is a direction from the back to the stomach of the golfer 7, and the Y-axis is parallel to the ground plane and goes from the ball hitting point to the target point. Direction.

<1-2-1. Measurement process>
In the measurement process, the golf club 4 with the sensor unit 1 is swung by the golfer 7. The golf club 4 swung in the measurement process is one of the two test clubs. However, as will be described in detail later, when it is determined that the first selected test club is not suitable for the golfer 7, another golf club is also swung. These test clubs are different types of golf clubs. In the present embodiment, one test club is a professional golf club (hereinafter referred to as a professional model club), and the other is a golf club suitable for average users (hereinafter referred to as a golf club). Average model club). In the present embodiment, the professional model club is heavier than the average model club. Which test club is to be swung in the measurement process is determined based on the preference and experience of the golfer 7.

Subsequently, as described above golf club 4 grip acceleration a _x during the swing operation, a _y, a _z, grip angular velocity ω _x, ω _y, ω _z and grip geomagnetism m _x, m _y, m _z of the measurement data Is measured by the sensor unit 1. This measurement data is transmitted to the fitting device 2 via the communication device 10 of the sensor unit 1. On the other hand, on the fitting device 2 side, the acquisition unit 24 </ b> A receives this via the communication unit 25 and stores it in the storage unit 23. In this embodiment, at least time-series measurement data from the address to the impact is measured.

  Note that the golf club swing operation generally proceeds in the order of address, top, impact, and finish. The address means an initial state in which the head 41 of the golf club 4 is arranged near the ball, as shown in FIG. 4A, and the top means the golf club 4 from the address, as shown in FIG. 4B. Is taken back and the head 41 is swung up most. As shown in FIG. 4 (C), the impact means a state at the moment when the golf club 4 is swung down from the top and the head 41 collides with the ball, and the finish is as shown in FIG. 4 (D). It means a state in which the golf club 4 is swung forward after impact.

In the measurement process, it is preferable that the golf club 4 described above is tried a plurality of times, preferably five times or more. In this case, an average value of the measurement data can be calculated and used for subsequent calculations. Also, in order to remove abnormal values due to miss shots, measurement errors, etc., the standard deviation σ of the measurement data is calculated, and the measurement data of all trial hits are preferably within the average value ± 1.65σ, more preferably the average value ± 1. It is preferable to obtain measurement data that falls within 28σ. Then, in order to perform the check, the control unit 24 calculates the standard deviation σ of the measurement data. When the value of σ does not satisfy the above conditions, a message for requesting addition or re-measurement is displayed. You may make it display on the part 21. FIG. It should be noted that, instead of the average value of the measurement data itself, an average value of machining values calculated based on the measurement data (for example, average shoulder torque T _AVE , arm energy E _AVE and head speed V _h described later) is calculated. May be. Similarly, when calculating the average value of the processed values, it is possible to check the reliability of the data based on the standard deviation σ.

<1-2-2. First conversion step>
Hereinafter, the first conversion step for converting the measurement data of the xyz local coordinate system into the value of the XYZ global coordinate system will be described with reference to FIG. Specifically, first, the acquisition unit 24A performs grip accelerations a _x , a _y , a _z , grip angular velocities ω _x , ω _y , ω _z and grips in the xyz local coordinate system stored in the storage unit 23. geomagnetic m _x, m _y, reads the measurement data of the time series of m _z (step S1).

Next, based on the time-series measurement data in the xyz local coordinate system read out in step S1, the grip behavior deriving unit 24B derives the times t _i , t _t , and t _a of impact, top, and address. (Step S2). In this embodiment, the impact time t _i is first derived, the top time t _t is derived based on the impact time t _i , and the address time t _a is derived based on the top time t _t .

Specifically, the increment of the grip angular velocity ω _x per sampling period Δt is the threshold value 3
The time that first exceeds 00 deg / s is set as the temporary impact time. Then, the time at which the increment of the grip angular velocity ω _x per sampling period Δt exceeds 200 deg / s from the time when a predetermined time is passed from the time of the temporary impact to the time of the temporary impact is detected. The time t _i is set.

Next, the time before the impact time t _i and when the grip angular velocity ω _y is switched from negative to positive is specified as the top time t _t . The address time t _a is calculated according to the flowchart of FIG.

In subsequent step S3, the grip behavior deriving unit 24B calculates a posture matrix N (t) at time t from the address to the impact. Assume that the posture matrix is expressed by the following formula. The posture matrix N (t) is a matrix for converting the XYZ global coordinate system at time t to the xyz local coordinate system.

The meanings of the nine components of the posture matrix N (t) are as follows.
Component a: Cosine of the angle formed by the X axis of the global coordinate system and the x axis of the local coordinate system Component b: Cosine of the angle formed by the Y axis of the global coordinate system and the x axis of the local coordinate system Component c: Overall Cosine of angle formed by Z axis of coordinate system and x axis of local coordinate system Component d: Cosine of angle formed by X axis of global coordinate system and y axis of local coordinate system Component e: Y of global coordinate system Cosine of angle between axis and y axis of local coordinate system Component f: cosine of angle between Z axis of global coordinate system and y axis of local coordinate system Component g: X axis of global coordinate system and local Cosine of angle between z axis of coordinate system Component h: Cosine of angle between Y axis of global coordinate system and z axis of local coordinate system Component i: Z axis of global coordinate system and z of local coordinate system Here, the vector (a, b, c) represents a unit vector in the x-axis direction, and the vector (d, e, f) represents a single vector in the y-axis direction. Represents a vector, the vector (g, h, i) represents a unit vector in the z-axis direction.

Further, the attitude matrix N (t) can be expressed by the following equation according to the concept of Euler angles in the ZYZ system. Here, φ, θ, and ψ are rotation angles around the Z axis, the Y axis, and the Z axis.

In calculating the posture matrix N (t) to the impact from the address, first, the address of the time t _a at the posture matrix N (t _a) is calculated. Specifically, φ and θ at the time of address are calculated according to the following equations. Note that the following expression uses that the golf club 4 is stationary at the time of addressing and only the gravity in the vertical direction is detected by the acceleration sensor 11. The grip accelerations a _x , a _y , and a _z in the following expressions are values at the time of addressing.

Subsequently, ψ at the time of addressing is calculated according to the following equation.
However, the values of m _xi and my _{y in} the above formula are calculated according to the following formula. Further, the following grip geomagnetic m _x in the formula, m _y, m _z is the value at the address.

From the above, to generate address phi, theta, [psi is grip acceleration a _x in the xyz local coordinate system, a _y, a _z and grip geomagnetism m _x, m _y, is calculated on the basis of the m _z. And these phi, theta, by substituting the value Expression 2 of [psi, Adoresu when the posture matrix N (t _a) is calculated.

Subsequently, the posture matrix N (t) from the address to the impact is calculated by updating the posture matrix N (t _a ) at the time of time at intervals of the sampling period Δt. More specifically, first, the attitude matrix N (t) is expressed by the following expression using four quaternion variables q ₁ , q ₂ , q ₃ , q ₄ (q ₄ is a scalar part).

Therefore, from Equations 1 and 7, the four quaternion variables q ₁ , q ₂ , q ₃ , and q ₄ can be calculated according to the following equations.

Now, the value of a~i that defines the address at the time of the attitude matrix N (t _a) is known. Therefore, according to the above formula, first, quaternion four variables q ₁ , q ₂ , q ₃ , and q ₄ at the time of address are calculated.

The quaternion q ′ after a lapse of a minute time from the time t is expressed by the following equation using the quaternion q at the time t.

The first-order differential equation representing the temporal change of the four quaternion variables q ₁ , q ₂ , q ₃ , and q ₄ is expressed by the following equation.

By using the equations (9) and (10), the quaternion at time t can be sequentially updated to the quaternion at the next time t + Δt. Here, the quaternion from the address to the impact is calculated. Then, by sequentially substituting the four quaternion variables q ₁ , q ₂ , q ₃ , and q ₄ from the address to the impact into the formula 7, the attitude matrix N (t) from the address to the impact is calculated. The

Next, in step S4, the grip behavior deriving unit 24B determines the grip acceleration a _x , a _y , a _{z in} the xyz local coordinate system from the address to the impact based on the posture matrix N (t) from the address to the impact. The time series data of the grip angular velocities ω _x , ω _y , and ω _z is converted into time series data in the XYZ global coordinate system. The converted grip accelerations a _X , a _Y , a _Z and grip angular velocities ω _X , ω _Y , ω _Z are calculated according to the following equations.

In subsequent step S5, the grip behavior deriving unit 24B integrates the time series data of the grip accelerations a _X , a _Y , and a _Z , thereby grip speeds v _X , v _{Y in} the XYZ global coordinate system from the address to the impact. , V _Z is derived. At this time, it is preferable to perform offset so that the grip speeds v _X , v _Y , and v _Z from the address to the impact are 0 m / s at the top. For example, offset at an arbitrary time t is grip velocity v _X at time t, v _Y, v from _Z, (grip velocity v _X at the top of the time _{t t, v Y, v Z} ) × t / (t t is performed by subtracting the -t _a).

<1-2-3. Second conversion step>
Hereinafter, the second conversion step for converting the behavior of the grip 42 in the XYZ global coordinate system calculated in the first conversion step into the behavior of the grip 42 in the swing plane P will be described. In the present embodiment, the swing plane P is defined as a plane that includes the origin of the XYZ global coordinate system and is parallel to the Y axis and the shaft 40 of the golf club 4 at the time of impact (see FIG. 7). In the second conversion step, the grip behavior deriving unit 24B projects the grip speeds v _X , v _Y , v _Z and the grip angular velocities ω _X , ω _Y , ω _Z in the XYZ global coordinate system into the swing plane P. v _pX , v _pY , v _pZ and grip angular velocities ω _pX , ω _pY , ω _pZ are calculated.

Specifically, based on the z-axis vector (g, h, i) included in the posture matrix N (t) representing the extending direction of the shaft 40, the golf player 7 is viewed from the front. E) Time series data of the inclination of the shaft 40 is calculated. Then, based on this time series data, the time when the shaft 40 is parallel to the Z axis when viewed from the positive direction of the X axis is specified, and this is set as the impact time t _i . Here, the impact time t _i does not always coincide with the impact time t _i already described. Subsequently, the inclination of the shaft 40 viewed from the negative Y-axis direction is calculated based on the z-axis vector (g, h, i) included in the posture matrix N (t _i ) at the time t _i of the impact. That is, an angle α ′ formed between the shaft 40 and the X axis viewed from the negative direction of the Y axis at the time of impact is calculated, and this is set as the swing plane angle.

When the swing plane angle α ′ is obtained, a projection transformation matrix A for projecting an arbitrary point in the XYZ global coordinate system onto the swing plane P can be calculated as follows. However, α = 90 ° −α ′.

Here, based on the above projective transformation matrix A, the grip speeds v _pX , v _pY , v _pZ and the grip angular velocities ω _pX , ω _pY , ω _pZ after the projective transformation from the address to the impact according to the following formula: Series data is calculated.

Note that the grip speed (v _pY , v _pZ ) obtained by the above calculation represents the grip speed (vector) in the swing plane P, and the grip angular speed ω _pX is an axis perpendicular to the swing plane P. It represents the angular velocity around. Here, the grip speed (scalar) in the swing plane P from the address to the impact is calculated according to the following equation.

Further, here, the inclination of the shaft 40 at the top in the swing plane P, which is necessary for the later calculation, is also calculated. Specifically, first, the z-axis vector (g, h, i) included in the top posture matrix N (t _t ) is projected into the swing plane P using the projective transformation matrix A according to the following equation. To do. Here, the projected vector is (g ′, h ′, i ′).

The vector (h ′, i ′) specified by the above expression is a vector representing the inclination of the shaft 40 at the top in the swing plane P. Therefore, the inclination β of the shaft 40 at the top in the swing plane P is calculated by substituting the above calculation results into the following equation.

<1-2-4. Shoulder behavior derivation process>
Hereinafter, while calculating FIG. 8, shoulder behavior derivation for deriving a pseudo shoulder behavior in the swing plane P based on the grip behavior (grip speed V _GE and grip angular velocity ω _pX ) in the swing plane P will be described below. The process will be described. In the present embodiment, the behavior of the golf club 4 is a double pendulum model in which the shoulder of the golfer 7 and the grip 42 (or the wrist of the golfer who holds the golf club 7) are nodes, and the arm of the golfer 7 and the golf club 4 are linked. Based on the analysis. However, the shoulder behavior is not directly measured, but is derived as a pseudo shoulder behavior based on the actually measured grip behavior. In the following, unless otherwise specified, the term “shoulder” simply means such a pseudo shoulder. The same applies to the pseudo “arm” defined as extending linearly between the pseudo shoulder and the grip 42 (wrist).

In specifying the shoulder behavior from the grip behavior, the double pendulum model according to the present embodiment is based on the following (1) to (5). FIG. 9 is a diagram conceptually illustrating the following preconditions.
(1) On the swing plane P, the grip 42 (wrist) moves circularly around the shoulder.
(2) On the swing plane P, the distance (radius) R between the shoulder and the grip 42 is constant.
(3) The shoulder does not move (but rotates) during the swing motion.
(4) On the swing plane P, the angle between the top arm and the golf club 4 is 90 °.
(5) The arm at the time of impact faces the lower side of the Z axis when viewed from the positive direction of the X axis.

Under the above assumption, the shoulder behavior deriving unit 24C calculates the movement distance D of the grip 42 from the top to the impact in the swing plane P (step S21). The moving distance D is derived by integrating the grip speed V _GE from the top to the impact.

  Subsequently, the shoulder behavior deriving unit 24C calculates the arm rotation angle γ from the top to the impact in the swing plane P (step S22). The rotation angle γ is calculated based on the inclination β of the shaft 40 at the top calculated in the second conversion step. Next, the shoulder behavior deriving unit 24C calculates the radius R = D / γ (step S23).

Then, the shoulder behavior deriving unit 24C calculates an angular velocity (an angular velocity of the arm) ω ₁ around the shoulder from the top to the impact in the swing plane P as the shoulder behavior according to the following equation. That is, the arm angular velocity ω ₁ is a value reflecting the measured grip velocity V _GE .
ω ₁ = V _GE / R

<1-2-5. Index calculation process>
Hereinafter, an index calculation process for calculating a swing index for determining the optimum club weight based on the behavior of the grip 42 and the behavior of the shoulder will be described with reference to FIG. In the present embodiment, average shoulder torque T _AVE , arm energy E _AVE and head speed V _h described later are calculated as swing indices.

Specifically, first, in step S31, the shoulder behavior deriving unit 24C integrates the angular velocity ω ₁ of the arm from the top to the impact, and calculates the rotational angle θ ₁ of the arm from the top to the impact. At this time, it is preferable to use trapezoidal integration. The rotation angle θ ₁ is defined as shown in FIG. 11, and the paper surface of FIG. 11 is equal to the swing plane P. In the following, the analysis proceeds based on a new XY coordinate system in the swing plane P shown in FIG. The X axis of the new XY coordinate system in the swing plane P is equal to the Y axis of the XYZ global coordinate system described above, and the Y axis of the new XY coordinate system is the Z axis of the XYZ global coordinate system within the swing plane P. Is the axis projected onto

Further, the shoulder behavior deriving unit 24C differentiates the angular velocity ω ₁ of the arm from the top to the impact, and calculates the angular acceleration ω ₁ ′ from the top to the impact. Next, the shoulder behavior deriving unit 24C calculates the position (X ₁ , Y ₁ ), the velocity (V _X1 , V _Y1 ) and the acceleration (A _X1 , A _Y1 ) of the center of gravity of the arm from the top to the impact. These values are calculated by substituting the calculation results described above into the following equations.

  Where r is the distance from the shoulder to the center of gravity of the arm. In the present embodiment, it is assumed that the center of gravity of the arm is at the center of the arm. Therefore, R = 2r.

Next, in step S <b> 32, the grip behavior deriving unit 24 </ b> B performs a calculation similar to that in step S <b> 31 also about the grip 42. That is, the grip angular velocity ω _pX from the top to the impact = the angular velocity ω ₂ of the golf club 4 around the grip 42 is integrated, and the rotation angle θ ₂ of the golf club 4 (shaft 40) around the grip 42 from the top to the impact is calculated. To do. Also at this time, it is preferable to use trapezoidal integration, and the rotation angle θ ₂ is defined as shown in FIG.

Subsequently, the grip behavior deriving unit 24B differentiates the angular velocity ω ₂ of the golf club 4 from the top to the impact, and calculates the angular acceleration ω ₂ ′ from the top to the impact. Next, the grip behavior deriving unit 24B calculates the position (X ₂ , Y ₂ ), speed (V _X2 , V _Y2 ), and acceleration (A _X2 , A _Y2 ) of the golf club 4 from the top to the impact. . These values are calculated by substituting the calculation results described above into the following equations.

  However, L is the distance from the grip 42 to the center of gravity of the golf club 4. The value of L is a specification of the golf club 4 and is set in advance.

Next, in step S33, the index calculation unit 24D substitutes the calculation result described above into the following formula to thereby obtain the restraining force R ₂ = (R _X2 , R _Y2 ) generated in the grip 42 from the top to the impact. calculate. The following formula is based on the balance of forces in the translational direction. Where m ₂ is the mass of the golf club and g is the acceleration of gravity. Further, m ₂ is the specification of the golf club 4 and is determined in advance.

In step S34, the index calculation unit 24D, by substituting the following equation calculation results described above, calculates the torque T ₁ and the torque T ₂ of the around the grip 42 of the shoulder around from the top to the impact.

However, I ₁ is the moment of inertia around the center of gravity of the arm, and I ₂ is the moment of inertia around the center of gravity of the golf club 4. In this embodiment, the moment of inertia I ₁ of the arm of the center of gravity around the center of gravity of the arm under the assumption that the center of the arm is calculated as ^{_{I 1 = m 1 r 2/}} 3. m ₁ is the mass of the arm. In the present embodiment, the mass m ₁ of the arm is appropriately determined in advance. For example, before starting the analysis, the weight of the golfer 7 is input, and the weight of the arm is automatically calculated by multiplying the input weight by a predetermined coefficient. I ₂ is the specification of the golf club 4 and is determined in advance.

Further, in the present embodiment, index calculation unit 24D calculates the value T _ti integrated over the interval of the torque T ₁ of the surrounding shoulder from the top to the impact. T _ti is the total amount of torque (hereinafter referred to as total shoulder torque) that is exerted around the shoulder of the golfer 7 from the top to the impact, and in this sense is an index that represents the torque around the shoulder during the swing operation. In addition, the index calculation unit 24D calculates an average shoulder torque (hereinafter referred to as an average shoulder torque) T _AVE = T _ti / t _i −t _c during a swing operation by dividing the torque T _ti by the time from the top to the impact. Is calculated. The average shoulder torque T _{AVE is} also an index representing the torque around the shoulder during the swing operation, and is one of the swing indices. In calculating the total shoulder torque T _ti , only the positive torque T ₁ may be integrated, or the average value of the torque T ₁ may be integrated.

In subsequent step S35, the index calculating unit 24D calculates the arm power (power) E ₁ ′ from the top to the impact based on the above-described calculation result. Specifically, E ₁ 'is the velocity vector of the shoulder and v _s, the velocity vector of the grip 42 as v _g, expressed according to the following equation. Here, R ₁ is a binding force generated on the shoulder. Further, v _s and v _g can be calculated by first-order differentiation of the shoulder position vector d _s and the grip 42 position vector d _g = d _s + (2X ₁ , 2Y ₁ ), respectively.

In this embodiment, since the shoulder does not move, v _s = (0, 0), and the arm power E ₁ ′ is calculated according to the following equation. The index calculation unit 24D calculates the arm power E ₁ ′ from the top to the impact by substituting the calculation result described above into the following equation.

In subsequent step S36, the index calculation unit 24D identifies the time t _c a work rate of the arm with subsequent top E ₁ 'turns from positive to negative, the workload of the arm from time t _t of the top to the time t _c E ₁ is calculated. The work E ₁ of the arm is calculated by integrating the work E ₁ ′ of the arm in the interval from time t t to _{t c} (see FIG. 12). The work amount E ₁ can be considered as an index representing the work amount (energy) exerted by the arm between the times t t to _{t c} , and in this sense, the work amount E ₁ is called arm energy during the swing operation. Can do.

In subsequent step S37, the index calculation unit 24D calculates a cock release timing t _r during the swing operation. The present inventors, through experimentation, head speed V _h at impact as the swing indicator, cook release timing t _r during the swing operation, and the arm energy E ₁ or the average work rate E _AVE = E ₁ / It was found that there is a correlation with t _c −t _t . The average power E _AVE is arm energy that is averaged or consumed per unit time during the swing motion. Therefore, here, in order to calculate the head speed V _h at impact, cock release timing t _r is calculated. In the present embodiment, the cock release timing t _r is the time t _t ~t section arms work rate at up to _c E ₁ 'is the time when the maximum is identified as a cook release timing t _r (see FIG. 12) .

In subsequent step S38, the index calculation unit 24D calculates the head speed V _h at the time of impact based on the cock release timing _tr and the arm energy E _AVE . Specifically, the head speed V _h at the time of impact is calculated according to the following formula. Note that k ₁ , k ₂ , and k ₃ are constants obtained by multiple regression analysis from the results of many experiments performed in advance, and are values that are stored in the storage unit 23 in advance. Thus, the index calculation process ends.
V _h = k ₁ · E _AVE + k ₂ · _tr + k ₃

<1-2-6. Optimal total weight determination process>
Hereinafter, the flow of the optimum total weight determining step for determining the optimum club weight will be described with reference to FIG.

  First, in step S40, the determination unit 24E determines the type of the test club that has been tried in the measurement process. If the professional model club has been trial-struck, the process proceeds to step S41, and if the average model club has been trial-struck, the process proceeds to step S51. Which test club has been tried is determined based on information input by the user via the input unit 22.

Note that the type of the test club is determined in step S40 because the area where the swing index is distributed differs depending on the type of the test club. More specifically, the inventors have 10 golfers who normally use professional model golf clubs (hereinafter referred to as professional model users) and golfers who normally use average model golf clubs (hereinafter referred to as averages). When 10 model users) tried to strike a professional model club and an average model club, and calculated the swing index, the results shown in FIG. 14 were obtained. The swing indices calculated in this experiment are the average shoulder torque T _AVE and the arm energy E _AVE , and specific values were calculated according to the same process as described above. And the present inventors discovered from the result shown in FIG. 14 that the space which shows a swing parameter | index is divided | segmented into a pro model area | region and an average model area | region, as shown in FIG. Note that the pro model region is a region where the swing index at the time of the swing motion by the professional model user is distributed, and the average model region is the region where the swing index at the time of the swing motion by the average model user is distributed. In the example of FIG. 14, the pro model region is a region where the average shoulder torque T _AVE ≧ 55 N · m and the arm energy E _AVE ≧ 150 N · m / s, and the average model region is the average shoulder torque T _AVE <55 N · m. M or arm energy E _AVE <150 N · m / s. However, these numerical values may differ depending on swing conditions such as the type of test club. In this experiment, as a professional model club, a driver of SRIXON (registered trademark) Z-525 (shaft is Miyazaki KENA Blue6 S-Flex, golf club weight is 315 g, balance is D2) manufactured by Dunlop Sports Co., Ltd. XXIO (registered trademark) 7 manufactured by Dunlop Sports Co., Ltd. (shaft is MP-700 R-Flex, golf club weight is 285 g, balance is D1) was used as the average model club. .

Next, step S41 and subsequent steps S42 to S45 will be described. Steps S41 to S45 are steps for determining an optimum club weight range (hereinafter referred to as an optimum weight zone) according to the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h . Here, as the values of T _AVE , E _AVE and V _h are larger, the optimum weight band is set to a larger value stepwise.

Specifically, in step S41, the determination unit 24E determines that the head speed V _h is 45 m / s or more, the average shoulder torque T _AVE is 85 N · m or more, and the arm energy E _AVE is 350 N · m / s. It is determined whether or not s or more (hereinafter referred to as professional condition 1). When professional condition 1 is satisfied, it is determined that 320 g or more is the optimum weight band. On the other hand, if the professional condition 1 is not satisfied in step S41, the process proceeds to step S42. In step S42, the determination unit 24E determines whether or not the average shoulder torque T _AVE is 75 N · m or more and the arm energy E _AVE is 250 N · m / s or more (hereinafter, professional condition 2). . When professional condition 2 is satisfied, it is determined that 310 g to 320 g are optimum weight zones. On the other hand, if the professional condition 2 is not satisfied in step S42, the process proceeds to step S43. In step S43, the determination unit 24E determines whether the average shoulder torque T _AVE is 70 N · m or more and the arm energy E _AVE is 175 N · m / s or more (hereinafter, professional condition 3). . And when professional condition 3 is satisfy | filled, it determines with 305g-315g being an optimal weight belt | band | zone. On the other hand, if the professional condition 3 is not satisfied in step S43, the process proceeds to step S44. In step S44, the determination unit 24E determines whether or not the average shoulder torque T _AVE is 55 N · m or more and the arm energy E _AVE is 120 N · m / s or more (hereinafter, professional condition 4). . And when professional condition 4 is satisfy | filled, it determines with 300g-310g being an optimal weight range. On the other hand, if the professional condition 4 is not satisfied in step S44, the process proceeds to step S45. In step S45, determination unit 24E determines whether the head velocity V _h greater than 38m / s (hereinafter, pro condition 5). When professional condition 5 is satisfied, it is determined that 290 g to 300 g is the optimum weight band. On the other hand, if the professional condition 5 is not satisfied in step S45, it is determined that the average model club is more appropriate than the professional model club. In response, the display control unit 24F displays a message on the display unit 21 indicating that the measurement process is to be performed again using the average model club, and the process proceeds to step S51.

On the other hand, step S51 and subsequent steps S52 and S53 are steps for determining an optimum weight band according to the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h . Here, as the values of T _AVE , E _AVE and V _h are larger, the optimum weight band is set to a larger value stepwise.

Specifically, in step S51, the determination unit 24E has a head speed V _h of 40 m / s or more, an average shoulder torque T _AVE of 65 N · m or more, and an arm energy E _AVE of 200 N · m / s. It is determined whether or not it is equal to or greater than s (hereinafter, average condition 1). When the average condition 1 is satisfied, it is determined that the professional model club is more appropriate than the average model club. In response to this, the display control unit 24F displays a message on the display unit 21 indicating that the measurement process is to be started again using the professional model club, and proceeds to step S41. On the other hand, when the average condition 1 is not satisfied in step S51, the process proceeds to step S52. In step S52, the determination unit 24E determines whether or not the average shoulder torque T _AVE is 55 N · m or more and the arm energy E _AVE is 150 N · m / s or more (hereinafter, average condition 2). . When average condition 2 is satisfied, it is determined that 290 g to 300 g is the optimum weight band. On the other hand, when the average condition 2 is not satisfied in step S52, the process proceeds to step S53. In step S53, the determination unit 24E determines that the head speed V _h is greater than 38 m / s, the average shoulder torque T _AVE is 48 N · m or more, and the arm energy E _AVE is 100 N · m / s or more (hereinafter, average). It is determined whether or not the condition 3) is satisfied. When the average condition 3 is satisfied, it is determined that 290 g to 300 g is the optimum club weight band. On the other hand, if the average condition 3 is not satisfied in step S53, it is determined that 290 g or less is the optimum club weight band.

The above steps S41 to S45 and steps S51 to S53 are based on the following knowledge. That is, in the above-described experiment described with reference to FIGS. 14 and 15, the head speed V _h and the optimum club weight that maximizes the flight distance were also calculated. The head speed V _h was calculated according to the same process as described above. On the other hand, regarding the optimum club weight, golf clubs having various weights were swung to the golfer, the weight of the golf club giving the maximum flight distance was specified, and this was set as the optimum club weight. More specifically, five types of golf clubs of 324 g, 316 g, 312 g, 308 g, and 300 g were swung to the professional model user, and golf clubs of 295 g, 290 g, and 285 g were swung to the average model user. 17 and 18 show the values of the head speed V _h and the optimum club weight obtained by this experiment.

The optimum club weight can also be calculated as the optimum club weight in the simulation according to the following algorithm. That is, the present inventors have found that there is a relationship as shown in FIG. 16 between the head speed V _h and the weight of the golf club 4. Specifically, as the golf club 4 to be swung becomes heavier, the golf club 4 cannot be shaken off by the golfer's force, and the head speed V _h becomes smaller. However, even if it is lighter than a certain level, the head speed V _h reaches its peak (see FIG. 16). This is because it cannot be shaken with a force more than a full swing. In other words, even if it is easy to swing lightly above a certain level, the head speed V _h does not increase when the golfer reaches the limit. Thus, the head speed V _h is a certain point (change point of FIG. 16) as a boundary, a proportional region head speed V _h is proportional to the weight of the golf club 4, head speed V _h is the weight of the golf club 4 Regardless, it is divided into constant areas that are generally constant. Further, in order to extend the flight distance, it is preferable that both the head speed V _h and the weight of the golf club 4 are large. Therefore, the weight of the golf club 4 corresponding to the above change point is the optimum club that maximizes the flight distance. It can be said that it is weight. Accordingly, by swinging a golf club 4 different weight golfers participating in the experiment, head speed V _h and the head speed relationship between the weight of the golf club 4 V _h - plotted in golf club weight plane. Then, a regression line in the proportional region and a regression line in the constant region are calculated, a change point that is an intersection thereof is obtained, and the golf club weight corresponding to the change point can be set as the optimum golf club weight.

As shown in FIGS. 17 and 18, it can be seen from the above experiment that the optimum club weight increases as the arm energy E _AVE and the average shoulder torque T _AVE both increase. As a result, the present inventors calculated the average shoulder torque T _{AVE -arm} energy E _{AVE -head} speed V _h space as shown in FIG. 19 for the professional model user and as shown in FIG. 20 for the average model user. It was discovered that the region showing the optimum weight band can be defined by dividing the region. However, for simplicity, in FIG. 19 and FIG. 20, the axis representing the head speed V _h is omitted, and the average shoulder torque T _{AVE -arm} energy E _AVE plane is shown. That is, in steps S41 to S45 and steps S51 to S53 described above, the points indicating the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h are plotted in any region in the T _AVE -E _AVE -V _h space. Depending on whether or not, it is a step for determining the optimum weight band. Note that the threshold values of T _AVE , E _AVE, and V _h used for the determinations in steps S41 to S45 and steps S51 to S53, in other words, the values indicating the boundaries of the divided areas shown in FIGS. Are stored in the storage unit 23 as correspondence data 28. That is, the correspondence relationship data 28 is data that defines the correspondence between the sizes of T _AVE , E _AVE, and V _h and the optimum weight band. In steps S41 to S45 and steps S51 to S53, the correspondence data 28 in the storage unit 23 is referred to, and the above determination is performed. In FIG. 2, the correspondence relationship data 28 is shown as data different from the fitting program 3, but may be incorporated in the program 3.

By the way, in the professional model area, when T _AVE ≧ 85 N · m and E _AVE ≧ 350 N · m / s, the optimum weight band is 320 g or more or 310 g depending on whether the head speed V _h is 45 m / s or not. It can be distributed to any of ~ 320g. As shown in the experimental results shown in FIG. 17, this is the optimum in the case of V _h ≦ 45 m / s in the vicinity of the boundary line of the region where T _AVE ≧ 85 N · m and E _AVE ≧ 350 N · m / s. This is because when the club weight is 316 g and V _h ≧ 45 m / s, the optimum club weight is 324 g or more.

<1-2-7. Optimal shaft determination process>
Hereinafter, the optimum shaft determination process for calculating the optimum shaft weight and the optimum bending rigidity suitable for the golfer 7 will be described.

First, the determination unit 24E determines an optimum shaft weight range (hereinafter, optimum shaft weight band) based on the optimum weight band. Specifically, as shown in Tables 1 and 2 below, the optimum shaft weight band is determined in advance for each optimum weight band for each model of the golf club 4, and the optimum shaft weight band is determined based on this. .

  Subsequently, the determination unit 24E determines the range of the optimum bending stiffness of the shaft 40 (hereinafter, the optimum bending stiffness band). Since the method for determining the optimum bending stiffness band is known (refer to Japanese Patent Laid-Open No. 2013-208366 if necessary), detailed description thereof is omitted here.

  When the optimal weight band, the optimal shaft weight band, and the optimal bending rigidity band are determined by the above steps, the determination unit 24E determines the weight of the entire golf club 4, the weight of the shaft 40, and the bending rigidity of various golf clubs whose known weights are known. A golf club belonging to the optimum weight zone, the optimum shaft weight zone and the optimum bending stiffness zone is identified from the clubs. The display control unit 24F displays on the display unit 21 the optimum weight band, the optimum shaft weight band, and the optimum bending rigidity band together with information indicating the type of the specified golf club. Thereby, the user can know the type of golf club suitable for the golfer 7, and can know the optimum weight band, the optimum shaft weight band, and the optimum bending rigidity band.

<2. Second Embodiment>
FIG. 21 shows the overall configuration of the fitting system 120 according to the second embodiment. The fitting system 120 is common in many respects to the fitting system 100 according to the first embodiment. Therefore, in the following, for the sake of simplicity, the description will focus on the differences between the two embodiments, the same reference numerals will be given to the same components, and the description will be omitted. The fitting system 120 is also a system that analyzes the swing motion of the golf club 4 based on the measurement data obtained by measuring the swing motion of the golf club 4 by the golfer 7, and is used for supporting the fitting of the golf club 4. The measurement of the swing motion is also performed by the sensor unit 1 attached to the grip 42 of the golf club 4 as in the first embodiment.

  The difference between the first and second embodiments is that an optimum swing MI decision step is mainly executed instead of the optimum total weight decision step. More specifically, in the first embodiment, the optimum club weight is determined as an optimum swinging ease index. However, in the second embodiment, instead of this, the swing inertia moment (hereinafter referred to as the golf club 4) suitable for a golfer is used. , The optimum swing MI) is determined. In the second embodiment, the optimum shaft determination step is also omitted.

The swing moment of inertia is the moment of inertia around the shoulder during a swing, and can be defined according to the following equation, for example. Where I _S is the swing moment of inertia.
I _S = I ₂ + m ₂ (R + L) ² + I ₁ + m ₁ (R / 2) ²
In addition, about each golfer 7, even if the golf club 4 changes, the weight of an arm is the same. Therefore, in this embodiment, for the sake of simplicity, the swing inertia moment is calculated according to the following equation, omitting the inertia moment corresponding to the rotation of the arm.
I _S = I ₂ + m ₂ (R + L) ²
Furthermore, in this embodiment, the swing inertia moment is calculated with the arm length R = 60 cm (constant). However, the value of the arm length R calculated in step S23 can be substituted for R in the above formula. Incidentally, m ₂ , I ₂ , and L, which are parameters for determining I _S, are the specifications of the golf club 4. Therefore, the swing moment of inertia in the present embodiment is also a specification of the golf club 4.

  As shown in FIG. 21, the fitting system 120 includes a fitting device 102 instead of the fitting device 2. The fitting device 102 has the same hardware configuration as the fitting device 2, but the fitting device 102 is installed with a fitting program 103 instead of the fitting program 3. Therefore, the control unit 24 reads and executes the fitting program 103 in the storage unit 23 to virtually acquire the acquisition unit 24A, the grip behavior deriving unit 24B, the shoulder behavior deriving unit 24C, the index calculating unit 24D, and the display control unit. In addition to operating as 24F, it can operate as the determination unit 124E. The determination unit 124E is a virtual unit that executes an optimal swing MI determination process that is a difference from the first embodiment. In addition, correspondence data 128 is stored in the storage unit 23 of the fitting apparatus 102 in place of the correspondence data 28 so that the optimum swing MI determination step can be executed. The correspondence data 128 is data indicating conditions for determining the optimal swing MI.

  Also in the second embodiment, as in the first embodiment, the first conversion process, the second conversion process, the shoulder behavior derivation process, and the index calculation process are sequentially performed, and then the optimum swing MI determination process is performed. Below, the measurement process and the optimal swing MI determination process, which are different from the first embodiment, will be described.

<2-1. Measurement process>
Also in the second embodiment, the measurement process is executed as in the first embodiment. However, in the second embodiment, instead of the two test clubs of the professional model club and the average model club, one test club with the sensor unit 1 is tried by the golfer 7. In other respects, the measurement steps according to the first embodiment and the second embodiment are the same. However, also in the second embodiment, it is possible to improve the fitting accuracy as in the first embodiment by arranging two test clubs.

<2-2. Optimal swing MI determination process>
Hereinafter, the flow of the optimum swing MI determination process will be described with reference to FIG. In the optimal swing MI determination step, the range of the optimal swing MI (hereinafter referred to as the optimal swing MI band) is determined according to the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h . Here, as the values of T _AVE , E _AVE, and V _h are larger, the optimum swing MI band is set to a larger value stepwise.

Specifically, in step S60, the determination unit 124E determines whether or not the average shoulder torque T _AVE is 77 N · m or more and the arm energy E _AVE is 180 N · m / s or more (hereinafter, condition 1). Determine whether. When the condition 1 is satisfied, the determination unit 124E determines whether or not the head speed V _h is 45 m / s or more (hereinafter, condition 2) (step S61). When condition 2 is satisfied, it is determined that 5600 g · cm ² or more is the optimal swing MI band, and otherwise, it is determined that 5590 to 5630 g · cm ² is the optimal swing MI band. On the other hand, when the condition 1 is not satisfied in step S60, the process proceeds to step S62. In step S62, the determination unit 124E determines whether or not the average shoulder torque T _AVE is 60 N · m or more and the arm energy E _AVE is 170 N · m / s or more (hereinafter, condition 3). When the condition 3 is satisfied, the determination unit 124E determines whether or not the head speed V _h is 45 m / s or more (hereinafter, condition 4) (step S63). If condition 4 is satisfied, it is determined that 5590-5630 g · cm ² is the optimal swing MI band, and otherwise, it is determined that 5510-5590 g · cm ² is the optimal swing MI band. On the other hand, when the condition 3 is not satisfied in step S62, the process proceeds to step S64. In step S64, the determination unit 124E determines whether the average shoulder torque T _AVE is 50 N · m or more and the arm energy E _AVE is 130 N · m / s or more (hereinafter, condition 5). If condition 5 is satisfied, it is determined that 5460 to 5510 g · cm ² is the optimal swing MI band, and otherwise, it is determined that 5480 g · cm ² or less is the optimal swing MI band.

The above steps S60 to S64 are based on the following knowledge. That is, an experiment similar to FIG. 14 was performed here. More specifically, the inventors calculated the swing index by causing 21 golfers to try the test club according to the present embodiment, and obtained the results shown in FIG. Note that the swing index calculated in this experiment is also the average shoulder torque T _AVE and arm energy E _AVE , and specific values were calculated according to the same steps as described above. In this experiment, a driver of SRIXON (registered trademark) Z-525 (shaft is Miyazaki Kosuma Blue 6 S-Flex, the weight of the golf club is 314 g, and the balance is D3) manufactured by Dunlop Sports Co., Ltd. is used as a test club. It was. In this experiment, the head speed V _h and the optimum swing MI that maximizes the flight distance were also calculated. The head speed V _h was calculated according to the same process as described above. On the other hand, with respect to the optimum swing MI, golf clubs having various swing inertia moments were swung by the golfer, the swing inertia moment of the golf club giving the maximum flight distance was specified, and this was set as the optimum swing MI. More specifically, five types of golf clubs having swing moments of inertia of 5650 g · cm ² , 5610 g · cm ² , 5550 g · cm ² , 5485 g · cm ² , and 5400 g · cm ² were swung. FIG. 23 shows values of the head speed V _h and the optimum swing MI obtained by this experiment.

Note that the optimum swing MI can also be calculated as an optimum swing MI in simulation according to the following algorithm. That is, the present inventors not only have a weight m _{2 of} the golf club 4 but also a swing index such as the swing inertia moment I _S of the golf club 4 and a swing index such as the head speed V _h. Has been found to have a general relationship as shown in FIG. Specifically, as the swing inertia moment I _S of the golf club 4 to be swung increases, the golf club 4 becomes harder to swing and the head speed V _h decreases. However, even if I _S is made smaller than a certain value, the head speed V _h reaches its peak (see FIG. 24). This is because it cannot be shaken with a force more than a full swing. In other words, even if easily fixed or smaller swing in the swing moment of inertia I _S, if reached the limit of golfers, head speed V _h can not stretch. Thus, the head speed V _h is a certain point (change point of FIG. 24) as a boundary, a proportional region head speed V _h is proportional to the swing moment of inertia I _S, head speed V _h is the swing moment of inertia I _S Regardless, it is divided into constant areas that are generally constant. Further, in order to extend the flying distance, since better head speed V _h is large, or the swing moment of inertia I _S of a golf club 4 corresponding to the change point is the optimum swing MI to the distance to the maximum It can be said. Therefore, golf clubs 4 having various swing inertia moments I _S are swung by golfers participating in the above experiment, and the relationship between the head speed V _h and the swing inertia moment I _S is represented by the head speed V _h -swing inertia moment I _S plane. Plot inside. Then, a regression line in the proportional region and a regression line in the constant region are calculated, a change point that is an intersection thereof is obtained, and the swing inertia moment I _S corresponding to the change point can be set as the optimum swing MI. Note that the same relationship is established between the grip end inertia moment _IG and the head speed V _h described in the third embodiment, and an optimum grip end MI described later is calculated in the same manner based on the relationship of FIG. Is possible.

As shown in FIG. 23, it can be seen from the above experiment that the optimum swing MI increases as the arm energy E _AVE and the average shoulder torque T _AVE both increase. As a result, the present inventors can define an area indicating the optimum swing MI band by dividing the average shoulder torque T _{AVE -arm} energy E _{AVE -head} speed V _h space as shown in FIG. I discovered that there is. However, for simplicity, in FIG. 25, the axis representing the head speed V _h is omitted, and the average shoulder torque T _{AVE -arm} energy E _AVE plane is shown. That is, the above-described steps S60 to S64 are performed according to which region in the T _AVE -E _AVE -V _h space the points indicating the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h are plotted. In this step, the optimum swing MI band is determined. Note that the threshold values of T _AVE , E _AVE, and V _h used in the determinations in steps S 60 to S 64, in other words, values indicating the boundaries of the divided areas shown in FIG. 25 are stored in the storage unit 23 as correspondence data 128. Stored. That is, the correspondence relationship data 128 is data that defines the correspondence between the magnitudes of T _AVE , E _AVE, and V _h and the optimum swing MI band. In steps S60 to S64, the correspondence data 128 in the storage unit 23 is referred to, and the above determination is performed. In FIG. 21, the correspondence relationship data 128 is shown as data different from the fitting program 103, but may be incorporated in the program 103.

By the way, when T _AVE ≧ 77 N · m and E _AVE ≧ 180 N · m / s (condition C1), or when T _AVE ≧ 60 N · m and E _AVE ≧ 170 N · m / s , when the condition is not satisfied C1, due whether head speed V _h is 45 m / s or more, the optimum swing MI band 5600 g · cm ² or more or 5590~5630g · cm ² or 5510~5590g · cm ² Cano It can be assigned to either. This is because, as shown in the experimental results shown in FIG. 23, in these regions, the optimum swing MI is different depending on whether V _h ≧ 45 m / s or not.

When the optimal swing MI band is determined by the above steps, the determination unit 124E specifies a golf club belonging to the optimal swing MI band from various golf clubs whose swing MI is known. In the storage unit 23, the swing inertia moment of the golf club or a value necessary for calculating the swing inertia moment is associated with information (manufacturer, model number, etc.) specifying a number of golf clubs in advance. Stores information indicating specifications including. The display control unit 24F displays the optimum swing MI band on the display unit 21 together with information indicating the type of the specified golf club. Thereby, the user can know the type of golf club suitable for the golfer 7 and can know the optimum swing MI band.
<3. Third Embodiment>
FIG. 26 shows the overall configuration of the fitting system 130 according to the third embodiment. The fitting system 130 is common in many respects to the fitting systems 100 and 120 according to the first and second embodiments, particularly the fitting system 120. Therefore, for the sake of simplicity, the following description will focus on the differences between the third embodiment and the second embodiment, and the same components as those in the first and second embodiments will be denoted by the same reference numerals. Is omitted. The fitting system 130 is also a system that analyzes the swing motion of the golf club 4 based on the measurement data obtained by measuring the swing motion of the golf club 4 by the golfer 7, and is used for the purpose of supporting the fitting of the golf club 4. The measurement of the swing motion is also performed by the sensor unit 1 attached to the grip 42 of the golf club 4 as in the first and second embodiments.

  The difference between the third embodiment and the second embodiment is that an optimum grip end MI determining step is mainly executed instead of the optimum swing MI determining step. More specifically, in the second embodiment, the optimum swing MI is determined as the optimum swing ease index. However, in the third embodiment, instead of this, the grip end moment of inertia of the golf club 4 suitable for a golfer ( Hereinafter, the optimum grip end MI) is determined. Also in the third embodiment, the optimum shaft determination step is omitted.

The grip end moment of inertia is the moment of inertia around the grip end, and is calculated according to the following equation in this embodiment. Where _IG is the grip end moment of inertia.
I _G = I ₂ + m ₂ L ²
Incidentally, m _2, I _2, L is a parameter which determines the I _G is the specifications of the golf club 4. Therefore, the grip end moment of inertia in the present embodiment is also a specification of the golf club 4.

  As shown in FIG. 26, the fitting system 130 includes a fitting device 202 instead of the fitting device 102. The fitting device 202 has the same hardware configuration as the fitting device 102, but in the fitting device 202, a fitting program 203 is installed instead of the fitting program 103. Therefore, the control unit 24 reads and executes the fitting program 203 in the storage unit 23 to virtually acquire the acquisition unit 24A, the grip behavior deriving unit 24B, the shoulder behavior deriving unit 24C, the index calculating unit 24D, and the display control unit. In addition to operating as 24F, it can operate as the determination unit 224E. The determination unit 224E is a virtual unit that executes an optimal grip end MI determination process that is a difference from the second embodiment. In addition, correspondence data 228 is stored in the storage unit 23 of the fitting device 202 in place of the correspondence data 128 so that the optimum grip end MI determination step can be executed. The correspondence data 228 is data indicating conditions for determining the optimum grip end MI.

  In the third embodiment, as in the second embodiment, the measurement process, the first conversion process, the second conversion process, the shoulder behavior derivation process, and the index calculation process are sequentially performed, and then the optimum grip end MI determination process is performed. The Below, the optimal grip end MI determination process which is a difference with 2nd Embodiment is demonstrated.

<3-1. Optimal grip end MI determination process>
As can be seen by comparing FIGS. 22 and 27, the optimum grip end MI determination step is the same as the optimum swing MI determination step except that the optimum swing ease index that is finally distributed is different. Therefore, although detailed description is omitted, as shown in FIG. 27, in the optimum grip end MI determination step, the optimum grip end is determined according to the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h. The range of MI (hereinafter, optimum grip end MI band) is determined. Further, as the values of T _AVE , E _AVE and V _h are larger, the optimum grip end MI band is set to a larger value stepwise.

Although it is clear if the definition equations of the grip end inertia moment and the swing inertia moment are compared, both can be converted into each other if the arm length R is known. Here, returning to FIG. 23, the legend on the right side of FIG. 23 also shows a value obtained by converting the optimum swing MI into the optimum grip end MI. From the above, it can be seen that the optimum grip end MI increases as the arm energy E _AVE and the average shoulder torque T _AVE both increase. Therefore, the average shoulder torque T _{AVE -arm} energy E _{AVE -head} speed V _h space can be divided into regions indicating the optimum grip end MI band as shown in FIG. However, for simplicity, in FIG. 28, the axis representing the head speed V _h is omitted, and the average shoulder torque T _{AVE -arm} energy E _AVE plane is shown. That is, in steps S70 to S74 of the optimum grip end MI determination step, the points indicating the average shoulder torque T _AVE , arm energy E _AVE, and head speed V _h are the same as in the second embodiment. T _AVE −E _AVE −V _h The optimum grip end MI band is determined in accordance with which region in the space is plotted. Note that the threshold values of T _AVE , E _AVE, and V _h used for the determinations in steps S 70 to S 74 of this embodiment, in other words, values indicating the boundaries of the divided areas shown in FIG. 28 are stored as correspondence data 228. It is stored in the unit 23. That is, the correspondence relationship data 228 is data that defines the correspondence relationship between the sizes of T _AVE , E _AVE, and V _h and the optimum grip end MI band. In steps S70 to S74 of the present embodiment, the correspondence relationship data 228 in the storage unit 23 is referred to, and the optimum grip end MI band is determined. In FIG. 26, the correspondence relationship data 228 is shown as data different from the fitting program 203, but may be incorporated in the program 203.

  When the optimum grip end MI band is determined by the above steps, the determination unit 224E specifies a golf club belonging to the optimum grip end MI band from among various golf clubs whose grip end moments of inertia are known. In the storage unit 23, it is necessary to calculate the grip end inertia moment or the grip end inertia moment of the golf club in advance in association with information (manufacturer, model number, etc.) specifying a large number of golf clubs. Information indicating specifications including various values is stored. The display control unit 24F displays the optimum grip end MI band on the display unit 21 together with information indicating the type of the specified golf club. Thus, the user can know the type of golf club suitable for the golfer 7 and can know the optimum grip end MI band.

  By the way, in FIG. 23, as described above, the results of an experiment in which 21 golfers actually tried a test club are plotted. Of the 21 swing data, 18 of the 18 people except for the 3 marked F are the optimum swing MI and the optimum grip end MI. The results belong to the optimal swing MI band and the optimal grip end MI band determined in the determination step. That is, it was confirmed that fitting was possible with a correct answer rate of 86%.

<4. Modification>
As mentioned above, although one Embodiment of this invention was described, this invention is not limited to the said embodiment, A various change is possible unless it deviates from the meaning. For example, the following changes can be made. Moreover, the gist of the following modifications can be combined as appropriate.

<4-1>
In the above-described embodiment, the sensor unit 1 including the acceleration sensor, the angular velocity sensor, and the geomagnetic sensor is used as the measurement device that measures the swing motion of the golfer 7, but the measurement device may have other configurations. . For example, the geomagnetic sensor can be omitted. In this case, the measurement data can be converted from the xyz local coordinate system to the XYZ global coordinate system by a statistical method. In addition, since such a method is a well-known technique (refer to Unexamined-Japanese-Patent No. 2013-56074 if needed), detailed description is abbreviate | omitted here. Alternatively, a three-dimensional measurement camera can be used as a measurement device. Since a method of measuring the behavior of a golfer, a golf club, or a golf ball with a three-dimensional measurement camera is also known, detailed description thereof is omitted here. If a 3D measurement camera is used, the process of converting the measurement data from the xyz local coordinate system to the XYZ global coordinate system can be omitted, and the grip behavior in the XYZ global coordinate system can be measured directly. can do.

<4-2>
The head speed V _h can be calculated by other methods than the statistical method based on the multiple regression equation described above. For example, you may calculate geometrically according to the following formula | equation. However, L _club is the length of the golf club which is the spec of the golf club.
Position vector (d _hX , d _hY ) of the head at the tip of the shaft 40
d _hX = 2X ₁ + L _club cos θ ₂
d _hY = 2Y ₁ + L _club sin θ ₂
Speed vector (V _hX , V _hY ) of the head at the tip of the shaft 40
V _hX = 2V _X1 -L _club ω ₂ sin θ ₂
V _hY = 2V _Y1 + L _club ω ₂ cos θ ₂
V _h = sqrt (V _hX ² + V _hY ² )

<4-3>
The optimal total weight determination process, the optimal swing MI determination process, and the optimal grip end MI determination process can be executed by the user, not the fitting devices 2, 102, 202. That is, the calculated swing index is displayed on the display unit 21 after the index calculation process is completed. Then, based on the correspondence data 28, 128, and 228 as shown in FIGS. 19, 20, 25, and 28, the user himself / herself determines the optimum weight zone, the optimum swing MI zone, and the optimum grip end MI zone. Also good. At this time, if the correspondence data 28, 128, 228 is displayed on the display unit 21 or a paper medium on which the correspondence data 28, 128, 228 is printed is prepared, the user's judgment becomes easy.

<4-4>
The swing index shown in the above embodiment is an exemplification, and various other indexes having a certain relationship (correlation) with the optimal swing ease index such as the optimal club weight, the optimal swing MI, and the optimal grip end MI are shown. Can be used. For example, not only the arm energy E _AVE and the average shoulder torque T _AVE but also the arm energy E ₁ and the shoulder torque T _ti can be used as swing indices. In addition, an index representing the energy of an arbitrary part or the torque of an arbitrary part exhibited by the golfer 7 during the swing operation can be used as the swing index. Further, the term “during swing operation” here is not limited to the time from the top to the impact or cock release timing described in the above embodiment, and there is a certain relationship (correlation) between the swing index and the optimum club weight. As long as it is, it can be said to mean any time or period during the swing operation. In addition, the following indices that are considered to have a certain relationship (correlation) with the energy or torque exerted by the golfer 7 can also be used as the swing indices. In addition, the indexes other than the following (7) to (9) tend to increase the optimal swing ease index as the value increases, and the indexes (7) to (9) are optimal as the value decreases. There is a tendency for the ease of swinging index to increase.
(1) Angle θ3 formed by the shaft 40 and the overall coordinate system Z-axis (below the grip) at the top time (see FIG. 29)
(2) Average value of angular velocity ω ₂ during swing motion (3) Maximum value of angular velocity ω ₂ from top to impact (4) Average value of grip speed V _GE from top to impact (5) From top Maximum value of the grip speed V _GE to impact (6) Distance D of grip 42 between top and impact D
(7) cock release timing t _r say the difference (in this case between the time of cock release timing t _r and the impact is, Hayamari the release speed of the cock angle θ4 the arm and the shaft 40 is made, the arm of the energy of the shaft 40 energy It can be defined as the timing when it starts to change.
(8) cock angle θ4 to cook release timing t arm and the shaft 40 at _r is made (see Figure 29)
(9) Down swing time, that is, time from top to impact (10) Integral value of torque T ₂ from time of top to time when torque T ₂ around grip 42 reverses positive and negative

<4-5>
In the above embodiment, the optimal swing such as the optimal club weight, the optimal swing MI, and the optimal grip end MI is selected according to the magnitudes of the three swing indexes (head speed V _h , arm energy E _AVE and average shoulder torque T _AVE ). Although the ease index is determined, the optimum swing ease index may be determined according to the size of one, two, four or more swing indices. For example, the optimal swingability index may be determined only according to the magnitude of any one or two indices selected from the three of arm energy E _AVE , average shoulder torque T _AVE, and head speed V _h. it can. The results of the experiments shown in FIGS. 17, 18 and 23 show that the arm energy E _AVE and the average shoulder torque T _AVE have a higher correlation with the optimal swingability index than the head speed V _h . Therefore, when the head speed V _h is used as a swing index, it is preferable to use it in combination with at least one of the arm energy E _AVE and the average shoulder torque T _AVE . Further, in addition to E _AVE , T _AVE , and V _h , an optimum swing ease index is determined according to the size of an arbitrary number of indices selected from the swing indices exemplified in the modified example 4-4. You can also.

<4-6>
In the above-described embodiment, the calculation of the optimal club weight, the optimal swing MI, and the optimal grip end MI is exemplified as the optimal swing ease index. However, an optimum value may be calculated for various other indexes representing the ease of swinging the golf club and used as the optimum swing ease index. For example, the inertia moment I ₂ around the center of gravity of the golf club is also correlated with the swing index, and the optimum value of the inertia moment I ₂ around the center of gravity of the golf club may be calculated as the optimum swing ease index. In the above embodiment, the optimum club weight, the optimum swing MI, or the optimum grip end MI is calculated as the optimum swing ease index. However, by calculating a plurality of optimum swing ease indices, The fitting may be performed based on the ease index. For example, all of the optimum club weight, the optimum swing MI, and the optimum grip end MI may be calculated, and a golf club meeting these three conditions may be searched from the database.

  Examples of the present invention will be described below. However, the present invention is not limited to the following examples.

The average shoulder torque T _AVE , arm energy E _AVE , and head speed V _h were calculated by letting a professional model user and five average model users each hit a professional model club and an average model club, respectively. The result shown in was obtained. The values of T _AVE and E _AVE were calculated from the above measurement process according to the same process as the index calculation process. Further, the value of V _h was measured by a BCS system. The BCS system used here is a measurement system disclosed in Japanese Patent Application Laid-Open Nos. 2012-170547 and 2012-170532 by the applicants. In addition, the optimum weight band was derived by the method according to the first embodiment, and the optimum club weight for maximizing the flight distance was derived by the same method as the experiment described in the first embodiment.

  From the above results, with regard to the professional model user, except for one of the five P3, the optimum weight band according to the present example and the optimum club weight by the experiment coincided. In addition, there was no significant difference between the two values for one person who did not agree. In addition, regarding the average model user, except for 1 of A4 among 5 people, the optimum weight band in the method according to the present example and the optimum club weight by the experiment coincided. In addition, there was no significant difference between the two values for one person who did not agree. Thereby, the high accuracy of the fitting method according to the first embodiment was confirmed.

1 Sensor unit (measuring device)
2,102,202 Fitting device 3,103,203 Fitting program 4 Golf club 7 Golfer 24A Acquisition unit 24B Grip behavior deriving unit 24C Shoulder behavior deriving unit 24D Index calculation unit (calculation unit)
24E, 124E, 224E determination unit 41 head 42 grip

Claims (10)

A fitting device for determining an optimal swingability index that is a golf club swingability index suitable for a golfer,
An acquisition unit for acquiring a measurement value obtained by measuring a swing motion of the test club by the golfer with a measurement device;
A calculation unit that calculates a swing index indicating a characteristic amount of the swing motion based on the measurement value;
A determination unit that determines the optimal swingability index according to the size of the swing index;
Fitting device. The calculation unit calculates a plurality of types of the swing index,
The determining unit determines the optimal swingability index according to the size of the plurality of types of swing indexes.
The fitting device according to claim 1. The swing index includes an index representing energy or torque exerted by the golfer during the swing operation, or an index correlated with these.
The fitting apparatus according to claim 1 or 2. The swing index includes at least one of arm energy of the golfer during the swing operation, torque around a shoulder of the golfer during the swing operation, and head speed.
The fitting device according to claim 3. The calculation unit calculates a head speed during the swing operation,
The determining unit determines the optimal swingability index according to the size of the head speed in addition to the size of at least one of the arm energy and the torque around the shoulder.
The fitting device according to claim 4. The determination unit determines the optimum swing ease index to a larger value as the swing index is larger or smaller.
The fitting apparatus in any one of Claim 1 to 5. For each type of the test club, a storage unit for storing correspondence data that defines a correspondence relationship between the swing index size and the optimum swing ease index size is further provided.
The determining unit determines the optimal swingability index by referring to the correspondence data in the storage unit according to the type of the test club.
The fitting apparatus in any one of Claim 1 to 6. The swing ease indicator includes at least one of a weight of the golf club, a moment of inertia of the golf club, and a moment of inertia around a shoulder of the golfer.
The fitting apparatus in any one of Claim 1 to 7. A fitting method for determining an optimal swingability index, which is a golf club swingability index suitable for a golfer,
Measuring the swing movement of the test club by the golfer with a measuring device;
Calculating a swing index indicating a characteristic amount of the swing motion based on a measurement value obtained by the measuring device;
Determining the optimal swingability index according to the size of the swing index,
Fitting method. A fitting program for determining an optimal swingability index that is a golf club swingability index suitable for a golfer,
Obtaining a measurement value obtained by measuring the swing motion of the test club by the golfer with a measuring device;
Calculating a swing index indicating a characteristic amount of the swing motion based on a measurement value obtained by the measuring device;
Causing the computer to execute a step of determining the optimal swingability index according to the size of the swing index;
Fitting program.