CN109030014B - Prediction method for subjective scoring of noise in vehicle during vehicle acceleration - Google Patents
Prediction method for subjective scoring of noise in vehicle during vehicle acceleration Download PDFInfo
- Publication number
- CN109030014B CN109030014B CN201810572704.4A CN201810572704A CN109030014B CN 109030014 B CN109030014 B CN 109030014B CN 201810572704 A CN201810572704 A CN 201810572704A CN 109030014 B CN109030014 B CN 109030014B
- Authority
- CN
- China
- Prior art keywords
- noise
- mark
- kurtosis
- subjective
- acceleration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a method for predicting subjective scoring of noise in a vehicle during vehicle acceleration, which comprises the following steps: step 1: collecting and storing signals of the sound pressure level of the acceleration noise and the sound pressure level of multi-order acceleration noise varying with the rotating speed of the engine in the acceleration process of the vehicle; step 2: calculating noise kurtosis according to the current engine speed and the current sound pressure level, and calculating normalized noise kurtosis and maximum normalized noise kurtosis through the noise kurtosis; and step 3: and calculating the noise subjective prediction scores corresponding to different rotating speeds in the vehicle acceleration process and the total value of the synthesized subjective prediction score in the whole acceleration process according to the noise kurtosis, and obtaining the acceleration noise evaluation condition according to a score standard table. The invention is based on real vehicle measurement, uses a determined objective evaluation method, greatly reduces uncertainty for traditional subjective evaluation, can quickly and accurately identify problems or compare nuances, and has strong reference.
Description
Technical Field
The invention relates to an acoustics/noise measurement method in the automobile field, in particular to a prediction method for subjective scoring of noise in an automobile when the automobile accelerates.
Background
At present, the popularization rate of automobiles is very high, the acceleration performance of the automobiles is better and better along with the continuous development of automobile technology, but in the acceleration process of the automobiles, the high-speed rotation of an engine inevitably generates noise; for noise evaluation in an automobile acceleration process, the existing method mainly carries out noise evaluation by subjectively driving a vehicle, the evaluation has high consumption of manpower and material resources, meanwhile, the evaluation result is influenced by a plurality of influence factors, so that the deviation is large, particularly, when the sample size is small, the deviation of the subjective evaluation method is large due to the subjective tendency of people, and therefore, the referenceability of the existing scoring method on problem identification or performance comparison is greatly reduced.
Disclosure of Invention
The invention aims to provide a method for predicting subjective scoring of noise in a vehicle during vehicle acceleration, which is based on real vehicle measurement and uses a determined objective evaluation method, so that the uncertainty is greatly reduced for the traditional subjective evaluation, problems can be quickly and accurately identified or subtle differences can be compared, and the reference is strong.
The invention is realized by the following steps:
a prediction method for subjective scoring of noise in a vehicle during vehicle acceleration comprises the following steps:
step 1: collecting and storing signals of the sound pressure level of the total acceleration noise beside the ears and the sound pressure level of multi-order acceleration noise changing along with the rotating speed of the engine in the acceleration process of the vehicle;
step 2: calculating noise kurtosis according to the current engine speed and the current sound pressure level, and calculating normalized noise kurtosis and maximum normalized noise kurtosis through the noise kurtosis;
and step 3: and calculating the noise subjective prediction scores corresponding to different rotating speeds in the vehicle acceleration process and the total value of the synthesized subjective prediction score in the whole acceleration process according to the noise kurtosis, and obtaining the acceleration noise evaluation condition according to a score standard table.
In the step 2, the method further comprises the following sub-steps:
step 2.1: a point is taken as A on any one noise curve in the sound pressure level of the medium acceleration noise and the sound pressure level of the multi-order acceleration noisen(n,Lpn) Wherein n is AnEngine speed of point, LpnA sound pressure level corresponding to An;
step 2.2: respectively to AnBilateral extension of dots, polypeptide, and nucleic acidkThe rotating speed is taken as two extended rotating speed points An-½△nk(n-½△nk,Lpn-½△nk)、An+½△nk(n+½△nk,Lpn+½△nk) The difference between the two extended rotation speed points is delta nkObtaining sound pressure levels of the corresponding rotating speeds of the two extended rotating speed points, wherein k is any rotating speed and is an integer;
step 2.3: taking Deltan when calculating rotating speed nkThe noise kurtosis of the time is calculated according to the following formula:
Fnk=[ Lpn-½(Lpn-½△nk + Lpn+½△nk)];
step 2.4: calculating the normalized noise kurtosis when the rotating speed is n, wherein the calculation formula is as follows:
Fnkr=Fnk/(△nk/△nr)
wherein Δ nrIs a normalized rotation speed;
step 2.5: calculating the maximum normalized noise kurtosis when the rotating speed is n, wherein the calculation formula is as follows:
Fnrmax=MAX(Fnrk)。
in the step 3, the method further comprises the following sub-steps:
step 3.1: calculating the subjective prediction score of the noise when the rotating speed corresponding to each noise curve in the step 2 is n:
when F is presentnrmax <At time 0:
Markn=10-0.5*| Fnrmax|,(10-0.5*| Fnrmax|≥0);
Markn=0,(10-0.5*| Fnrmax|<0);
when F is presentnrmaxWhen the ratio is more than or equal to 0:
Markn=10-Fnrmax,(10-Fnrmax≥0);
Markn=0,(10-Fnrmax<0);
step 3.2: and (3) calculating a synthesized subjective evaluation prediction score according to the following calculation formula:
wherein L ispn-OAMark as total noise sound pressure leveln-OAPredicting scores for subjective evaluation corresponding to the total noise sound pressure level; markn-morderFor an arbitrary order noise sound pressure level Lpn-morderCorresponding subjective evaluation prediction scores, wherein m is the order, and m belongs to (0, + ∞); sigma represents that m-order noise which participates in the operation arbitrarily is operated and summed in a formula;
step 3.3: mark subjective scoring corresponding to multiple rotating speedsntThe average value and the sample estimation standard deviation value of (a) are calculated to predict the total value of the score for the synthetic subjective evaluation in the whole acceleration process, and the calculation formula is as follows:
Markt=average(Markn1t,Markn2t,Markn3t,…)-2*stdev(Markn1t,Markn2t,Markn3t,…)。
according to the method, the degree of the noise causing the human auditory sense complaining caused by the noise changing along with the rotating speed is directly evaluated through analyzing the objective data of the accelerated noise in the vehicle, the normalized noise kurtosis is taken as a key physical quantity, the score of the change of the human auditory sense along with the rotating speed is further calculated, and finally the total evaluation score of the kurtosis complaining is calculated for the noise in the whole accelerating process.
The method can quickly obtain a result consistent with subjective evaluation by adopting the normalized noise kurtosis, can be used for evaluating the noise level of the whole vehicle along with the change of the rotating speed, can also transversely compare the quality of noise in different accelerated vehicles, provides a tool for efficiently identifying problems for the NVH performance development of the whole vehicle, and can be used for setting an optimization target.
The invention is based on real vehicle measurement, uses a determined objective evaluation method, greatly reduces uncertainty for traditional subjective evaluation, can quickly and accurately identify problems or compare nuances, and has strong reference.
Drawings
FIG. 1 is an acceleration noise curve diagram of the sound pressure level of the noise in the vehicle corresponding to the rotation speed in the prediction method for subjective scoring of the noise in the vehicle when the vehicle is accelerated according to the present invention;
FIG. 2 is a graph showing a prediction curve of the in-vehicle noise and the subjective score of the prediction method of the subjective score of the in-vehicle noise when the vehicle is accelerated according to the present invention;
fig. 3 is a flowchart of a method for predicting subjective scoring of noise in a vehicle when the vehicle is accelerated according to the present invention.
Detailed Description
The invention is further described with reference to the following figures and specific examples.
Referring to fig. 3, a method for predicting subjective score of noise in a vehicle during vehicle acceleration includes the following steps:
step 1: signals of the total acceleration noise sound pressure level and the multi-order acceleration noise sound pressure level which are weighted by the A at the ear in the acceleration process of the vehicle and change along with the rotating speed of the engine are collected and stored, and a graph of the noise in the vehicle is shown in an attached figure 1.
Step 2: and calculating the noise kurtosis according to the current engine speed and the current sound pressure level, and calculating the normalized noise kurtosis and the maximum normalized noise kurtosis through the noise kurtosis.
Step 2.1: a point is taken as A on any one of the noise curves of the total acceleration noise sound pressure level and the multi-order acceleration noise sound pressure leveln(n,Lpn) Wherein n (r/min) is AnEngine speed of point, Lpn(dB (A)) is the sound pressure level corresponding to An point.
Step 2.2: respectively to AnBilateral extension of dots, polypeptide, and nucleic acidkThe rotating speed is taken as two extended rotating speed points An-½△nk(n-½△nk,Lpn-½△nk)、An+½△nk(n+½△nk,Lpn+½△nk) The difference between the two extended rotation speed points is delta nkAnd obtaining sound pressure levels of the corresponding rotating speeds of the two extended rotating speed points, wherein k is any integer and is the code of any rotating speed.
Step 2.3: taking Deltan when calculating rotating speed nkThe noise kurtosis of the time is calculated according to the following formula:
Fnk=[ Lpn-½(Lpn-½△nk + Lpn+½△nk)]。
step 2.4: calculating the normalized noise kurtosis when the rotating speed is n, wherein the calculation formula is as follows:
Fnkr=Fnk/(△nk/△nr)
wherein Δ nrIs a normalized rotation speed.
Step 2.5: calculating the maximum normalized noise kurtosis when the rotating speed is n, wherein the calculation formula is as follows:
Fnrmax=MAX(Fnrk)。
and step 3: and calculating the noise subjective scores corresponding to different rotating speeds in the vehicle acceleration process and the total value of the synthesized subjective evaluation prediction scores in the whole acceleration process according to the noise kurtosis, and obtaining the acceleration noise evaluation condition according to a score standard table.
Step 3.1: calculating the subjective prediction score of the noise when the rotating speed corresponding to each noise curve in the step 2 is n:
when F is presentnrmax <At time 0:
Markn=10-0.5*| Fnrmax|,(10-0.5*| Fnrmax|≥0);
Markn=0,(10-0.5*| Fnrmax|<0);
when F is presentnrmaxWhen the ratio is more than or equal to 0:
Markn=10-Fnrmax,(10-Fnrmax≥0);
Markn=0,(10-Fnrmax<0)。
the above calculation method is applicable to various acceleration noise curves. To distinguish between total (OA) and different order (order) accelerated noise sound pressure levels and corresponding scores, the subscripts are used to distinguish and represent: l ispn-OAMark as total noise sound pressure leveln-OAPredicting scores for subjective evaluation corresponding to the total noise sound pressure level; markn-morderFor an arbitrary order noise sound pressure level Lpn-morderCorresponding subjective evaluation prediction score, m is order, m ∈ (0, + ∞);
Step 3.2: and (3) calculating the comprehensive effect of the total noise (OA) and each order noise (order), namely synthesizing the subjective evaluation prediction score, wherein the calculation formula is as follows:
and sigma represents that any m-order noise participating in the operation is subjected to operation in the formula and summed. Mark to distinguish between noise subjective prediction scores calculated from single curvesnThe subscript increases by t.
Finally, the score can be plotted as a curve of the variation of the rotation speed, and the curve of the noise in the vehicle and the subjective score prediction is shown in the attached figure 2.
Step 3.3: mark subjective scoring corresponding to multiple rotating speedsntThe average value and the sample estimation standard deviation value of the total score of the synthetic subjective prediction in the whole acceleration process are calculated, and the calculation formula is as follows:
Markt=average(Markn1t,Markn2t,Markn3t,…)-2*stdev(Markn1t,Markn2t,Markn3t,…)
wherein n is taken to be all the rotation speeds tested in the rotation speed section.
Example 1:
taking normalized rotation speed delta nr=300 r/min, the interval of the revolution speed of the upspeed engine is 25r/min, and the following calculated revolution speed is n =2100r/min and Δ nk1Normalized noise kurtosis of =50r/min, 50r/min being the step size.
As can be seen from table 3, at a rotational speed of n =2100r/min, the sound pressure level is 67.54, the rotational speeds at the two spread rotational speed points are 2125r/min and 2075r/min, respectively, and the corresponding sound pressure levels are 67.99 and 65.36, respectively.
AnNoise kurtosis F of a pointnk=[ Lpn-½(Lpn-½△nk + Lpn+½△nk)]=67.54-½*(67.99+65.36)=0.865。
AnGrouping of pointsNormalized noise kurtosis Fnkr=Fnk/(△nk/△nr)=0.865/(50/300)=5.19。
Taking the following steps according to the step length of 50 r/min:
△nk2=100 r/min, corresponding Fnkr=4.85
△nk3=150 r/min, corresponding Fnkr=3.66
…
△nk20=1000 r/min, F corresponding theretonkr=1.57
20 normalized noise kurtosis F can be calculatednkrTo obtain the maximum normalized noise kurtosis Fnkrmax=5.19。
Due to FnkrmaxNot less than 0 and 10-FnkrmaxThe rotation speed is more than or equal to 0, so the subjective prediction score of the rotation speed of 2100r/min is as follows:
Markn=10-Fnrmax=10-5.19=4.81。
since this score is a subjective prediction score corresponding to OA, it is noted Mark2100-OA=4.81。
Referring to table 3, corresponding to the corresponding order noise at the rotation speed of 2100r/min, the subjective prediction scores of 2 th order, 4 th order, 6 th order and 8 th order are calculated in the same way:
Lp2100-2order=66.90dB(A),Mark2100-2order=4.06;
Lp2100-4order=49.09dB(A),Mark2100-4order=8.74;
Lp2100-6order=50.31dB(A),Mark2100-6order=6.40;
Lp2100-8order=40.01dB(A),Mark2100-8order=9.37。
the subjective noise prediction evaluation at the rotating speed of 2100r/min comprises the following steps:
the calculation at other rotational speeds is the same as the above calculation method.
In the whole acceleration process, the total value of the combined subjective evaluation and prediction score is calculated according to the following formula;
Markt=average(Markn1,Markn2,Markn3,…,Markn20)-2*stdev(Markn1,Markn2,Markn3,…,Markn20)=9.32-1.15*2=7.02。
for the present example, the subjective evaluation was performed by 10 testers at a rotation speed of 2100r/min and the whole acceleration process, and the actual subjective scores are shown in table 1 below:
personnel code | 2100r/min actual subjective assessment score | Actual total subjective evaluation score |
Number 1 | 4 | 7.5 |
|
5 | 7 |
No. 3 | 4.5 | 6.5 |
|
4 | 7 |
Number 5 | 5 | 6.5 |
|
4 | 6.5 |
No. 7 | 4.5 | 7 |
Number 8 | 5 | 6.5 |
Number 9 | 5 | 7 |
|
4.5 | 7 |
Mean value | 4.55 | 6.85 |
Standard deviation of | 0.44 | 0.34 |
Wherein, the scoring criteria are shown in table 2:
comparing the subjective evaluation scores, wherein the difference value between the 2100r/min synthesized subjective prediction score (4.26) obtained by calculation through the method and the mean value (4.55) of the actual subjective evaluation score is 0.29, the value is smaller than the standard deviation of the actual subjective evaluation score, the standard deviation is that the standard uncertainty is 0.44, and the synthesized subjective prediction score is equal to the actual subjective score in terms of statistics, so that the predicted value is used for representing the actual subjective score, and the method has considerable scientific, reference and actual application values.
Similarly, the same is true for the evaluation of noise in the whole acceleration process, the difference value between the total value of the synthetic noise subjective prediction score (7.02) and the total value of the actual subjective score (6.85) is 0.17, the value is smaller than the standard deviation of the total value of the actual subjective evaluation score, the standard deviation is 0.34, and the synthetic noise subjective prediction score is equal to the total value of the actual subjective score in a statistical sense, so that the actual subjective score is represented by a predicted value, and the method has considerable scientific, reference and practical application values.
The present invention is not limited to the above embodiments, and any modifications, equivalent replacements, improvements, etc. within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Table 3: relation table of engine speed, total acceleration noise sound pressure level (OA-level) and multi-order acceleration noise sound pressure level (m order)
Lp/dB(A) | Lp/dB(A) | Lp/dB(A) | Lp/dB(A) | Lp/dB(A) | |
Engine speed/r/min | OA-level | 2 order | 4 order | 6 order | 8 order |
1400 | 54.76 | 47.54 | 49.57 | 36.06 | 37.24 |
1425 | 56.96 | 51.97 | 49.07 | 39.75 | 39.03 |
1450 | 56.97 | 51.97 | 48.11 | 39.98 | 40.43 |
1475 | 56.98 | 51.93 | 46.36 | 40.96 | 42.60 |
1500 | 57.31 | 52.63 | 43.12 | 41.32 | 48.52 |
1525 | 57.45 | 52.99 | 41.34 | 41.13 | 49.04 |
1550 | 57.59 | 53.49 | 40.05 | 40.65 | 49.46 |
1575 | 57.61 | 53.70 | 40.67 | 40.89 | 45.20 |
1600 | 57.54 | 53.59 | 42.23 | 41.35 | 45.19 |
1625 | 57.31 | 52.85 | 44.63 | 40.78 | 45.98 |
1650 | 57.25 | 52.37 | 45.49 | 37.52 | 45.05 |
1675 | 57.28 | 52.07 | 45.85 | 35.36 | 44.28 |
1700 | 57.44 | 52.00 | 45.93 | 34.41 | 44.81 |
1725 | 57.84 | 52.39 | 45.93 | 35.46 | 45.78 |
1750 | 58.14 | 52.77 | 46.51 | 36.68 | 46.11 |
1775 | 58.54 | 53.19 | 47.23 | 38.99 | 43.96 |
1800 | 58.67 | 53.29 | 47.56 | 40.91 | 41.98 |
1825 | 58.73 | 53.30 | 47.70 | 42.70 | 39.83 |
1850 | 58.70 | 53.17 | 47.41 | 43.65 | 39.73 |
1875 | 58.50 | 52.65 | 46.76 | 44.66 | 41.70 |
1900 | 58.35 | 52.25 | 46.03 | 44.44 | 42.79 |
1925 | 58.44 | 52.53 | 45.61 | 44.03 | 43.15 |
1950 | 58.82 | 53.71 | 45.66 | 48.59 | 43.17 |
1975 | 60.97 | 58.34 | 47.17 | 50.80 | 41.69 |
2000 | 61.92 | 59.85 | 47.59 | 51.41 | 40.00 |
2025 | 62.70 | 61.01 | 47.80 | 51.59 | 39.52 |
2050 | 63.48 | 62.08 | 48.13 | 51.70 | 39.94 |
2075 | 65.36 | 64.43 | 48.67 | 51.33 | 39.68 |
2100 | 67.54 | 66.90 | 49.09 | 50.31 | 40.01 |
2125 | 67.99 | 67.39 | 49.09 | 48.09 | 40.13 |
2150 | 68.37 | 67.81 | 49.23 | 47.28 | 40.03 |
2175 | 68.72 | 68.17 | 49.43 | 47.30 | 39.96 |
2200 | 69.18 | 68.64 | 50.32 | 47.85 | 39.49 |
2225 | 69.68 | 69.14 | 52.04 | 48.85 | 39.71 |
2250 | 69.93 | 69.37 | 52.70 | 49.49 | 41.38 |
2275 | 70.16 | 69.59 | 53.71 | 50.07 | 42.25 |
2300 | 70.22 | 69.64 | 53.91 | 50.01 | 42.46 |
2325 | 70.29 | 69.69 | 54.18 | 50.41 | 43.41 |
2350 | 70.40 | 69.76 | 54.60 | 50.56 | 43.61 |
2375 | 70.00 | 69.07 | 55.86 | 51.20 | 43.60 |
2400 | 69.80 | 68.81 | 55.85 | 50.91 | 43.44 |
2425 | 69.66 | 68.61 | 55.94 | 49.27 | 43.55 |
2450 | 69.37 | 68.22 | 56.00 | 47.56 | 44.18 |
2475 | 68.89 | 67.57 | 55.83 | 46.43 | 43.53 |
2500 | 68.06 | 66.44 | 55.35 | 45.60 | 44.32 |
2525 | 67.70 | 65.96 | 55.28 | 45.93 | 45.42 |
2550 | 67.51 | 65.70 | 55.22 | 45.80 | 46.52 |
2575 | 67.31 | 65.44 | 55.16 | 45.58 | 46.73 |
2600 | 67.05 | 65.06 | 54.90 | 46.04 | 47.48 |
2625 | 66.93 | 64.87 | 55.31 | 46.19 | 48.07 |
2650 | 67.00 | 64.90 | 56.03 | 45.62 | 48.43 |
2675 | 67.09 | 64.95 | 56.61 | 45.63 | 48.67 |
2700 | 67.21 | 65.02 | 57.27 | 45.58 | 47.98 |
2725 | 67.35 | 65.11 | 57.92 | 45.71 | 47.57 |
2750 | 67.49 | 65.19 | 58.62 | 46.28 | 47.20 |
2775 | 67.74 | 65.33 | 59.66 | 45.94 | 46.28 |
2800 | 68.03 | 65.44 | 60.83 | 46.00 | 45.34 |
2825 | 68.11 | 65.45 | 61.17 | 45.66 | 44.64 |
2850 | 68.19 | 65.44 | 61.54 | 45.36 | 45.87 |
2875 | 68.23 | 65.41 | 61.77 | 44.79 | 46.98 |
2900 | 68.25 | 65.27 | 62.13 | 44.48 | 47.65 |
2925 | 68.09 | 64.89 | 62.30 | 43.20 | 49.73 |
2950 | 68.00 | 64.73 | 62.23 | 42.78 | 50.24 |
2975 | 67.82 | 64.43 | 62.01 | 41.97 | 51.09 |
3000 | 67.43 | 63.85 | 61.40 | 40.68 | 52.09 |
3025 | 67.12 | 63.37 | 60.78 | 41.13 | 52.78 |
3050 | 66.80 | 62.87 | 60.01 | 41.44 | 53.06 |
3075 | 66.50 | 62.35 | 59.26 | 41.45 | 53.16 |
3100 | 66.22 | 61.78 | 58.44 | 42.19 | 52.67 |
3125 | 66.04 | 61.31 | 57.92 | 42.62 | 51.99 |
3150 | 65.74 | 60.07 | 57.28 | 42.80 | 50.45 |
3175 | 65.70 | 59.54 | 57.27 | 43.22 | 49.90 |
3200 | 65.72 | 58.87 | 57.83 | 43.97 | 49.13 |
3225 | 65.84 | 58.39 | 58.61 | 45.11 | 48.93 |
3250 | 66.00 | 58.24 | 59.23 | 45.99 | 48.85 |
3275 | 66.20 | 58.34 | 59.77 | 46.69 | 49.26 |
3300 | 66.37 | 58.59 | 60.19 | 47.21 | 49.45 |
3325 | 66.77 | 59.54 | 60.83 | 48.11 | 50.68 |
3350 | 66.98 | 60.13 | 61.07 | 48.48 | 51.47 |
3375 | 67.29 | 61.00 | 61.29 | 49.04 | 51.81 |
3400 | 67.63 | 61.94 | 61.42 | 49.39 | 52.00 |
3425 | 68.04 | 63.02 | 61.39 | 50.31 | 52.03 |
3450 | 68.35 | 63.81 | 61.27 | 50.82 | 51.80 |
3475 | 68.66 | 64.56 | 61.08 | 51.28 | 51.26 |
3500 | 69.21 | 65.79 | 60.64 | 51.83 | 50.40 |
3525 | 69.55 | 66.49 | 60.31 | 52.21 | 50.16 |
3550 | 69.80 | 67.00 | 60.00 | 52.30 | 49.79 |
3575 | 69.96 | 67.33 | 59.75 | 52.38 | 49.69 |
3600 | 70.23 | 67.83 | 59.28 | 52.11 | 49.40 |
3625 | 70.52 | 68.36 | 58.62 | 52.08 | 49.01 |
3650 | 70.76 | 68.77 | 57.69 | 51.61 | 49.19 |
3675 | 70.83 | 68.87 | 57.00 | 51.27 | 49.96 |
3700 | 70.84 | 68.86 | 56.61 | 51.13 | 50.07 |
3725 | 70.80 | 68.77 | 56.17 | 50.59 | 50.44 |
3750 | 70.71 | 68.53 | 55.70 | 50.44 | 50.55 |
3775 | 70.52 | 68.08 | 55.33 | 50.43 | 50.22 |
3800 | 70.36 | 67.66 | 55.15 | 50.91 | 50.32 |
3825 | 70.15 | 67.05 | 55.12 | 50.65 | 51.23 |
3850 | 69.93 | 66.34 | 55.25 | 51.23 | 52.39 |
3875 | 69.70 | 65.42 | 55.63 | 51.62 | 54.74 |
3900 | 69.66 | 65.15 | 55.78 | 51.94 | 55.76 |
3925 | 69.65 | 64.98 | 55.90 | 52.18 | 57.22 |
3950 | 69.69 | 64.95 | 56.06 | 52.43 | 57.70 |
3975 | 69.83 | 65.13 | 56.28 | 52.66 | 58.26 |
4000 | 70.15 | 65.72 | 56.54 | 53.07 | 58.67 |
Claims (1)
1. A method for predicting subjective scoring of noise in a vehicle during vehicle acceleration is characterized by comprising the following steps: the method comprises the following steps:
step 1: collecting and storing signals of the sound pressure level of the total acceleration noise beside the ears and the sound pressure level of multi-order acceleration noise changing along with the rotating speed of the engine in the acceleration process of the vehicle;
step 2: calculating noise kurtosis according to the current engine speed and the current sound pressure level, and calculating normalized noise kurtosis and maximum normalized noise kurtosis through the noise kurtosis;
step 2.1: a point is taken as A on any one of the noise curves of the total acceleration noise sound pressure level and the multi-order acceleration noise sound pressure leveln(n,Lpn) Wherein n is AnEngine speed of point, LpnIs AnThe corresponding sound pressure level;
step 2.2: respectively to AnTwo-sided spread 1/2 Δ n of pointskThe rotating speed is taken as two extended rotating speed points An-1/2△nk(n-1/2△nk,Lpn-1/2△nk)、An+1/2△nk(n+1/2△nk,Lpn+1/2△nk) The difference between the two extended rotation speed points is delta nkObtaining sound pressure levels of the corresponding rotating speeds of the two extended rotating speed points, wherein k is any rotating speed and is an integer;
step 2.3: taking Deltan when calculating rotating speed nkThe noise kurtosis of the time is calculated according to the following formula:
Fnk=[Lpn-1/2(Lpn-1/2△nk+Lpn+1/2△nk)];
step 2.4: calculating the normalized noise kurtosis when the rotating speed is n, wherein the calculation formula is as follows:
Fnkr=Fnk/(△nk/△nr)
wherein Δ nrIs a normalized rotation speed;
step 2.5: calculating the maximum normalized noise kurtosis when the rotating speed is n, wherein the calculation formula is as follows:
Fnrmax=MAX(Fnrk);
and step 3: calculating noise subjective prediction scores corresponding to different rotating speeds in the vehicle acceleration process and a total value of the synthesized subjective prediction scores in the whole acceleration process according to the noise kurtosis, and obtaining the acceleration noise evaluation condition according to a score standard table;
in the step 3, the method further comprises the following sub-steps:
step 3.1: calculating the subjective prediction score of the noise when the rotating speed corresponding to each noise curve in the step 2 is n:
when F is presentnrmax<At time 0:
Markn=10-0.5*|Fnrmax|,10-0.5*|Fnrmax|≥0;
Markn=0,10-0.5*|Fnrmax|<0;
when F is presentnrmaxWhen the ratio is more than or equal to 0:
Markn=10-Fnrmax,10-Fnrmax≥0;
Markn=0,10-Fnrmax<0;
step 3.2: and (3) calculating a synthesized subjective evaluation prediction score according to the following calculation formula:
wherein L ispn-OAMark as total noise sound pressure leveln-OAPredicting scores for subjective evaluation corresponding to the total noise sound pressure level; markn-morderFor an arbitrary order noise sound pressure level Lpn-morderCorresponding subjective evaluation prediction scores, wherein m is the order, and m belongs to (0, + ∞); sigma represents that m-order noise which participates in the operation arbitrarily is operated and summed in a formula;
step 3.3: mark subjective scoring corresponding to multiple rotating speedsntMean and sample estimation ofAnd calculating a total score of the synthesized subjective evaluation prediction in the whole acceleration process by the standard deviation value, wherein the calculation formula is as follows:
Markt=average(Markn1t,Markn2t,Markn3t,…)-2*stdev(Markn1t,Markn2t,Markn3t,…)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810572704.4A CN109030014B (en) | 2018-06-06 | 2018-06-06 | Prediction method for subjective scoring of noise in vehicle during vehicle acceleration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810572704.4A CN109030014B (en) | 2018-06-06 | 2018-06-06 | Prediction method for subjective scoring of noise in vehicle during vehicle acceleration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109030014A CN109030014A (en) | 2018-12-18 |
CN109030014B true CN109030014B (en) | 2020-11-27 |
Family
ID=64612142
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810572704.4A Active CN109030014B (en) | 2018-06-06 | 2018-06-06 | Prediction method for subjective scoring of noise in vehicle during vehicle acceleration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109030014B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110296847B (en) * | 2019-04-17 | 2021-06-01 | 开沃新能源汽车集团有限公司 | Noise testing method for new energy passenger car |
CN111751119B (en) * | 2020-06-10 | 2022-11-11 | 中国第一汽车股份有限公司 | Automobile acceleration sound quality evaluation method based on sound order frequency characteristics |
CN111780979B (en) * | 2020-07-01 | 2022-06-24 | 浙江吉利汽车研究院有限公司 | Sound quality scoring method and sound quality scoring system for engine |
CN112525544B (en) * | 2020-11-12 | 2022-07-01 | 东风汽车集团有限公司 | Comparison method of vehicle acceleration performance parameters and related equipment |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393246A (en) * | 2011-10-22 | 2012-03-28 | 重庆长安汽车股份有限公司 | Noise evaluation method for automobile generator under finished automobile state |
CN103323101A (en) * | 2013-06-04 | 2013-09-25 | 格特拉克(江西)传动系统有限公司 | Method of evaluating transmission gear noise quantificationally and objectively by means of pitch |
CN103471709A (en) * | 2013-09-17 | 2013-12-25 | 吉林大学 | Method for predicting noise quality of noise inside passenger vehicle |
CN104697786A (en) * | 2015-03-17 | 2015-06-10 | 重庆长安汽车股份有限公司 | Objective evaluation method of transmission whine noise |
CN105021275A (en) * | 2015-07-15 | 2015-11-04 | 吉林大学 | Time domain dynamic characteristic analysis-based intra-vehicle sound quality prediction method |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4676111A (en) * | 1986-01-13 | 1987-06-30 | Donaldson Company, Inc. | Method and apparatus for sound testing mufflers |
JP3033342B2 (en) * | 1992-05-22 | 2000-04-17 | 横河電機株式会社 | Cavitation detector |
KR100938217B1 (en) * | 2008-03-11 | 2010-01-22 | 아주대학교산학협력단 | Drivability measurement and analysis system |
DE102009004607A1 (en) * | 2009-01-14 | 2009-08-27 | Daimler Ag | Air noises analyzing method for passenger car, involves determining interval sound level from sound pressure spectrum for frequency interval, and comparing interval sound level with desired value depending on measuring point |
CN105300706A (en) * | 2015-10-22 | 2016-02-03 | 华晨汽车集团控股有限公司 | NVH performance test method for automobile acoustic package |
CN107796620B (en) * | 2017-10-13 | 2019-09-10 | 重庆长安汽车股份有限公司 | A method of evaluation vehicle warming-up operating condition NVH performance |
CN107664579A (en) * | 2017-11-20 | 2018-02-06 | 哈尔滨东安汽车发动机制造有限公司 | Transmission gear, which is uttered long and high-pitched sounds, engages order extraction and decision method |
CN108020427B (en) * | 2017-11-21 | 2020-04-24 | 朱乐慈 | Pure electric vehicle gear shifting quality evaluation method based on GA-BP neural network |
-
2018
- 2018-06-06 CN CN201810572704.4A patent/CN109030014B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393246A (en) * | 2011-10-22 | 2012-03-28 | 重庆长安汽车股份有限公司 | Noise evaluation method for automobile generator under finished automobile state |
CN103323101A (en) * | 2013-06-04 | 2013-09-25 | 格特拉克(江西)传动系统有限公司 | Method of evaluating transmission gear noise quantificationally and objectively by means of pitch |
CN103471709A (en) * | 2013-09-17 | 2013-12-25 | 吉林大学 | Method for predicting noise quality of noise inside passenger vehicle |
CN104697786A (en) * | 2015-03-17 | 2015-06-10 | 重庆长安汽车股份有限公司 | Objective evaluation method of transmission whine noise |
CN105021275A (en) * | 2015-07-15 | 2015-11-04 | 吉林大学 | Time domain dynamic characteristic analysis-based intra-vehicle sound quality prediction method |
Also Published As
Publication number | Publication date |
---|---|
CN109030014A (en) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109030014B (en) | Prediction method for subjective scoring of noise in vehicle during vehicle acceleration | |
CN111521406B (en) | High-speed wind noise separation method for passenger car road test | |
CN111832462B (en) | Frequency hopping signal detection and parameter estimation method based on deep neural network | |
CN110186556B (en) | New energy automobile motor bench test noise evaluation method | |
CN115062678B (en) | Training method of equipment fault detection model, fault detection method and device | |
CN106611599A (en) | Voice recognition method and device based on artificial neural network and electronic equipment | |
CN113378927A (en) | Clustering-based self-adaptive weighted oversampling method | |
CN111413925A (en) | Machine tool fault prediction method based on sound signals | |
CN112131662A (en) | Passenger car wind noise subjective evaluation objective quantification method | |
CN111367174A (en) | Linear quadratic form control improvement method based on convolutional neural network vibration identification | |
CN101163239A (en) | Novel vector quantization inceptive code book generating method | |
CN108631817B (en) | Method for predicting frequency hopping signal frequency band based on time-frequency analysis and radial neural network | |
CN107292275B (en) | Frequency domain division human face feature recognition method and system | |
CN1420486A (en) | Voice identification based on decision tree | |
CN114329806B (en) | BP neural network-based engine rotor bolt assembly quality assessment method | |
CN111667836B (en) | Text irrelevant multi-label speaker recognition method based on deep learning | |
Ko et al. | Data augmentation and squeeze-and-excitation network on multiple dimension for sound event localization and detection in real scenes | |
CN110532512B (en) | Equipment failure mode diagnosis characteristic parameter analysis method based on big data analysis | |
CN110781814A (en) | Signal classification method, device and medium based on Gaussian mixture neural network model | |
Ueno et al. | A proposal of regression hybrid modeling for combining random forest and x-means methods | |
CN111462479A (en) | Traffic flow prediction method based on Fourier-recurrent neural network | |
CN108460405A (en) | A kind of image latent writing analysis Ensemble classifier optimization method based on deeply study | |
CN114925721B (en) | Signal sorting method based on U-Net segmentation network | |
CN107451330B (en) | Accelerated transient simulation method for random telegraph noise in circuit | |
CN110555370B (en) | Channel effect inhibition method based on PLDA factor analysis method in underwater target recognition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |