

Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Problem sets from the msc systems engineering course, focusing on modeling and simulation. Topics include calculating averages and standard deviations using monte carlo methods, simulating ellipsoids and moment of inertia, testing gaussian distributions, and simulating radioactive decay and binomial processes. Students are expected to write programs and perform simulations to understand the concepts.
Typology: Exercises
1 / 2
This page cannot be seen from the preview
Don't miss anything!


Modeling & Simulation MSc Systems Engineering, First Semester, 2004
6.1 Use a random number generator of your choice to generate about ten thousand numbers between zero and one. Calculate the average and standard deviation after developing about 10 classes in the interval of equal width. Compare the simulated average with expected mean of 0.5.
6.2 The equation of an ellipsoid centered on the origin is given by the following:
1
2 2 2 ^
c
z b
y a
x.
Simulate the solution for volume using Monte Carlo procedure. The ellipsoid can be inscribed in a solid with six rectangular faces.
6.3 The moment of Inertia of a sphere about z-axis is defined as:
Find the value of Iz for a given (say 2, using Monte Carlo procedure).
6.4 Write a program to test Gaussian distribution of random numbers. How you are going to establish that the simulated distribution is a normal distribution?
6.5 Simulate a model of a radioactive element which decays into a radioactive daughter which subsequently decays into a stable element. Use initial number of parent isotope as 6000 and no daughter elements. Assume the pp = 0.07 and pd = 0.02. Let your domain and range be [0, 300] and [0, 6000] respectively. Explain the following systems through simulations a. How the results change when pd is varied from 0.5pp to 0.2pp? b. What happens when pd = pp (say pp = 0.01)? c. What will be the effect on results when pd = 5pp?
6.6 A drunk person starts off at a lamppost and takes 36 steps each of length 2 and
each in a random direction where 0 o^ 360 o. Use Monte Carlo method to find his most probable distance from the lamppost after 36 steps. Compare it with correct answer, 12. What is the average distance after 100 steps? Also check the result at 1000, and 10000 steps.
6.7 Consider a binomial process in a system with the probability of success equal to 0.1 and probability of failure equal to 0.9. Calculate the variance and standard deviation of number of successes after 100, 1000 and 10000 trials respectively.
6.8 Suppose that you are given with a probability distribution function of the following form in the domain [0, 1]: p(x) a( 1 x^2 ), 0 x 1.
It is zero elsewhere. Find the cumulative probability distribution function. Using uniform random numbers and inverse transform method find the distribution of x.
6.9 Find central difference approximation for a general point indexed as (i, j) to the following partial differential model:
G(x,y) f(x,y) x
G(x,y) x
2
2 2
2
6.10 Develop a normal distributed noise coupled with a sinusoidal function of the following form: f (t) Asinwt Bn(t) , in the range 1 t 10. Then use it in Simulink after developing a block model for first, second and third order filters to study the effect of filters. Compare efficiencies of each filter and perform parametric studies over weighting factors.
6.11 A random variable X has a normal distribution if the probability distribution function is
, x.
(x ) fX (x) exp
2
2 2
only generation from N(0, 1). The inverse transform method cannot be applied to the normal distribution, therefore, suggest some alternative procedure to be employed.
6.12 A two-dimensional picture has normal random noise in it. Design an algorithm that will satisfactorily restore the image and reduce the noise to an acceptable level. Formulate the model and algorithm using MATLAB-Simulink.
6.13 Write down a function program in MATLAB language to test the Hoshen- Kopelman algorithm for cluster counting in a three by three and four by four sized system.
6.14 Simulate for an 8 x 8 lattice based on Ising model with periodic boundary conditions using Metropolis Monte Carlo method. Assume system is at a given temperature T. Plot the specific heat, internal energy and order parameter (M) as a function of T. Assume only nearest neighbor interactions.
6.15 Consider an L x L square lattice with L = 8, 16 and 32 size. Develop a computer program to study the percolation. Compute percolation probability for site percolation as a function of probability (p) of sites thrown on empty lattice. How you will estimate the percolation threshold?
Docsity.com