This paper describes the current research status of the λ-tunable Wavelength-Division Multiplexing/Time-Division-Multiplexing-Passive Optical Network(WDM/TDM-PON). This optical access system realizes new cost-effective operations from the network operator′s point of view. The main feature of the proposed system is that it aggressively exploits wavelength tunability. We classify the system configurations and applications. Then, we describe key technologies, namely newly developed λ-tuning burst-mode transceivers and a widely applicable Dynamic Wavelength and Bandwith Allocation(DWBA) algorithm. Finally, we summarize our research and discuss several issues related to the λ-tunable Wavelength-Division Multiplexing/Time-Division-Multiplexing-Passive Optical Network(WDM/TDM-PON).

To solve the problems of complicated structure and rate limiting in conventional format conversion system, 40 Gbps all-optical format conversion from return-to-zero to non-return-to-zero modulation format converter was proposed and demonstrated with a single semiconductor optical amplifier and an optical band-pass filter. The format converter consists of a single semiconductor optical amplifier which is acted as a nonlinear element to broaden the spectrum of input signal and the optical band-pass filter which is used to extract the special spectrum from the broadened spectrum. By adopting the ultra-fast semiconductor optical amplifier model, the 40 Gbps return-to-zero to non-return-to-zero format conversion was demonstrated with simulation. The proof of the principle experiment at 40 Gbps was demonstrated by using the test semiconductor optical amplifier and optical band-pass filter converter. Experimental results show that the power penalty induced by the pattern converter is about 2 dB, which is well coincidence with simulated results. The proposed scheme is robust and potential for applications in future optical networks.

For making a polarization control device with stable wideband, a novel holey birefringence fiber with stable beatlength in broadband was proposed. Some specific asymmetrical structures were introduced into the cladding by adjusting a pair of the special holes near the core, the shape and size of holes in background and the lattice pitches. The polarization beat length is not sensitive to disturbances of these asymmetric structure parameters. The polarization beat length sensitivity to wavelength was reduced effectively, and the polarization beat length sensitivity to the errors of the cladding structures was also reduced through the mutual compensation of birefringence with different trends to suppress non-linear changes. The research results showed that the operating bandwidth was up to 200 nm near the 1 550 nm wavelength window if the relative change of polarization beat length must to be less than ±4%, while the error tolerance of the clad structure geometry parameter was reached 0.1 μm. A more feasible process structure parameters was provided for the actual drawing the holey fiber with broadband stable polarization beat length.

In large scale optical networks, great challenge are brought by the fast growing of multicast services. Multicast lightpath-tree may fail to work in data plane due to the low OSNR (optical signal to noise ratio) value, even though the signaling procedure is successful in the control plane. That is called “false signaling-success” in control plane. To solve this problem, an OSNR Constrained Multicast Lightpath Provisioning Mechanism (OC-MLPM) was proposed. The OC-MLPM makes advantage of the OSNR monitoring function in each node, and is able to allocate the wavelength by judging the OSNR value carried in signaling message in order to make sure lightpath could be established smoothly in both data plane and control plane for multicast service. Simulation results show that the OC-MLPM is able to solve the “false signaling-success” problem, and to reduce the blocking probability effectively as well.

Radar common intra-pulse characteristics are sensitive to signal-noise ratio, thus they are not adaptable to the complex radar systems. To solve this problem, a method was proposed which applied the multi-wavelet packets characteristics of bispectrum to sort unknown complicated radar signals under the large-scale signal-noise ratio condition with a high sorting rate. The bispectrum of received signals was extracted and predigested to two dimensions characteristic. Then, multi-wavelet packets were used to extract characteristics from two dimensions of the bispectrum and the best characteristics were selected as the sorting parameters. The best characteristics of eight classes radar emitter signals were extracted, then kernelized fuzzy c-means was used to cluster and sort signals. The experiment results demonstrate that the characteristics of eight typical radar emitter signals extracted by this method have good performance on noise-resistance and clustering with the large-scale signal-noise ratio.

Aimed at the problems that the existing denoising methods suppress noise incompletely and blur the details of image, an image denoising method using kernel fuzzy C-means clustering and regularization on sparse model was proposed. Firstly, the image was divided into equal pieces and kernel fuzzy C-means clustering algorithm was used for clustering the similar image pieces, thereby ensuring image pieces in the same class share the same sparse denoising model. Then, the global dictionary trained by images from the classical image library was selected as the initial dictionary to adapt to the various characteristics of image very well. Next, a 1/2 norm regularization constraint condition was imposed and sparse decomposition of image pieces in the same class under the dictionary was achieved, which made decomposition coefficients sparser. Finally, the update of dictionary was completed by improved K-singular value decomposition algorithm, and image pieces with the largest difference from the original sparse model were selected to replace the redundancy atoms of the uapdated dictionary. Thus, noise in the image was suppressed effectively. Experimental results show that, compared with denoising method based on wavelet combining with nonlinear diffusion, denoising method based on constant dictionary, denoising method of optimal directions and K-singular value decomposition denoising method, the proposed method can remove noise of the image more effectively and preserve the details of the image and improve the visual effect better.

Based on projective geometry, the relation between the moving camera and still scene was analyzed, and a fast estimation algorithm was proposed based on the optical center coordinate system. On this coordinate in which the camera was always static, the scene and the targets were moving instead. The problem of moving target detection on a moving camera was transformed into distinguishing independent movement of targets from the global motion of scenes on a still platform and a concise form of optical flow constraint was deduced. This framework compensated the random movement according to the optical flow field which was got from KLT. Compared with the original algorithm, the computing speed increased about 10 times while ensuring accuracy and robustness. Experiments were performed based on actual image data, and the results show that this model can reduce the computation time extreamly under the premise of ensuring the accuracy and robustness, and be pratical on moving target detection with a moving camera without construction.

In order to enhance the spatial resolution of multispectral images and avoid severe spectral distortion, the amplified multispectral images and the original panchromatic image were decomposed by non-subsampled contourlet transform with the same layers. A 5×5 sliding window around the fusion pixel of the matrix in the low frequency sub-band of the decomposed images was selected to compute the local correlation coefficient and fourth-order correlation coefficient. If the calculated local correlation coefficient was greater than the value of fourth-order correlation coefficient, the center pixel of the matrix of the high frequency coefficients of multispectral images was replaced by high frequency coefficients of panchromatic image because the position of the ground objects were similar on the spectral characteristics. Otherwise, the high frequency coefficients of multispectral image were kept unchanged. At last, the fused image was reconstructed from the joint of approximation coefficients of multispectral image and detail coefficients obtained from the fusion procedure with inverse non-subsampled contourlet transform. The experimental simulation was conducted on Landsat 7 remote sensing images and the subjective and objective evaluation index of the fused results were given by comparing with the existing similar latest literature. The experimental results indicate that the proposed approach has good effects both on enhancing the original multispectral image and maintaining good spectral characteristics.

A novel interference-based multiple-image encryption method with high encryption capacity was proposed by position and wavelength multiplexing. For encryption, the original images to be encrypted were divided into two groups, and each group of target images was analytically hidden into two phase only masks by position multiplexing. Subsequently the obtained four phase only masks were further merged into two diffractive phase elements by wavelength multiplexing. The encryption process did not need iterative algorithm and should be implemented digitally, while the decryption process could be performed both digitally and optically. The encryption capacity of the system was evaluated with correlation coefficient, and it was shown that the capacity was greatly improved due to the introduction of both position and wavelength multiplexing. In addition, the robustness of the proposal against noise and occlusion was also investigated, and the results indicated that the proposal is of high robust against noise attack but a little vulnerable to occlusion attack.

Strehl ratio is often used as a method for evaluation of the Fourier telescope imaging quality. This method computes the correlation of the reconstructed image and the ideal image to evaluate the image quality, which is mainly applicable in the computer simulation analysis. However, the object to be observed may be unknow in the real world, which means there is no ideal image for use. Based on laboratory imaging data, the feasibility of the non-reference image evaluation method applied to Fourier telescopy is proposed and analyzed. The gray mean grads (GMG) method and Laplacian Sum(LS) are adopted to describe how much detail information the reconstructed image has. For the same object, Employ Strehl ratio and GMG, LS to evaluate the simulated image and the laboratory imaging result, respectively. Statistical results show that the two kinds of evaluations are identical in changing tendency. It proves the GMG and LS methods are applicable in Fourier telescopy as a non-reference image evaluation.

The directional accuracy and wavelength resolution is highly dependent on the laser spot center location accuracy for the super wide-angle laser warning system based on grating diffraction. The spot features of this system was analyzed. As the imaged spots became tilt and its intensity was similar to Gaussian distribution with the laser arrival angle increasing, a new improved Gaussian fitting algorithm was proposed to locate the spot center. Then the accuracy and stability were verified by simulation and experiment test. The simulation results show that when the angle between long axis of imaged spot and x axis increases from 0° to 80° with standard error of Gaussian noise 0.01, mean errors and standard error of this method are smaller than 0.005 pixel and 0.02 pixel, respectively. The experiment test shows that after pre-process of frame subtraction and Gaussian filtering, results of the proposed method are close to weighted centroid and Gaussian fitting if the long axis of imaged spot is parallel to x or y axis, and both mean error and standard error for laser incoming direction of this proposed method are smaller than the other two methods when the imaged spot is tilted.

A novel 2×2 polymer electro-optic switch was optimized and simulated based on two-group five-serial-coupled microring resonator and their cross-coupling with two channel waveguides. Detailed structure, theory and formulation were available to characterize the output powers of the switch. For realizing boxlike spectrum (drop port) as well as low crosstalk and insertion loss, resonance order and coupling gaps were optimized. Analysis and simulation on output powers and output spectrum indicated that, a switching voltage of about 4.0 V was desired to realize the exchange between cross-state and bar-state; the crosstalk under cross state and bar state were extremely low as －66 dB and －54.7 dB, respectively; the insertion losses under the two operation states were 2.34 dB and 0.24 dB, respectively. For the drop port, the 10%~90% rise time and fall time were estimated to be 15 ps and 90 ps, respectively, under the operation of 1 GHz switching operation.The bending radius of each electro-optic polymer microring was as small as 19.45 μm, leading to ultra-compact size of only 0.407 mm in both length and width, which was nearly 1/10 of the length of Mach-Zehnder interferometer or directional coupler type polymer electro-optic switches. Consequently, because of small footprint size and extremely low crosstalk, this switching configuration can be densely integrated onto optoelectronic chips, and thus it possesses potential applications in optical signal control in optical networks-on-chip.

Two types of photodetectors including conventional P+/N-EPI/BN+ photodetector and multi-finger P+/N-EPI/BN+ photodetector were implemented in a standard 0.5 μm Biplor, CMOS and DMOS process with different sizes. The performance and parameter of photodetectors were simulated and optimized. Results of the characteristics of photodetectors simulations and tests were presented, which showed that the multi-finger structure P+/N-EPI/BN+ photodetector can improve the responsibility at 650 nm and decrease the junction capacitance. The large-area multi-finger structure P+/N-EPI/BN+ photodetector was chosen for the monolithic integration with a trans-impedance amplifier and a preamplifier. A proposed monolithic optoelectronic integrated receiver was fabricated in 0.5 μm Biplor, CMOS and DMOS technology for 650 nm plastic optical fiber communication. The measurements of the proposed receiver were done. The receiver achieved a sensitivity of -15 dBm with the bit-error-rate of 10－9 at 160 Mb/s pseudo random binary sequence signal for 650 nm input light. A clear eye diagram was demonstrated for 160 Mb/s pseudo random binary sequence signal. These indicate that the proposed photodetector can be employed in the receiver chip of high-speed plastic-optical-fiber-based fast ethernet system for broadband access network application.

Based on the method of the circuit, N-order series-coupled microring resonator filters was equivalent to a baseband LC ladder network. By solving the element parameters of LC ladder network, the coupling coefficients of microring filters can be obtained. Using this method, microring filters with Chebyshev response were designed. To prove the reliability of this method, it was compared with the traditional coupled mode theory method and the ideal Chebyshev response. Finally, the effects of coupling coefficients and the ring number on the output characteristics of Chebyshev microring filters were discussed. Numerical simulation indicated that the passband flatness was worsened with the deviation of coupling coefficients increasing, and at the same time, the bandwidth of the passband would increase with the coupling coefficients. And a square passband shape was achieved for the larger ring number and the faster transition attenuation.

The visible micro-imaging technology cannot observe the chip internal structure. In order to sovle this issue, a silicon-based semiconductor chip laser infrared microscopy imaging system was designed. The characteristic of 1 064 nm wavelength of infrared laser that has a certain penetration depth for silicon material was applied for the design. A microscope objective with numerical aperture of 0.42 was adopted in this system. The multimode silica fiber vibration method was used to eliminate the speckle noise. The track pitch of CD-RW disk image was observed by microscope infrared imaging. The results showed that the system resolution can reach 1.6 μm and close to the theoretical value. So it can achieve internal structure observation of static RAM of which the thickness of the chip is 70 μm.

Mid-infrared lasers have a lot of potential applications, such as laser medicine, spectroscopy and infrared countermeasures. In order to obtain a compact and potable mid-infrared laser source, a diode-pumped fiber laser emitting at 2.8 μm from a heavily erbium-doped ZrF4-BaF2-LaF3-AlF3-NaF double-clad fiber was demonstrated. A fiber-coupled laser diode centered at 975 nm was coupled into the inner cladding as the pumping source by an achromatic coupling system consisted of two aspheric lenses. The laser cavity was composed of a protected gold mirror butted against the rear end of the active fiber, and the other end with 4% Fresnel reflection. The 2.8 μm lasing was achieved when the launched pump power was higher than 0.37 W. The experimental results indicated that the center wavelength and linewidth of free running were 2.785 μm and 0.9 nm, respectively. Moreover, the maximum output power of the fiber laser working at the transverse-fundamental-mode was 0.98W, corresponding to the slop efficiency of 17%. High efficient watt-level single mode mid-infrared fiber laser was obtained using the efficient energy transfer upconversion processes between Er-Er ions.

Luminescence from Bi and Bi/Er doped aluminum calcium germanate glasses were investigated. The glasses were prepared by conventional melting-quenching technique. The near infrared luminescence of Bi-doped calcium aluminum germinate glasses can be divided into two Gaussian peaks which are located at 1 265 and 1 420 nm separately. The X-ray photoelectron spectroscopy spectrum indicates the valence state of Bi ions in the calcium aluminum germinate glasses is the mixed. Bi3+ ions could be reduced to low valence states bismuth ions. The experiments support the assumption that the near-infrared luminescence is connected with a low valence state of Bi, may be Bi+. Er3+ was introduced to adjust the 1420 nm band of Bi-related fluorescence. With the increment of Er3+ concentration, near-infrared fluorescence relative intensity of Bi ions is decreasing, which implies the energy transfer from Bi ions to Er3+.

A series of orange-red phosphors Ca3Y2－2x(Si3O9)2∶2xSm3+ were synthesized by conventional solid-state reaction, and their photoluminescence properties were investigated. The X-ray diffraction (XRD) pattern shows that Ca3Y2(Si3O9)2∶Sm3+ crystals are the pure phase. The excitation spectra contain the characteristic excitation of Sm3+. The emission spectra of Ca3Y2(Si3O9)2∶Sm3+ phosphors exhibit three main peaks assigned to the 4G5/2→6HJ/2(J=5,7,9) transitions of Sm3+ under 402 nm excited radiation, the dominating emission peaks at 565 nm, 604 nm, 651 nm. The decay time of 4G5/2 level in Sm3+ was measured by the time resolved spectrum. The luminescence intensity firstly increases with increasing of Sm3+ concentration, and then decreases, the emission reaches the maximum intensity at x=0.02, and the concentration quenching mechanism is the electric dipole-electric dipole interaction.

Exfoliated graphite is a potential military passive interfering material, which can attenuate infrared and millimeter-wave simultaneously. Exfoliated graphites with different expanding volumes were prepared by two-step intercalation, in which the mixture of nitric acid and phosphoric acid, nitric acid and acetic acid were employed as intercalating agent, and potassium permanganate was used as oxidant. The microstructures of exfoliated graphites with different expanding volumes were analyzed by scanning electron microscope. The millimeter-wave attenuation performances of exfoliated graphites with different expanding volumes were measured by a static method. The results show that the exfoliated graphites′ sheets become thinner with the increasing of expanding volume, the interlayers are fully opened and there are more tiny graphite flakes with size distribution between several microns and several ten microns. And the millimeter-wave attenuation performance enhances with the increasing of expanding volume. When the expanding volume increases from 233 mL/g to 450 mL/g, the 3 mm wave attenuation value rises from 5.7 dB to 8.68 dB.

Space combining beam of Terahertz was designed in the paper. The optical intensity of single THz beam, double THz beams (2×1) and four THz beams (2×2) were calculated in the condition of spatially combining images. Through the obtained images, the relationship between the receiving distance, the Gaussian beam peak separation distance and combining intensity were obtained, and the effect of phase difference on combining beam was analyzed. The results show that THz peak intensity is enhanced by near four and sixteen multiples in the conditions of 2×1 and 2×2 arrangements. The optimal THz beam receiving distances of 2×1 and 2×2 arrangements are found through calculation and images, and they are 0.325 m and 0.223 m. Gaussian beam peak separation distance and peak intensity are nonlinearly inverse relationship, which has no maximum or minimum points. The image of 2×1 arrangement was achieved when there is phase difference between two THz beams, and the minimum point of THz intensity of π was attained in the image. So it can be concluded that the phase difference can be decreased but it can not be avoided.

Steady-state spatial soliton pairs resulted from the combination of pyroelectric and photovoltaic effect were predicted for unbiased photovoltaic photorefractive crystals. The numerical solutions of pyroelectric photovoltaic soliton pairs were deduced theoretically based on the incoherent coupling wave equations. The intensity envelopes of soliton pairs were obtained by finite difference method. The results show that the bright-bright soliton pair can be established in self-defocusing photovoltaic crystals when the positive pyroelectric field is large enough compared with photovoltaic field. For the dark-dark and grey-grey soliton pairs, the negative pyroelectric field can enhance the self-defocusing nonlinear effect and reduce the full width at half maximum of soliton pairs. The full width at half maximum of soliton pairs can be controlled by adjusting the value of the pyroelectric field. Moreover, the pyroelectric photovoltaic soliton pairs can degenerate into pyroelectric soliton pairs or photovoltaic soliton pairs under a certain conditions. The characteristics of photorefractive crystal LiNbO3 can transfer from self-defocusing to self-focusing by virtue of the pyroelectric effects, which can help to control the characteristics of soliton pairs.

The size of a single CCD cannot meet the requirements of modern optical system. To enlarge the field of optical system, a reflector based butting system was adopted. Aiming at the problems of vignetting and decline of modulation transfer function caused by butting, a reflector based butting system which has 6 reflectors was investigated. By analyzing the feature of vignetting, a sub-image pixel overlapping method was proposed to reduce the vignetting. At the same time, vignetting at edge was restored, and the modulation transfer function was improved by compositive point spread function. The results show that though the cut-off frequency has not been increased, the value of modulation transfer function is improved obviously.

Based on the thermal control system of Space solar telescope, the temperature of primary mirror in deep-space solar observatory was simulated under new orbit and larger view. Based on the theory of thermal-structural-optical integrated analysis, thermal deformation of the primary mirror under supported condition was calculated. And the deformed surface point of the primary mirror is fitted using parabolic polynomial and Zernike formula. The Root Mean Square(RMS) of deformed surface was greater than λ/40(λ=533 nm). It indicates the former thermal control system cannot meet its 0.1″ diffraction limit quality. Owing to the former thermal control system meets the demands of optical design difficultly, the supporting system of mirror is suggested to be adjusted in order to decrease thermal deformation. The integrated analysis of primary mirror is useful for optimizing the thermal control system and supporting system of Deep-space solar observatory.

A new descent camera which used for deep space explorartion was designed. This camera can explore different distance targets by opening and closing the first optical lens. With the first optical lens in the optical path, the objective was focused for objects at distances ranging from 1.5 m to infinity. In this mode, the camera can take full-color pictures (spectral bandwidth 420~900 nm) of the cytaster surface. With it removed, the camera was focused at an object distance of 30 cm. In this mode, three lights of different wavelengths which were blue-green, red, and NIR were used independently to illuminate the close range objects and the camera can take multispectral pictures of the explored objects. This camera have a focal ratio f/5.4 and project a square field of view 49.5°×49.5° on the CCD. When at the infinity mode (with the first optical lens), the focal length is 12.9 mm and when at the close-up mode (without the first optical lens), the focal lenth is 12.5 mm. The first optical lens of the descent camera is a meniscus lens and it almost have no focal power. So the system focal changes lightly no matter with or without this lens. By this way, the drawback of a traditional two-way switch system that detected objects easily be lost when modes shift can be easily overcome. The design result show that the Modulation Transfer Function is higher than 0.6 in full-field at 35 lp/mm. The distortion of the system in full-field is better than 1.2%. Tilt and decenter tolerance of the first lens also have a reasonable result.

It is difficult to fabricate and test convex aspherical surfaces, especially for steep convex aspheical surfaces. Based on the analysis of the conventional null test, a simple test system with cemented doublet for two steep convex ellipsoids was discussed. The wavefront of the both test system was less than 1/30λ RMS (Root Meat Square, λ=632.8 nm). The convex aspheric lens had their successful application in practical optical system. The result showed that the null test for convex aspheric lens discussed in this paper is valid. Compared with the traditional test method, the cemented doublet offsets the spherical aberration generated by the sphere surface of the lens. So the proposed test method decreases the amount of the spherical aberration generated by the aspherical lens, with fewer elements, and reduces the system′s alignment difficulty. It was demonstrated to be an effective, high-precision testing method.

The topological transformation of vortex beams using cylindrical lens is experimentally and numerically demonstrated. The intensity distributions and phase structures of vortex beams with different topological charges were analyzed after the transformation by a cylindrical lens, which were generated using the computing generated holographic method. The results show that the transformation of the cylindrical lens results in the sign change of the topological charges of the vortex beams. In addition, the higher order vortex beams decay into a beam with multiple singularities with charge-one. These observations were analyzed based on the linear superposition of higher order laser modes and Gouy phase delay. Finally, the experiment results were verified by the numerical simulations.

Good security is indispensable to any quantum secret communication protocol. However, information leakage has been a great security threat to quantum dialogue. In order to solve the problem, a quantum dialogue protocol without information leakage via cavity QED was proposed, which made full use of the evolution law of atoms in cavity QED. The proposed quantum dialogue protocol avoided the information leakage problem by using the measurement correlation property after entanglement swapping between two Bell states via cavity QED, and could securely exchange 4 bits secret messages per round communication. The results show that the proposed protocol is able to detect the active attacks from the outside eavesdropper through security checking, such as the intercept-resend attack, the measurement-resend attack and the entanglement-and-measurement attack. Therefore, it has good security towards both the information leakage problem and the active attacks from the outside eavesdropper.

The atom fidelity and degree of entanglement were calculated in the system consisting of M two-level atoms and M multi-mode fields by use of complete quantum theory and number computing method. Through number computing and analyzing, it was found that the atom fidelity and the atom entangled degree present always oscillatory with opposite in phase if initial atoms are in separate state, and both of their frequency gradually increase with the intensity of light field increasing. The atom fidelity tends to 0 and degree of entanglement tends to 1 under strongly coupling between atoms and light fields. The small fidelity is disadvantage of information transfer but high-fidelity atom information transfer can still obtained by controlling appropriately interaction time.

Performance of optical systems is closely related to the scattering characteristics of optical mirror surface. According to Mie theory, the bidirectional reflectance distribution function of optical mirror surface under four cases of mirror surface cleanliness levels such as clean, basic clean (cleanliness level 200), slight pollution (cleanliness level 500) and severe pollution (cleanliness level 750) was simulated and computed by MATLAB software when incident laser wavelengths are the visible light (0.632 8 μm), near infrared light (1.053 μm) and thermal infrared light (10.6 μm), respectively. The results show that the bidirectional reflectance distribution function of mirror surface is affected by mirror surface roughness and clean lines level, which will make a very important impact on properties of optical system. If wavelength is shorter, bidirectional reflectance distribution function will be greater, and scattering becomes more complicated. Useful help will be provided to estimate the pollution of optical mirror and clean the mirror surface in optical engineering.

On the basis of the method of moments and segments sinusoidal basis functions, the impedance matrix of plasmonic nano-rod antennas can be obtained by solving the Green′s function of the discretized current nodes in closed form formulations. Compared to the method of moments with different basis functions, this method can be utilized for reducing the size of matrix equations. The simulated results show that the polarization current and scattered fields can be computed accurately by employing a small-size matrix, which leads to a fast analysis for the scattering characteristic and resonant modes of nano-rod. The accuracy of the method was successfully demonstrated by comparing the simulated results with finite difference time domain technique, which is more time-consuming than the present method. Especially when dealing with the oblique incidence problem, the advantages become more apparently. The method is effective for the calculated models in this paper, and provides a fast and efficient electromagnetic analysis way for other nano-rod antennas and nanotube devices in different shapes when doing their scattering characteristics simulation.