 α= 11 deg and impact angle β= -5.4 deg; which incidentally reproduces the middle part of the PA sweep at 610 MHz) and G-2 (α=175.4 deg and β=- 6.9 deg; geometry derived from RVM fit to 325 MHz PA sweep). Pulsar radiation pattern simulated with both the geometries reproduces the average profile as well as the observed features in the drift pattern quite well. However, G-2 fits the PA sweep much better. We report that the peaks of the emission from the trailing and leading outer regions, as a function of the pulse number, are offset by a constant interval, P5~9P1. We also report a phase locked relation (PLR) between the inner and outer drift regions for PSR B0818-41. A new technique is introduced by us for resolving aliasing, using this constant offset (P5 ~ 9P1) between the peak emission from the leading and trailing outer regions. From the result of this technique, we propose that the subpulse drifting for PSR B0818 -41 is most likely first order aliased, and the corresponding carousel rotation period 4 =10 s. This implies that PSR B0818 -41 has the fastest known carousel. The drift pattern in the inner and outer rings are always phase locked for PSR B0818 -41. This could be a signiﬁcant constraint for the theoretical models of pulsar radio emission, and favors a pan magnetospheric emission mechanism. We observe frequent nulling for PSR B0818-41. We calculate a nulling fraction ~30% at 325 MHz for this pulsar. Lengths of neighboring nulls and bursts are found to be independent. For the inner drift region, our investigations bring out the fact that the nature of the transitions from burst to null are different from the transitions from null to burst. Switching off of pulsar radiation during nulling for PSR B0818 -41 is not abrupt, but is gradual, whereas the transitions from null to burst are found to be rather abrupt for the inner drift region. This effect is not prominent in the outer drift regions. Although, the inner region of the last active pulses before nulls are dimmer, the ﬁrst active pulses after nulls outshines the normal ones. The intensity of the inner region is maximum for the average proﬁle from the ﬁrst active pulse immediately after the nulls and then gradually goes down. This is consistent with the behavior of the individual nulls described above. However, this is not the case for the leading and trailing outer regions. The average profiles from the first active pulse immediately after the nulls follows similar shape as the normal profile but shows an increased intensity (in the form of a bump) in the inner region which is not present in the normal average profile. In addition, the leading and the trailing peaks appear to be of similar intensity, while trailing peak is significantly more intense for the normal profile. The average profiles from the pulses immediately after the nulls are wider than the normal profile. The average profiles of the first active pulses after the nulls are drastically similar between two epochs of observations. This is a very unique result which is not reported for any other pulsar so far and may imply that the phenomenon of nulling is associated with some systematic energy re-distribution in the pulsar magnetosphere. In depth study of PSR B0826-34 PSR B0826-34 is a pulsar with one of the widest known profile. The earlier studies of this pulsar (Durdin et al. (1979), Biggs et al. (1985) and Gupta et al. (2004)) have brought out some unique properties : strong evolution of the average profile with frequency, apparent nulling for 70% of time and a remarkable subpulse drift property- multiple curved drift bands with frequent changes and sign reversals of drift rate. We studied PSR B0826 -34 using the GMRT, simultaneously at 303 and 610 MHz, and individually at 157, 325, 610 and 1060 MHz. Detailed investigation of PSR B0826- 34 are reported in Chapter 5 of this thesis. Some of the interesting new results from our work are, As a natural out-come of the simultaneous dual frequency observations, we obtain an accurate DM value, equal to 52.2(6) pc/cm3, for this pulsar. Unlike most no