合肥网站维护,公司建立网站的步骤,网站建设费是什么费用,刷赞网站推广qq关注作者了解更多
我的其他CSDN专栏
过程控制系统
工程测试技术
虚拟仪器技术
可编程控制器
工业现场总线
数字图像处理
智能控制
传感器技术
嵌入式系统
复变函数与积分变换
单片机原理
线性代数
大学物理
热工与工程流体力学
数字信号处理
光电融合集成电路… 关注作者了解更多
我的其他CSDN专栏
过程控制系统
工程测试技术
虚拟仪器技术
可编程控制器
工业现场总线
数字图像处理
智能控制
传感器技术
嵌入式系统
复变函数与积分变换
单片机原理
线性代数
大学物理
热工与工程流体力学
数字信号处理
光电融合集成电路技术
电路原理
模拟电子技术
高等数学
概率论与数理统计
数据结构
C语言
模式识别原理
自动控制原理
数字电子技术
关注作者了解更多
资料来源于网络如有侵权请联系编者
目录
实验一 离散时间信号与系统的时域分析
一、实验目的
二、实验原理
三、实验过程记录
四、结果分析 实验二 离散信号 Z 变换
一、实验目的
二、实验原理
三、离散系统的频率特性
实验三 IIR 滤波器的设计与信号滤波
、实验目的
2.实验原理
、实验步骤及内容
实验代码
6实验报告要求
实验四 用窗函数法设计 FIR 数字滤波器
、实验目的
、实验原理
、实验步骤及内容
6实验要求
实验五 用 FFT 作谱分析
1实验目的
、实验原理
、实验步骤及内容
实验数据记录
6实验报告要求
实验 6 综合实验
1实验目的
、实验内容 三选一
3实验报告要求 实验一 离散时间信号与系统的时域分析 1
实验二 离散信号 Z 变换 7
实验三 双线性变换法设计 IIR 数字滤波器 19
实验四 用窗函数法设计 FIR 数字滤波器 19
实验五 用 FFT 作谱分析 19
实验六 综合实验 25 实验一 离散时间信号与系统的时域分析 一、实验目的 1掌握离散时间信号与系统的时域分析方法。 2掌握序列傅氏变换的计算机实现方法利用序列的傅氏变换对离散信号、系统及系 统响应进行频域分析。 3熟悉理想采样的性质了解信号采样前后的频谱变化加深对采样定理的理解。 二、实验原理 实现2离散时间系统的单位脉冲响应3连续时间信号的采样绘制单位脉冲序列绘制单位阶跃序列绘制理想采样序列绘制正弦序列分析特性 三、实验过程记录 1单位脉冲序列 clear all;
clc;
n1-5;n25;n00;
nn1:n2;%自变量
x[nn0];%因变量利用大括号判断里面的表达式是否成立如果成立则函数有效有数值
%否侧函数无效没有数值
stem(n,x,filled);%利用filled绘制实心火柴梗
axis([n1,n2,0,1.1*max(x)]);%绘制坐标系
xlabel((n));ylabel(幅度x(n));
title(单位脉冲序列);
2单位冲激序列 clear all
n 40;
x ones(1,n);
xn 0:n-1;
subplot(211);
stem(xn,x);
grid on
axis([-1 51 0 1.1]);
title(单位阶跃序列u(n))
ylabel(u(n));
xlabel(n);
x [zeros(1,10),1,ones(1,29)];
subplot(212);
stem(xn,x);
grid on
axis([-1 51 0 1.1]);
title(单位跃迁序列u(n-10))
ylabel(u(n-10));
xlabel(n); 3正弦序列
clear all
n 0:59;
x sin(pi/7*n);
stem(n,x);
ylabel(x(n));
xlabel(n);
title(正弦序列)
axis([0 40 -1.5 1.5]);
grid on; 4理想采样序列
%绘制理想采样信号序列
%给定信号参数
A 200;
alpha 50*pi;
T 1/10000;
Omega 260*pi;
n 0:1:100;
%计算信号在不同时刻的值
x A.*exp(-alpha.*n.*T).*cos(Omega.*T.*n);
figure %绘制理想采样信号序列
stem(x,g) %绘制单位抽样序列
%修改时间坐标值保证信号0时刻的值为1添加图的坐标含义标题等信息
set(gca, XTick, 1:101);
set(gca,XTickLabel,0:100);
title((d) Desired sampled signal)
grid on
xlabel(Time index n)
ylabel(Amplitude)
axis([1 101 -200 200]) %规范坐标轴范围
hold on
%给定信号参数
A 100;
alpha 0.4;
T 1/1000;
Omega 700*pi;
n 0:100;
%计算信号在不同时刻的值
xx A.*exp(-alpha.*n.*T).*cos(Omega.*T.*n);
%绘制修改参数后的理想采样信号序列
stem(xx,b) %绘制单位抽样序列 改变参数A1a 0.4, W 700p 产生理想采样信号序列 x n( ) 。
代码
%给定信号参数
A 1;
alpha 0.4;
T 1/10000;
Omega 700*pi;
n 0:1:100;
%计算信号在不同时刻的值
x A.*exp(-alpha.*n.*T).*cos(Omega.*T.*n);
figure %绘制理想采样信号序列
stem(x,g) %绘制单位抽样序列
%修改时间坐标值保证信号0时刻的值为1添加图的坐标含义标题等信息
set(gca, XTick, 1:101);
set(gca,XTickLabel,0:100);
title((d) Desired sampled signal)
grid on
xlabel(Time index n)
ylabel(Amplitude)
axis([1 101 -200 200]) %规范坐标轴范围
hold on
%给定信号参数
A 100;
alpha 0.4;
T 1/1000;
Omega 700*pi;
n 0:100;
%计算信号在不同时刻的值
xx A.*exp(-alpha.*n.*T).*cos(Omega.*T.*n);
%绘制修改参数后的理想采样信号序列
stem(xx,b) %绘制单位抽样序列
结果截图 5.实指数序列 clear all
n 0:20;
x1 1.3.^n;
x2 (-1.3).^n;
x3 0.8.^n;
x4 (-0.8).^n;
subplot(221)
stem(n,x1,fill)
grid on;
ylabel(x1(n));
xlabel(n);
title(x1 1.6^n)
subplot(222)
stem(n,x2,fill)
grid on;
ylabel(x2(n));
xlabel(n);
title(x2 (-1.6)^n)
subplot(223)
stem(n,x3,fill)
grid on;
ylabel(x3(n));
xlabel(n);
title(x1 0.8^n)
subplot(224)
stem(n,x4,fill)
grid on;
ylabel(x4(n));
xlabel(n);
title(x1 1-0.8^n) 四、结果分析
分别绘制了单位脉冲序列单位阶跃序列正弦序列理想采样的指数正弦序列改变参数后的理想采样正弦序列实指数序列等图像 实验二 离散信号 Z 变换 一、实验目的
1、熟悉离散信号 Z 变换的原理及性质
2、熟悉常见信号的 Z 变换 3、了解正/反 Z 变换的 MATLAB 实现方法
4、了解离散信号的 Z 变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系
5、了解利用 MATLAB 实现离散系统的频率特性分析的方法
二、实验原理 实验过程记录
一、正/反 Z 变换 代码 clc;
clear;
close all;
%%
% Example 1 for the Z transform
syms n % Define the temporal independent variable represented by the mathmatical symbols.
f 0.5^n; % Define the sequence in the time domain.
Fz ztrans(f); % Perform z transform based on the function ztrans.
disp(Fz); % Print the output.
%%
% Example 2 for the Z transform
syms z % Define the independent variable in the complex frequency domain.
Fz 2*z/(2*z - 1); % Define the z transform in the complex frequency domain.
f iztrans(Fz); % Perform inverse z transform based on the function iztrans.
disp(f); % Print the output
%%
% Example 3 for the Z transform
syms n % Define the independent variable
hn kroneckerDelta(n, 1) kroneckerDelta(n, 2) kroneckerDelta(n,3); %Define the sequence in the time domain
disp(hn)
Hz ztrans(hn); % Perform the z transform
disp(Hz)
Hz simplify(Hz); % Simplify the expression of the H(z)
disp(Hz) 运行结果 z/(z - 1/2) (1/2)^n kroneckerDelta(n - 1, 0) kroneckerDelta(n - 2, 0) kroneckerDelta(n - 3, 0) 1/z 1/z^2 1/z^3 (z^2 z 1)/z^3 三、离散系统的频率特性 代码 %Example 4
a [1,0]; % Define the coefficients of the denominator polynomial
b [1,-0.5]; % Define the coefficients of the numerator polynomial
[H,omega] freqz(b,a,10); % Compute the frequency response of the transfer function
disp(H)
disp(omega)
%%
%Example 5
a [1,0]; % Define the coefficients of the denominator polynomial
b [1,-0.5]; % Define the coefficients of the numerator polynomial
[H,omega] freqz(b,a,200,whole); % Compute the frequency response of the transfer function
disp(H)
disp(omega)
Ha abs(H); % Amplitude
Ph angle(H); % Phase
运行结果 求点代码如上
绘图代码
figure
subplot(2,1,1)
plot(omega,Ha)
hold on
set(gca, XTick, 0:pi/10:2*pi);
set(gca,XTickLabel,0:pi/10:2*pi);
title((a) Amplitude-frequency response)
grid on
xlabel(omega)
ylabel(Amplitude)
axis([0 2*pi 0.5 1.5])
grid on
subplot(2,1,2)
plot(omega,Ph)
set(gca, XTick, 0:pi/10:2*pi);
set(gca,XTickLabel,0:pi/10:2*pi);
title((a) Phase-frequency response)
grid on
xlabel(omega)
ylabel(Phase)
axis([0 2*pi -0.6 0.6])
grid on 运行结果 三、实验内容 代码
syms z; % Define the independent variable
F1z (z^2z1)/(z^2z-2);
f1 iztrans(F1z); % Perform the inverse z transform
disp(f1)
运行结果
(-2)^n/2 - kroneckerDelta(n, 0)/2 1 题目二
代码
syms z;Flz(z^4z^3z^2z1)/(z^4);fliztrans(Flz);disp(fl)
结果
kroneckerDelta(n - 1, 0) kroneckerDelta(n - 2, 0) kroneckerDelta(n - 3, 0) kroneckerDelta(n - 4, 0) kroneckerDelta(n, 0) 题目三
代码
syms z;Flz(2*z^26*z12)/(z^4);fliztrans(Flz);disp(fl)
结果
2*kroneckerDelta(n - 2, 0) 6*kroneckerDelta(n - 3, 0) 12*kroneckerDelta(n - 4, 0)5.实验要求 六实验感想
通过Z变换和Z反变换熟悉了MATLAB的使用 实验三 IIR 滤波器的设计与信号滤波
、实验目的 熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法。
掌握数字滤波器的计算机仿真方法。
通过观察对实际心电图信号的滤波作用获得数字滤波的感性知识
2.实验原理 、实验步骤及内容 实验代码
用双线性变换法设计Butterworth低通数字滤波器要求通带频率为 通带波纹小于1dB阻带频率为 幅度衰减大于15dB采样频率为Fs100Hz。
clc;
%写出已知条件设置采样点数
wp0.25*pi;ws0.4*pi;
rp1;as15;
fs2000;Nn128;
%由数字指标映射到模拟滤波器指标
WPwp*fs;WSws*fs;
%利用巴特沃斯滤波器的设计方法设计一个模拟滤波器
[N,Wn]buttord(WP,WS,rp,as,s);
[z,p,k]buttap(N);
[B,A]zp2tf(z,p,k);
%转换到数字滤波器
[b,a]lp2lp(B,A,Wn);%从一个低通原型变换为具有不同截止频率的低通滤波器
[bz,az]impinvar(b,a,fs);%根据冲激响应不变法得到数字滤波器
freqz(bz,az,Nn,fs);
WS2*fs*tan(ws/2);
WP2*fs*tan(wp/2);
[bz,az]bilinear(b,a,fs);
freqz(bz,az,Nn,fs);[H,w]freqz(bz,az);
figure;
subplot(211);
plot(w/pi,abs(H));
subplot(212);
plot(w/pi,angle(H));
运行结果 6实验报告要求 7实验感想 通过IIR滤波器的设计熟悉了MATLAB的使用 实验四 用窗函数法设计 FIR 数字滤波器
、实验目的 掌握用窗函数法设计 FIR 数字滤波器的原理和方法。
熟悉线性相位 FIR 数字滤波器特性。
了解各种窗函数对滤波特性的影响
、实验原理
、 、实验步骤及内容 实验代码
1设计滤波器采用汉明窗长度41阶数40采样频率2000hz
clear all; close all; clc;
% 滤波器长度
N41;
%采样频率
fs2000;%各种滤波器的特征频率
fc_lpf200;
fc_hpf200;
fp_bandpass[200 400];
fc_stop[200 400];%以采样频率的一般对频率归一化
wn_lpffc_lpf*2/fs;
wn_hpffc_hpf*2/fs;
wn_bandpassfp_bandpass*2/fs;
wn_stopfc_stop*2/fs;%采用fir1函数设计FIR滤波器
b_lpffir1(N-1,wn_lpf);
b_hpffir1(N-1,wn_hpf,high);
b_bandpassfir1(N-1,wn_bandpass,bandpass);
b_stopfir1(N-1,wn_stop,stop);%求幅频响应
m_lpf20*log(abs(fft(b_lpf)))/log(10);
m_hpf20*log(abs(fft(b_hpf)))/log(10);
m_bandpass20*log(abs(fft(b_bandpass)))/log(10);
m_stop20*log(abs(fft(b_stop)))/log(10);% 设置频率响应的横坐标单位为hz
x_f0:(fs/length(m_lpf)):fs/2;% 单位脉冲响应
subplot(4,2,1);stem(b_lpf);xlabel(n);ylabel(h(n));legend(lpf);
subplot(4,2,3);stem(b_hpf);xlabel(n);ylabel(h(n));legend(hpf);
subplot(4,2,5);stem(b_bandpass);xlabel(n);ylabel(h(n));legend(bandpass);
subplot(4,2,7);stem(b_stop);xlabel(n);ylabel(h(n));legend(stop);% 幅频响应
subplot(4,2,2);plot(x_f,m_lpf(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(lpf)
subplot(4,2,4);plot(x_f,m_hpf(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(hpf)
subplot(4,2,6);plot(x_f,m_bandpass(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(bandpass)
subplot(4,2,8);plot(x_f,m_stop(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(stop);
仿真脉冲响应和幅频响应图 6实验要求 7.实验感想
通过FIR滤波器设计熟悉了MATLAB的使用 实验五 用 FFT 作谱分析
1实验目的 进一步加深 DFT 算法原理和基本性质的理解因为 FFT 只是 DFT 的一种快速算法所 以 FFT 的运算结果必然满足 DFT 的性质
熟悉 FFT 算法原理及子程序的应用。
掌握用 FFT 对连续信号和时域离散信号进行频谱分析的基本方法。了解可能出现的分 析误差和原因以便在实际中正确应用 FFT。
、实验原理 、实验步骤及内容 实验数据记录
%《数字信号处理(第4版》第10章实验2程序exp2a.m
% 西安电子科技大学出版社出版 高西全 丁玉美 合著 2016年
% 时域采样理论验证
Tp64/1000;
%产生M长采样序列x(n)
% Fs1000;T1/Fs;
Fs1000;T1/Fs;
MTp*Fs;n0:M-1;
A444.128;alphpi*50*2^0.5;omegapi*50*2^0.5;
xntA*exp(-alph*n*T).*sin(omega*n*T);
XkT*fft(xnt,M); %M点FFT[xnt)]
ynxa(nT);
subplot(3,2,1);tstem(xnt,yn);box on;title((a) Fs1000Hz);
k0:M-1;fkk/Tp;
subplot(3,2,2);plot(fk,abs(Xk));title((a) T*FT[xa(nT)],Fs1000Hz);
xlabel(f(Hz));
ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))])
% Fs300;T1/Fs;
Fs300;T1/Fs;
Mceil(Tp*Fs);n0:M-1;
xntA*exp(-alph*n*T).*sin(omega*n*T);
XkT*fft(xnt,M); %M点FFT[xnt)]
ynxa(nT);
subplot(3,2,3);tstem(xnt,yn);box on;title((b) Fs300Hz);
k0:M-1;fkk/Tp;
subplot(3,2,4);plot(fk,abs(Xk));title((b) T*FT[xa(nT),Fs300Hz);
xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))])
% Fs200;T1/Fs;
Fs200;T1/Fs;
Mceil(Tp*Fs);n0:M-1;
xntA*exp(-alph*n*T).*sin(omega*n*T);
XkT*fft(xnt,M); %M点FFT[xnt)]
ynxa(nT);subplot(3,2,5);tstem(xnt,yn);box on;title((c) Fs200Hz);
k0:M-1;fkk/Tp;
subplot(3,2,6);plot(fk,abs(Xk));title((c) T*FT[xa(nT),Fs200Hz);
xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))])先对下列序列进行谱分析 选择FFT的变换区间N为8和16两种情况进行频谱分析。分别打印其幅频特性曲线并进行对比、分析和讨论。
对于第一个序列 对于第2个序列 对于第3个序列 对下列周期序列进行频谱分析 选择FFT的变换区间N为8和16两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线并进行对比、分析和讨论。
对于第4个序列 对于第5个序列 对模拟周期信号进行频谱分析 选择采样率Fs64Hz对变换区间N163264三种情况进行谱分析。分别打印其幅频特性并进行分析和讨论。
对于第6个序列 6实验报告要求 内容总结
图1a和1b说明x1(n)R4(n)的8点DFT和16点DFT分别是x1(n)的频谱函数的8点采样和16点采样图2a和图3a表明x3(n)x2(n3)8R8(n)的8点DFT的模相等但是当N16时x3(n)与x2(n)不满足位移关系所以导致2b和3b的模不相同。 x4(n)的周期为8所以N8和N16是其周期的整数倍得到正确的单一频率正弦波的频谱仅在0.25Π处有1跟单一谱线如图4a和4b所示。 X5(n)的周期为16所以N8时得到的频谱不正确如5a所示。当N16时为其周期的整数倍故能得到正确的频谱在0.25Π和0.125Π处有两根单一谱线如图5b所示。 x6(n)有三个频率成分f14hzf28hzf310hz。所以x6(n)的周期为0.5采样频率Fs64Hz 16f1 8f2 6.4f3。变换区间N16时观察时间Tp 16T 0.25s,不是x6(t)的整数倍周期所以所得频谱不正确如图6a所示。变换区间N32,64时观察时间Tp0.5s1s是x6(t)的整数周期所以所得频谱正确如图6b和6c所示。图中3根谱线正好位于4Hz,8Hz,10Hz处。变换区间N64时频谱幅度是变换区间N32时2倍. 通过以上实验进一步论证了当周期信号的频谱是离散谱只有用整数倍周期的长度做FFT得到的离散谱才能代表周期信号的频谱。 实验 6 综合实验 1实验目的
1了解语音信号的获取及基本处理
2给定设计指标能选用多种方法设计出符合指标要求的滤波器。
3掌握各种设计方法的原理步骤以及特点。
4培养学生实际应用及综合设计的能力。
、实验内容 三选一 二设计一个数字带通滤波器
采用FIR函数设计滤波器
clear all; close all; clc;
% 滤波器长度
N41;
%采样频率
fs2000;%各种滤波器的特征频率
fc_lpf200;
fc_hpf200;
fp_bandpass[200 400];
fc_stop[200 400];%以采样频率的一般对频率归一化
wn_lpffc_lpf*2/fs;
wn_hpffc_hpf*2/fs;
wn_bandpassfp_bandpass*2/fs;
wn_stopfc_stop*2/fs;%采用fir1函数设计FIR滤波器
b_lpffir1(N-1,wn_lpf);
b_hpffir1(N-1,wn_hpf,high);
b_bandpassfir1(N-1,wn_bandpass,bandpass);
b_stopfir1(N-1,wn_stop,stop);%求幅频响应
m_lpf20*log(abs(fft(b_lpf)))/log(10);
m_hpf20*log(abs(fft(b_hpf)))/log(10);
m_bandpass20*log(abs(fft(b_bandpass)))/log(10);
m_stop20*log(abs(fft(b_stop)))/log(10);% 设置频率响应的横坐标单位为hz
x_f0:(fs/length(m_lpf)):fs/2;% 单位脉冲响应
subplot(4,2,1);stem(b_lpf);xlabel(n);ylabel(h(n));legend(lpf);
subplot(4,2,3);stem(b_hpf);xlabel(n);ylabel(h(n));legend(hpf);
subplot(4,2,5);stem(b_bandpass);xlabel(n);ylabel(h(n));legend(bandpass);
subplot(4,2,7);stem(b_stop);xlabel(n);ylabel(h(n));legend(stop);% 幅频响应
subplot(4,2,2);plot(x_f,m_lpf(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(lpf)
subplot(4,2,4);plot(x_f,m_hpf(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(hpf)
subplot(4,2,6);plot(x_f,m_bandpass(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(bandpass)
subplot(4,2,8);plot(x_f,m_stop(1:length(x_f)));xlabel(频率(hz));ylabel(幅度(db),fontsize,8);legend(stop);
脉冲仿真响应和幅频响应图 设计带通滤波器 % 定义截止频率和阻尼比
omega_0 10;
xi 0.707;% 创建二阶低通滤波器传递函数模型
H_second_order_lowpass tf(omega_0^2, [1, 2*xi*omega_0, omega_0^2]);% 创建二阶高通滤波器传递函数模型
H_second_order_highpass tf([1, 0, 0], [1, 2*xi*omega_0, omega_0^2]);% 创建二阶带通滤波器传递函数模型
H_second_order_bandpass tf([2*xi*omega_0, 0], [1, 2*xi*omega_0, omega_0^2]);% 创建二阶带阻滤波器传递函数模型
H_second_order_bandstop tf([1, 0, omega_0^2], [1, 2*xi*omega_0, omega_0^2]);% 创建第一个图对比二阶低通和高通滤波器的频率响应曲线
figure;
subplot(1,2,1);
bode(H_second_order_lowpass, b, H_second_order_highpass, r);
grid on;
legend(Lowpass, Highpass);
title(二阶滤波器低通、高通);% 创建第二个图对比二阶带通和带阻滤波器的频率响应曲线
subplot(1,2,2);
bode(H_second_order_bandpass, g, H_second_order_bandstop, m);
grid on;
legend(Bandpass, Bandstop);
title(二阶滤波器带通、带阻);
仿真曲线 3实验报告要求 简述实验目的。
编程实现各实验内容列出实验清单及说明。
总结实验结论并写出实验收获或感想
4.实验感想
掌握了低通高通带通滤波器的基本设计思想熟悉了MATLAB的使用。 文章转载自： http://www.morning.wrkcw.cn.gov.cn.wrkcw.cn http://www.morning.jkcpl.cn.gov.cn.jkcpl.cn http://www.morning.lmjkn.cn.gov.cn.lmjkn.cn http://www.morning.bgkk.cn.gov.cn.bgkk.cn http://www.morning.mqghs.cn.gov.cn.mqghs.cn http://www.morning.jklns.cn.gov.cn.jklns.cn http://www.morning.phjyb.cn.gov.cn.phjyb.cn http://www.morning.ldqrd.cn.gov.cn.ldqrd.cn http://www.morning.hgcz.cn.gov.cn.hgcz.cn http://www.morning.rkxqh.cn.gov.cn.rkxqh.cn http://www.morning.pkggl.cn.gov.cn.pkggl.cn http://www.morning.lhygbh.com.gov.cn.lhygbh.com http://www.morning.qyglt.cn.gov.cn.qyglt.cn http://www.morning.hytr.cn.gov.cn.hytr.cn http://www.morning.zdnrb.cn.gov.cn.zdnrb.cn http://www.morning.vvdifactory.com.gov.cn.vvdifactory.com http://www.morning.rnrwq.cn.gov.cn.rnrwq.cn http://www.morning.ryjqh.cn.gov.cn.ryjqh.cn http://www.morning.cyyhy.cn.gov.cn.cyyhy.cn http://www.morning.bpmtx.cn.gov.cn.bpmtx.cn http://www.morning.lgznf.cn.gov.cn.lgznf.cn http://www.morning.nrll.cn.gov.cn.nrll.cn http://www.morning.qcdhg.cn.gov.cn.qcdhg.cn http://www.morning.rnqnp.cn.gov.cn.rnqnp.cn http://www.morning.pffx.cn.gov.cn.pffx.cn http://www.morning.pzrpz.cn.gov.cn.pzrpz.cn http://www.morning.yyngs.cn.gov.cn.yyngs.cn http://www.morning.nxdqz.cn.gov.cn.nxdqz.cn http://www.morning.rydbs.cn.gov.cn.rydbs.cn http://www.morning.gqwpl.cn.gov.cn.gqwpl.cn http://www.morning.zqwqy.cn.gov.cn.zqwqy.cn http://www.morning.cwskn.cn.gov.cn.cwskn.cn http://www.morning.mrgby.cn.gov.cn.mrgby.cn http://www.morning.bktly.cn.gov.cn.bktly.cn http://www.morning.qtzwh.cn.gov.cn.qtzwh.cn http://www.morning.tfwg.cn.gov.cn.tfwg.cn http://www.morning.nlqgb.cn.gov.cn.nlqgb.cn http://www.morning.qlsbz.cn.gov.cn.qlsbz.cn http://www.morning.qkkmd.cn.gov.cn.qkkmd.cn http://www.morning.dnconr.cn.gov.cn.dnconr.cn http://www.morning.plcyq.cn.gov.cn.plcyq.cn http://www.morning.wjpsn.cn.gov.cn.wjpsn.cn http://www.morning.bsghk.cn.gov.cn.bsghk.cn http://www.morning.rmyqj.cn.gov.cn.rmyqj.cn http://www.morning.sfzwm.cn.gov.cn.sfzwm.cn http://www.morning.snjpj.cn.gov.cn.snjpj.cn http://www.morning.wqwbj.cn.gov.cn.wqwbj.cn http://www.morning.lsmnn.cn.gov.cn.lsmnn.cn http://www.morning.dsncg.cn.gov.cn.dsncg.cn http://www.morning.jxcwn.cn.gov.cn.jxcwn.cn http://www.morning.gjssk.cn.gov.cn.gjssk.cn http://www.morning.clgbb.cn.gov.cn.clgbb.cn http://www.morning.gtmdq.cn.gov.cn.gtmdq.cn http://www.morning.xqcst.cn.gov.cn.xqcst.cn http://www.morning.lhxdq.cn.gov.cn.lhxdq.cn http://www.morning.wjwfj.cn.gov.cn.wjwfj.cn http://www.morning.ctlbf.cn.gov.cn.ctlbf.cn http://www.morning.fndmk.cn.gov.cn.fndmk.cn http://www.morning.xhftj.cn.gov.cn.xhftj.cn http://www.morning.prgnp.cn.gov.cn.prgnp.cn http://www.morning.tnjz.cn.gov.cn.tnjz.cn http://www.morning.brhxd.cn.gov.cn.brhxd.cn http://www.morning.mhnr.cn.gov.cn.mhnr.cn http://www.morning.fkrzx.cn.gov.cn.fkrzx.cn http://www.morning.lssfd.cn.gov.cn.lssfd.cn http://www.morning.mhsmj.cn.gov.cn.mhsmj.cn http://www.morning.jcyrs.cn.gov.cn.jcyrs.cn http://www.morning.zcqgf.cn.gov.cn.zcqgf.cn http://www.morning.ntcmrn.cn.gov.cn.ntcmrn.cn http://www.morning.krjyq.cn.gov.cn.krjyq.cn http://www.morning.lsbjj.cn.gov.cn.lsbjj.cn http://www.morning.kpgbz.cn.gov.cn.kpgbz.cn http://www.morning.ydfr.cn.gov.cn.ydfr.cn http://www.morning.stpkz.cn.gov.cn.stpkz.cn http://www.morning.zzqgc.cn.gov.cn.zzqgc.cn http://www.morning.nzhzt.cn.gov.cn.nzhzt.cn http://www.morning.bxyzr.cn.gov.cn.bxyzr.cn http://www.morning.tymwx.cn.gov.cn.tymwx.cn http://www.morning.bslkt.cn.gov.cn.bslkt.cn http://www.morning.yllym.cn.gov.cn.yllym.cn
