TAILIEUCHUNG - Random Numbers part 5

if (n != nold) { If n has changed, then compute useful quantien=n; ties. oldg=gammln(en+); nold=n; } if (p != pold) { If p has changed, then compute useful quantipc=; ties. plog=log(p); pclog=log(pc); pold=p; } sq=sqrt(*am*pc) | 296 Chapter 7. Random Numbers if n nold en n oldg gammln en nold n if p pold pc plog log p pclog log pc pold p sq sqrt am pc do do angle PI ran1 idum y tan angle If n has changed then compute useful quantities. If p has changed then compute useful quantities. The following code should by now seem familiar rejection method with a Lorentzian comparison function. em sq y am while em em en Reject. em floor em Trick for integer-valued distribution. t sq y y exp oldg-gammln em -gammln en-em em plog en-em pclog while ran1 idum t Reject. This happens about times per devi- bnl em if p pp bnl n-bnl return bnl ate on average. Remember to undo the symmetry transformation. See Devroye 2 and Bratley 3 for many additional algorithms. CITED REFERENCES AND FURTHER READING Knuth . 1981 SeminumericalAlgorithms 2nd ed. vol. 2 of TheArtofComputerProgramming Reading MA Addison-Wesley pp. 120ff. 1 Devroye L. 1986 Non-Uniform Random Variate Generation New York Springer-Verlag . 2 Bratley P. Fox . and Schrage . 1983 A Guide to Simulation New York SpringerVerlag . 3 . Generation of Random Bits The C language gives you useful access to some machine-level bitwise operations such as left shift . This section will show you how to put such abilities to good use. The problem is how to generate single random bits with 0 and 1 equally probable. Of course you can just generate uniform random deviates between zero and one and use their high-order bit . test if they are greater than or less than . However this takes a lot of arithmetic there are special-purpose applications such as real-time signal processing where you want to generate bits very much faster than that. One method for generating random bits with two variant implementations is based on primitive polynomials modulo 2. The theory of these polynomials is beyond our scope although and will give you small tastes of it . Here Sample page from NUMERICAL RECIPES IN C THE .

• Kỹ thuật lập trình
• tư liệu
• tài liệu lập trình
• thủ thuật lập trình
• ngôn ngữ C#
• ngôn ngữ lập trình
• thành phần cấu tạo vật liệu
• phân loại vật liệu compozit
• công nghệ chế tạo vật liệu compozit
• sản phẩm chế tạo từ vật liệu compozit
• phương hướng phát triển
• Tài liệu điện tử
• Vật liệu điện tử
• Linh kiện điện tử
• Vật liệu linh kiện
• Mạch điện tử
• Đề tài nghiên cứu khoa học
• Tổng hợp
• đặc trưng vật liệu nano
• Vật liệu nano cacbon
• Nano cacbon từ vỏ cua
• Xử lý môi trường
• Bài giảng Y học
• Alcaloid và Dược liệu chứa alcaloid
• Phân lập alc từ dược liệu
• Phân loại alcaloid
• Tính chất lý hóa
• Cấu trúc hóa học
• Tự động học
• Điều khiển học
• Điều khiển tự động
• Lý thuyết hệ thống
• Tài liệu tự động hóa