Simulation and Optimization in Materials TechnologyView this Special Issue
Research Article | Open Access
Integrating the TRIZ and Taguchi's Method in the Optimization of Processes Parameters for SMT
SMT is an assembly technology for core circuit board parts. Unless process parameters are effectively controlled, poor solderability may result in a decline in product quality. This study looks at an SMT manufacturing process in a multinational company. First, the TRIZ contradiction matrix is revised to investigate the association between the 39 parameters in the contradiction matrix and 13 parameters that influence the unevenness of solder paste in the solder paste printing process. Expert verification is then used to screen the key factors affecting the quality of SMT, which are then combined with Taguchi's method to identify the optimal parameter set influencing the thickness of SMT solder paste. Results. TRIZ identifies squeegee pressure, ejection speed, squeegee speed, and squeegee angle as the four parameters with the greatest influence on SMT solder paste thickness. Taguchi's method is used to identify the optimum levels set for the experimental factors and carry out confirmation experiments. The S/N ratio improved from 21.732 db to 26.632 db, while the mean also improved from the current 0.163 mm to 0.155 mm, close to the target value of 0.15 mm. The results show that applying TRIZ and Taguchi's method for the purpose of product improvement is feasible.
Today, product cycles are short and undergo rapid changes. Consumer electronic products, such as smart phones, notebook computers, and digital cameras, play an important role in our everyday lives. As consumers become more demanding, there is an increasing emphasis on lighter, thinner, high quality, and low-priced products that can be delivered to the consumer quickly. In addition, due to global competition, businesses are increasingly demanding rigorous quality standards to meet customer’s demand [1, 2].
To meet the demand for light, thin, short, and small electronic products, print circuit boards (PCB) have evolved from single-layer to multilayer boards, with a consequent reduction in their size. As a result, most traditional plated-through hole parts are gradually falling out of use. An effective solution to this is to attach electronic components onto the PCB, resulting in the development of surface mount technology (SMT). However, unless process parameters are effectively controlled, poor solderability may result in a decline in product quality. For example, if insufficient solder paste is deposited at the PCB printing stage, the strength of the solder joints may be inadequate or an empty solder phenomenon may occur. However, if too much solder paste is deposited, bridges may form between the solder joints, leading to short circuits.
Yang et al.  point out that SMT has become the primary tool for PCB assembly. An SMT quality problem may cause a significant production loss, and SMT quality issues are crucial to competition between companies. In the research on optimizing stencil printing, Tsai  offers a comparative perspective, finding that in the SMT manufacturing process 60% of stencil quality issues are caused by defects in solder paste printing. In their research using Taguchi’s method to improve the parameters of solder paste printer processes, Lai and Wang  use optimized printing parameters, leading to an improvement of approximately 20% in solder paste printing quality. Recent studies on solder paste printing quality by Tsai  and Lai and Wang  have both produced significant improvements in solder paste printing quality. However, they fail to explain how they selected their experimental variables. Many researchers, including Yildiz [6, 7], Durgun and Yildiz , and Yildiz and Solanki , have applied Taguchi’s method or algorithms to plan experiments and parameter optimization in the past few years.
Andersen and Fagerhaug  argue that the TRIZ contradiction matrix is a tool for eliminating root causes of problems in the engineering process and can be used to identify engineering parameters that cause engineering problems or conflicts. This study randomly selected 30 sets of data from the process-site and measured and calculated the process capability index as . This figure does not meet Kane’s  recommended minimum acceptable process capability index of , showing that there is still significant room for improvement in the quality of the solder paste printing process. Therefore, this study attempts to combine TRIZ and Taguchi’s method to carry out empirical research. It is hoped that a predictive model of the manufacturing process outcomes can be used to identify optimized process parameters.
2. Literature Review
TRIZ is an abbreviation for the “theory of inventive problem solving” in the Russian language. The theory was developed by the Soviet inventor Altshuller  in 1946 with the aim of solving different types of contradictions.
Altshuller  developed the concept of the “technical contradiction” after observing a large number of invention proposals. The technical contradiction arises because improving one technical system parameter impacts negatively on another one. For example, making a product lighter by reducing its thickness also makes the product more prone to being easily damaged . However, using better materials increases the cost of manufacturing the product. When confronted with a technical contradiction, a contradiction matrix can normally be used to analyze the problem. Thus, the contradiction matrix is the key analytical tool of TRIZ, as well as a basic technique for technical deconstruction and innovation.
The aim of TRIZ is to avoid conflicts between different elements. Hence, Altshuller identified 39 engineering parameters that often produce technical contradictions. In the matrix, each cell indicates the principles used to resolve these contradictions. The matrix provides a fast and simple way to find solutions to technical contradictions. The matrix is a 39 × 39 matrix. To resolve the contradictions, Altshuller  proposed a set of 40 inventive principles. For each contradiction, several principles are suggested to resolve the contradiction. Selecting relevant inventive principles can help the user to produce a solution to the contradictions.
Loh et al.  argue that TRIZ is a method of knowledge extraction that can be applied systematically to resolve issues in the area of innovation and improvement. Nakagawa  argues that TRIZ is an advanced methodology, which uses a contradiction matrix and 40 inventive principles to provide clear answers to complex problems.
Figure 1 shows an example of a contradiction matrix. In the example, for the improving feature of the area of a moving object (5), the worsening feature is the volume of the moving object (7). From the contradiction matrix, we can identify four inventive principles (4, 7, 14, and 17) to resolve the contradiction. However, if we choose the improving and worsening parameters of productivity (39), we find that the contradiction matrix shows a blank cell, indicating the absence of inventive principles that can provide a solution to the problem. For information on the contradiction matrix and inventive principles, please see http://www.innovation-triz.com/TRIZ40/.
2.2. Taguchi’s Method
Taguchi’s method is derived from traditional experimental design methods. This method was developed by Genichi Taguchi in 1949. When applied in designing communication systems, it enables the number of experiments to be reduced and at the same time identifies problems in such systems . The method advocates using the orthogonal array function combined with a simple function evaluation process to achieve improvements in the manufacturing process and product design .
Taguchi et al.  point out that Taguchi’s method uses parameter design to improve quality; that is, for product target functions that need improvement, the factors and levels that affect the target functions are identified. Subsequently, orthogonal arrays are used to determine the configuration of the experimental factors and the number of experiments in order to obtain the same information provided by a full factorial experiment with a smaller number of experiments, analyzing a small amount of experimental data to increase product quality effectively.
Zhang et al.  and Yildiz  point out that the main tools in Taguchi’s method are orthogonal arrays and ratios, emphasizing the importance of quality issues in the product or during the design and manufacturing process. A tolerance design is used to reduce product performance variation. ANOVA is then used to find the level of influence of each significant factor. The tolerance of each significant factor is set based on the costs of each factor, ensuring that product quality variation is minimal and achieving the most appropriate level of quality, thereby realizing robust design objectives.
Yildiz  points out that Taguchi’s method applies engineering knowledge to experiment planning, focusing on solutions to achieve objectives. As the experimental results obtained by Taguchi’s experimental design methods have a high degree of reproducibility, the configuration of experimental factors is straightforward, the number of experiments required is reduced, and the method of analysis is simple and easy to understand. Because of these advantages, Taguchi’s method has been applied in many industries to improve and optimize design parameters. For instance, Su et al.  use Taguchi’s dynamic approach to increase optical whiteness effectively. Hong  uses Taguchi’s method to identify important factors in market segmentation. In addition, Celani de Souza et al.  have shown that Taguchi’s method can improve the quality of dialysis. Yildiz  uses an immune algorithm and Taguchi’s method to design a new design optimization framework.
Taguchi’s method specifies the objective function as a certain signal-to-noise ratio ( ratio). According to this method, the ratio for different quality characteristics can be separated into larger-the-better (LTB), nominal-the-best (NTB), and smaller-the-better (STB). This study is concerned with improving the uniformity of solder paste thickness in the solder paste printing process. Therefore, NTB is applied out of these quality characteristics. The calculation of the ratio is shown in the following formula (1): where is the sample mean and is the sample standard deviation. The variability characteristic is inversely proportional to the ratio. This means that a larger ratio corresponds to a more robust system.
3. Research Methods
TRIZ and Taguchi’s method are alternative experimental design methods used by enterprises to develop new products and improve product quality. The two approaches are often used separately. The TRIZ contradiction matrix enables technical variables that influence quality characteristics to be quickly identified. Although such an approach is able to identify inventive principles, it can only help users to speculate about solutions. The two methods show a wide variation in their experimental efficiency and additivity. Taguchi’s method can identify a more optimal value from preset factor levels. However, these variables do not necessarily have a significant effect on quality characteristics.
Therefore, this study attempts to combine TRIZ and Taguchi’s method, screening variables that have a significant influence on quality characteristics by linking factors that cause uneven solder paste in the solder paste printing process to the 39 engineering parameters in the TRIZ contradiction matrix, before applying Taguchi’s method to the screened variables to identify the optimum process parameter set. The steps are as follows.(1)Determine experimental variables: list the factors affecting the evenness of solder paste during the solder paste printing process, integrate the 39 engineering parameters in the TRIZ contradiction matrix to create a correlation table and produce a ranking, screening the variables that have a significant effect on quality characteristics.(2)Design and run experiment: use the orthogonal array function from Taguchi’s method for experiment design, the number of repetitions, and conduct the experiment; calculate the ratio and mean. Process managers, engineers, and quality assurance personnel select the levels based on the analysis and discussion of variables that have a significant influence and their range.(3)Optimization analysis: a two-phase optimization analysis is carried out on experimental data to identify the optimum combination and predict the optimum model for the manufacturing process.(4)Maximize the ratio: this study uses the delta value and percent combination suggested by scholars to select important influence factors. We identified experiment factors with a delta value greater than the average effect value and a pooled error smaller than 15% as important influence factors and the point closest to the target level as the optimum level of influence factors.(5)Move the mean closer to the target value: at this stage, the selection has no effect on the ratio. However, factors that have a significant effect on the mean are adjustable, moving the mean closer to the target value.(6)Forecast optimization: applying the additive model, the expected ratio and estimated mean for the optimal combination condition are calculated based on the set levels of the significant influence factors.(7)Confirmation tests: confirmation experiments are run on the optimum combination produced using Taguchi’s method. This result is then compared with the predicted results using Taguchi’s method to confirm the improvement in results.
4. Case Study
4.1. Determine Experiment Variables
This study uses the SMT solder paste printing process for PCB production by a multinational company as a case study to investigate the uniformity of solder paste application. The main purpose of solder paste in PCB is to fix parts to the PCB to ensure that the product functions normally. Solder paste printing is the first stage of the SMT manufacturing process. A stencil and solder paste printer squeeze are used to insert the solder paste onto corresponding pads in the PCB through holes in the board. After removing the stencil, the solder paste is left on the pad in the correct shape, completing the printing process, as shown in Figure 2.
Table 1 shows 13 factors that directly influence the thickness of solder paste during the SMT solder paste printing process used in the production of printed circuit boards as selected by process managers, engineers, and quality assurance personnel in the plant. These are then combined with 39 technical parameters in the TRIZ contradiction matrix to produce a correlation table showing the relationship between the 13 factors and 39 technical parameters. Finally, we screen the squeegee angle, squeegee pressure, squeegee speed, and ejection speed as the four variables that have a significant influence on quality characteristics.
4.2. Design and Run the Experiment
We carry out the experiment on the squeegee angle, squeegee pressure, squeegee speed, and ejection speed as the four variables that have a significant influence on quality characteristics. The experimental design uses the (34) orthogonal array function proposed by Taguchi’s method. The process managers, engineers, and quality assurance personnel select the levels based on the analysis and discussion of their range. Because this study is aimed at improving the uniformity of solder paste thickness during the solder paste printing process, the quality characteristic applied is that of nominal-the-best (NTB). The calculation of the S/N ratio is shown in formula (1).
Based on the influence factor and level settings, nine sets of parameter values are input into the solder paste printer to carry out actual PCB printing. The experiment was repeated four times for each run order and the actual solder paste thickness data was recorded. The target value for solder paste thickness is set at 0.15 mm, with an upper limit of 0.20 mm and a lower limit of 0.10 mm. Table 2 shows the (34) orthogonal array and results for uniformity.
4.3. Optimization Analysis
Carry out optimization analysis on the experimental data and predict the optimal model for the experimental process.
4.3.1. Maximize the Ratio
Analysis of the influence of each experiment factor (, , , and ) on the uniformity was performed with an response table, using a Minitab 16 software package.
Table 3 shows the orthogonal array and associated experimental results for uniformity with calculated ratios. The response table for uniformity is presented in Table 3. It shows the calculated ratios for each level of experimental factors. The experimental factor that has the strongest influence is determined according to the value of delta as shown in Table 3. The value of delta equals the difference between maximum and minimum ratios for a particular experimental factor. The higher the value of delta, the more influential the experimental factor. The experiment factors and their interactions are sorted in relation to the values of delta.
Following Lee , Sheu , and Yang’s  recommendations, this study defines important influence factors as experimental factors that have a value of delta greater than the average effect size. It can be seen from Table 3 that the strongest influence was exerted by squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ).
The ANOVA procedure was used to investigate which design parameters significantly affect quality characteristics. The procedure is performed by separating the total variability of the ratios into contributions by each of the design parameters and the errors. The total variability of the ratio is measured by the sum of the squared deviations from the total mean ratio.
Yildiz , Hong , Hsiang and Lin , and Su and Yeh  use percent contribution to carry out decision making. Percent contribution is the pure sum to squares (purse SS) for each factor as a ratio of the total sum to squares (total SS). When the pooled error percentage () 15%, it can be assumed that no important factors have been omitted from the experiment; in other words, pooled error factors are not significant and can be neglected.
An examination of the calculated percent contribution for all experiment factors also shows a very high influence of factor , factor , and factor on the ratios (Table 4).
First, we find that the percent contributions of squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ) are 34.14%, 31.17%, and 19.81%, respectively. Second, the pooled error is less than 15% (14.88%). Therefore, we can assume that no important factors were missed. By analyzing the experiment factors and ANOVA, this study shows that squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ) have significant effects on the contraction rate.
Based on the above analysis, we are able to determine that squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ) are important influence factors in the present study.
The plots for ratios are shown in Figure 3. The optimal levels for each experiment factor can be easily determined from these graphs in accordance with Taguchi’s “nominal-the-best” performance characteristic. Response graphs show a variation in the ratios when the setting of the experimental factors is changed from one level to another. Figure 3 suggests that the optimum condition for the uniformity is the combination of , , and levels of the respective experiment factors.
4.3.2. Moving the Mean Closer to the Target Value
At this stage, we select appropriate adjustment factors (with no effect on the ratio, but a significant effect on the mean) to move the mean closer to the target value. Table 5 shows the calculated mean for each level of experiment factors. Table 5 shows that squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ) should be selected as factors that have a significant influence on the mean. However, the above three factors also have a significant influence on the ratio. Therefore, we do not select any adjustment factors in this step.
Examination of the calculated percent contribution for all experiment factors also shows a very high influence of factor , factor , and factor on the means: see Table 6.
First, we find that the percent contributions of squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ) are 44.29%, 23.96%, and 23.96%, respectively. From the aforementioned analysis, we are able to determine that squeegee pressure (factor ), squeegee angle (factor ), and ejection speed (factor ) have significant effects on the contraction rate for means.
Based on the previous discussion, this study sets the optimal factor level to .
4.3.3. Forecast Optimization
Based on the previous discussion, this study sets the optimal factor level to . Therefore, the additive model is used to estimate the expected ratio and mean under optimal conditions.
The mean ratio for the nine experiments is , and the predicted ratio under optimal conditions is
Similarly, the mean observation for the nine experiments is , and the predicted mean under optimal conditions is
4.4. Confirmation Tests
The confirmation experiment under optimal conditions produced 25 individual values and 5 ratios. To confirm that optimal conditions can be reproduced, it is necessary to estimate the confidence interval of the and mean. Su  recommends using the formula proposed by Ross  to calculate the confidence interval in order to facilitate confirmation experiments. Some scholars [35, 36] suggest that the ratio should be within ±3 db of the optimal value as the basis for assessing the confidence interval in reproduced experiments. Therefore, combining the aforementioned recommendations, this study uses Ross’ proposed formula to calculate the confidence interval. However, when the calculated confidence interval exceeds ±3 db, the level is set at ±3 db: where CI is the confidence interval, is value at the significance level, , is the significance level, is the degree of freedom of the pooled error, is the variance of the pooled error, is the effective number of observations, and is the number of repeated confirmation experiments.
The ratio and mean confidence intervals for confirmation experiments under optimal conditions are as ; therefore the ratio confidence interval is ±3 db,
The average of the five confirmation experiments is as follows: ratio = 26.632 db and mean = 0.155 mm. These two values fall entirely within the corresponding confidence intervals, indicating that the experiment achieved an improvement. The comparison of values before improvement (current), forecast optimization, and confirmation experiments is shown in Table 7.
Confirmation experiments are run on the optimum combination produced using Taguchi’s method. This result is then compared with the predicted results using Taguchi’s method to confirm the improvement in results.
In this study, the ratio improved from 21.732 db to 26.632 db, while the mean also improved from the current 0.163 mm to 0.155 mm, indicating a reduced variation. In addition, the mean was closer to the target value (0.15 mm), showing that applying TRIZ and Taguchi’s method for improving the uniformity of solder paste thickness in the solder paste printing process is feasible.
This study revised the TRIZ contradiction matrix to investigate the correlation between contradiction matrix parameters and parameters that directly influence the uneven thickness of solder paste in the solder paste printing process, screening the squeegee pressure, ejection speed, squeegee speed, and squeegee angle as the key parameters affecting the quality of SMT solder paste thickness. This is an innovative approach that is empirically shown to be feasible.
Taguchi’s method is used to establish an optimal parameter set from the experimental data, with the prediction error rate reaching the required accuracy and delivering real improvements in process capability and product quality. These improvements can help lower the defect rate and reduce production costs, while shortening delivery times and increasing customer satisfaction. These results may help Taiwan’s SMT assembly factories to increase product quality, explore further different machines and productivity factors, and compare different level parameters to produce even better process parameters for realizing additional quality improvements.
- W.-P. Wang, “Evaluating new product development performance by fuzzy linguistic computing,” Expert Systems with Applications, vol. 36, no. 6, pp. 9759–9766, 2009.
- A. R. Yildiz and K. Saitou, “Topology synthesis of multicomponent structural assemblies in continuum domains,” Journal of Mechanical Design, vol. 133, no. 1, Article ID 011008, 2011.
- T. Yang, T.-N. Tsai, and J. Yeh, “A neural network-based prediction model for fine pitch stencil-printing quality in surface mount assembly,” Engineering Applications of Artificial Intelligence, vol. 18, no. 3, pp. 335–341, 2005.
- T.-N. Tsai, “Modeling and optimization of stencil printing operations: a comparison study,” Computers and Industrial Engineering, vol. 54, no. 3, pp. 374–389, 2008.
- Y.-S. Lai and T. H. Wang, “Optimal design towards enhancement of board-level thermomechanical reliability of wafer-level chip-scale packages,” Microelectronics Reliability, vol. 47, no. 1, pp. 104–110, 2007.
- A. R. Yildiz, “Hybrid taguchi-harmony search algorithm for solving engineering optimization problems,” International Journal of Industrial Engineering: Theory Applications and Practice, vol. 15, no. 3, pp. 286–293, 2008.
- A. R. Yildiz, “Optimization of cutting parameters in multi-pass turning using artificial bee colony-based approach,” Information Sciences, vol. 220, pp. 399–407, 2013.
- I. Durgun and A. R. Yildiz, “Structural design optimization of vehicle components using Cuckoo search Algorithm,” Materials Testing, vol. 54, no. 3, pp. 185–188, 2012.
- A. R. Yildiz and K. N. Solanki, “Multi-objective optimization of vehicle crashworthiness using a new particle swarm based approach,” International Journal of Advanced Manufacturing Technology, vol. 59, no. 1-4, pp. 367–376, 2012.
- B. Andersen and T. Fagerhaug, Root Cause Analysis: Simplified Tools and Techniques, American Society for Quality, Milwaukee, Wis, USA, 2006.
- V. E. Kane, “Process capability indices,” Journal of Quality Technology, vol. 18, no. 1, pp. 41–52, 1986.
- G. S. Altshuller, Creativity as an Exact Science, Gordon & Breach Science, New York, NY, USA, 1984.
- M. Sawaguchi, “A studay of systematic innovation based on an analysis of big hits,” in Proceedings of the ETRIA TRIZ Future Conference, pp. 143–151, Bergmo, Italy, 2010.
- G. Altshuller, 40 Principles: TRIZ Keys to Technical Innovation, Technical Innovation Center, 1998.
- H. T. Loh, C. He, and L. Shen, “Automatic classification of patent documents for TRIZ users,” World Patent Information, vol. 28, no. 1, pp. 6–13, 2006.
- T. Nakagawa, “Overall dataflow structure for creative problem solving in triz/usit,” in Proceedings of the 7th TRIZCON Conference, 2005.
- G. E. P. Box, J. S. Hunter, and W. G. Hunter, Statistics for Experimenters: An Introduction to Design, Data Analysis, and Model Building, John Wiley & Sons, New York, NY, USA, 2nd edition, 2010.
- A. R. Yildiz, “A new hybrid particle swarm optimization approach for structural design optimization in automotive industry,” Journal of Automobile Engineering, vol. 226, pp. 1340–1351, 2012.
- G. Taguchi, S. Chowdhury, and Y. Wu, Taguchi's Quality Engineering Handbook, John Wiley & Sons, New York, NY, USA, 2004.
- J. Z. Zhang, J. C. Chen, and E. D. Kirby, “Surface roughness optimization in an end-milling operation using the Taguchi design method,” Journal of Materials Processing Technology, vol. 184, no. 1-3, pp. 233–239, 2007.
- A. R. Yildiz, “A new hybrid artificial bee colony algorithm for robust optimal design and manufacturing,” Applied Soft Computing, vol. 13, no. 5, pp. 2906–2912, 2013.
- A. R. Yildiz, “Hybrid Taguchi-differential evolution algorithm for optimization of multi-pass turning operations,” Applied Soft Computing Journal, vol. 13, pp. 1433–1439, 2012.
- C.-T. Su, C.-M. Lin, and C.-C. Hsu, “Optimization of the optical whiteness ratio for flexible display by using Taguchi's dynamic approach,” IEEE Transactions on Semiconductor Manufacturing, vol. 25, no. 1, pp. 2–15, 2012.
- C.-W. Hong, “Using the Taguchi method for effective market segmentation,” Expert Systems with Applications, vol. 39, no. 5, pp. 5451–5459, 2012.
- H. J. Celani de Souza, C. B. Moyses, F. J. Pontes et al., “Molecular assay optimized by Taguchi experimental design method for venous thromboembolism investigation,” Molecular and Cellular Probes, vol. 25, no. 5-6, pp. 231–237, 2011.
- A. R. Yildiz, “A new design optimization framework based on immune algorithm and Taguchi's method,” Computers in Industry, vol. 60, no. 8, pp. 613–620, 2009.
- H. H. Lee, Quality Management, Tsang Hai Book Publishing, Taipei, Taiwan, 2009.
- S. H. Sheu, The Management and Control of Quality, Future Career Publishing Corporation, Taipei, Taiwan, 2007.
- S. F. Yang, Quality Management, Hwa Tai Publishing, Taipei, Taiwan, 2006.
- A. R. Yildiz, “Comparison of evolutionary based optimization algorithms for structural design optimization,” Engineering Applications of Artificial Intelligence, vol. 26, no. 1, pp. 327–333, 2013.
- S.-H. Hsiang and Y.-W. Lin, “Optimization of the extrusion process for magnesium alloy sheets using the fuzzy-based Taguchi method,” Arabian Journal for Science and Engineering, vol. 34, no. 1, pp. 175–185, 2009.
- C.-T. Su and C.-J. Yeh, “Optimization of the Cu wire bonding process for IC assembly using Taguchi methods,” Microelectronics Reliability, vol. 51, no. 1, pp. 53–59, 2011.
- C. T. Su, Quality Engineering, Chinese Society for Quality, Taipei, Taiwan, 2009.
- P. J. Ross, Taguchi Techniques for Quality Engineer, McGraw-Hill, New York, NY, USA, 1996.
- Y. M. Chen, Taguchi Methods, Tingmao, Taipei, Taiwan, 2012.
- K. Tatebayashi, Introduction Taguchi Method, Union of Japanese Scientists and Engineers, Tokyo, Japan, 2004.
Copyright © 2013 Yung-Tsan Jou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.