This paper presents an algorithm of identifying parametric uncertainty by way of an interval model. For a given set of frequency response data from an uncertain linear SISO system of which the upper and the lower bounds of both magnitude and phase responses are represented, the proposed algorithm consists of two main parts: first, the nominal model is identified by using Least Square Estimation (LSE), and then an interval model is constructed by expanding the extremal properties of interval systems, so that tightly enclose the given envelopes within those of interval model. Two numerical examples are given to demonstrate and verify the developed algorithm. The identified interval model can be used for evaluating the worst case performance and stability margins against parametric uncertainty by using some extremal properties on interval systems.