• 제목/요약/키워드: time variability

검색결과 1,143건 처리시간 0.024초

한국 동해에서의 해면기압, 해수면온도와 해상풍의 경험적 직교함수 분석 (Empirical Orthogonal Function Analysis of Surface Pressure, Sea Surface Temperature and Winds over the East Sea of the Korea (Japan Sea))

  • 나정열;한상규;서장원;노의근;강인식
    • 한국수산과학회지
    • /
    • 제30권2호
    • /
    • pp.188-202
    • /
    • 1997
  • 본 연구에서는 한국 동해에서의 바람장의 공간적인 분포특성을 정량적으로 파악하기 위하여 해면기압과 해상풍 자료에 EOF 분석방법을 적용하여 모드별로 고유함수의 공간적인 분포특성과 시간계수를 계산하였고 시간계수의 스펙트럼을 구하여 바람장의 공간 및 시간변동성 및 주기적인 특성을 파악하였다. 평균기압의 분포는 서고동저형의 등압선 배치를 나타내며, 기압의 최대 표준 편차는 블라디보스톡 부근에 위치하고 일본연안으로 갈수록 최저 표준편차를 갖는다. 블라디보스톡 부근은 $6.5^{\circ}C$의 해수면 온도의 최대 표준편차를 갖는 지역이기도 하다. 동서성분 해상풍의 제1모드는 분산의 $47.3\%$ 제2모드는 $30\%$를 반영하며 제1모드는 블라디보스톡 해역에서 최대 변화를 보이며 1년 및 6개월의 주기를 갖는다. 특히 제2모드에서 6개월의 주기가 탁월하다. 해상풍의 남북성분의 연평균 분포는 동해전역에서 북풍계열의 바람이 탁월하며, 특히 블라디보스톡 부근에서 최대 2.2m/s 이상의 북풍이 존재하고 3.8m/s의 표준편차 최대치가 존재하며 시간변동이 크다. 해수면 수온의 공간분포 역시 블라디보스톡 해역에서 최대치가 나타난다. 해상풍과 해면 기압의 결합 직교함수의 모드별 특성은 개별 직교함수에 비해 동서성분 해상풍의 제1모드의 기여도가 증가 하였으며, 제2모드의 기여도는 감소하는 양상을 띤다. 이로 인해 동서성분 해상풍과 해면기압의 결합 직교함수 (WU-PR)의 주기성은 해면 기압의 영향으로 제1모드에서는 1년 및 6개월의 주기가 탁월하며, 제2모드에서는 3개월의 계절변화로 해상풍의 개별 직교함수의 주기와 다르게 나타난다. 해수면 온도의 직교함수까지 포함한 3개의 결합 직교함수의 EOF분석에서는 제2모드에서 해상풍의 공간분포는 해수면 온도에 좌우되는 양상을 보인다. 스펙트럼 분석결과는 제1모드에서 1년 주기가 탁월하며 제2모드는 3개월 및 4개월의 계절변화가 나타난다.

  • PDF

지표피복변화 연구에서 최적의 공간스케일의 문제 : 가나 북부지역의 사바나 지역을 사례로 (Optimal Spatial Scale for Land Use Change Modelling : A Case Study in a Savanna Landscape in Northern Ghana)

  • 닉반드기슨;폴플렉;박수진
    • 대한지리학회지
    • /
    • 제40권2호
    • /
    • pp.221-241
    • /
    • 2005
  • 토지이용 및 지표피복변화 (Land Use and land Cover Changes, LUCC)는 지구환경변화의 원인으로 중요한 연구대상이 되고 있다. LUCC는 복잡한 사회적, 경제적, 정치적 상호작용속에서 다양한 시$\cdot$공간적 스케일에서 발생하게 된다. 따라서 LUCC를 모델화하기 위해서는 LUCC를 야기시키는 원인(driving forces)과 제한요인(constraints)들의 시$\cdot$공간적인 다양성을 이해하는 작업이 선행되어야 한다. 특히, 특정 지역에서 나타나는 LUCC의 동인을 파악하기 위해서는 스케일에 따른 그 특성의 변화를 이해하는 것이 급선무이다. 이 연구는 가나(Ghana) 북부지역의 사바나 지역을 대상으로 지난 15년간 나타난 지표피복변화의 공간적인 다양성을 파악한 뒤, 공간적 스케일을 달리하면서 나타나는 LUCC의 원인을 분석하였다. 이 과정을 통해 사바나 지역에서 LUCC 과정을 모형화하기 위한 최적의 공간적인 스케일을 규명하고자 하였다. 연구지역은 지난 15년간 인구증가의 결과로 농업생산활동이 급격하게 증가한 지역이다. 연구지역에서 나타나는 지표피복변화의 정도는 LANDSAT 위성영상에서 추출한 NDVI들을 다변량 통계분석기법을 이용하여 정량화하였다. 그리고 지표피복변화의 원인을 스케일별로 파악하기 위한 도구로 다축척 계층분석기법(multi-scale hierarchical adaptive model)을 개발$\cdot$제안하였다. 개발된 기법은 지표피복의 변화정도와 원인이 될 수 있는 공간변수들간의 상관성을 공간적인 스케일을 달리하면서 순차적으로 계산해낼 수 있는 기법이다. 이 연구에서 지표피복변화의 원인으로는 '도로에서부터의 거리', 하천으로부터의 거리', '지형특성' 의 세가지 변수를 사용하였다. 지표피복 변화정도와 위의 세가지 변수들간의 상관관계는 공간적인 범위가 10$\times$10km 이하인 경우에 높게 나타났다. 하지만 공간범위가 그 이상이 될 경우에는 그 내부에서 나타나는 다양성으로 인해 통계적인 상관성이 현격하게 낮아지는 것을 관찰할 수 있었다. 이러한 결과는 지역 및 국가 단위의 환경변화모델에서 모델의 공간적인 구성범위가 일정한 수준을 넘으면, 그 내부에서 발생하고 있는 다양성이 급격하게 증가하여 지표피복변화의 원인과 결과를 정확하게 파악하기 힘들게 된다는 것을 의미한다. 10$\times$10km의 공간적인 범위는 농업생산이 위주가 되는 사바나 지역에서는 주로 개별 마을이 차지하고 있는 공간적인 범위와 대체적으로 일치한다. 따라서 사바나 지역에서 나타나는 지표피복변화의 다양성을 고려하면서 보다 정확하게 모형화하기 위해서는 마을단위에서 나타나는 지표피복변화과정이 최소의 모델단위가 되어야 함을 시사한다.

주의력결핍 ${\cdot}$ 과잉운동장애의 Serotonin계의 개체발생적인 과정과 정신병리와의 상호관계에 관한 연구 (STUDY ON THE RELATIONSHIP BETWEEN THE ONTOGENETIC PROCESSES AND PSYCHOPATHOLOGY IN ATTENTION-DEFICIT HYPERACTIVITY DISORDER)

  • 정영;신성웅;황준원;신민섭;조수철
    • Journal of the Korean Academy of Child and Adolescent Psychiatry
    • /
    • 제12권2호
    • /
    • pp.165-178
    • /
    • 2001
  • 주의력결핍 ${\cdot}$ 과잉운동장애(이하 ADHD라 기술함)의 생화학적인 원인을 규명하고 개체발생적인 과정에 대한 이해를 하기 위하여 46명의 ADHD아동들과 18명의 발달성언어장애 아동들을 대상으로 혈장내 의 5-HT와 5-HIAA 함량을 측정하고 개체발생적인 과정에 대한 연구를 시행한 결과 다음과 같은 결과를 얻었다. 1) 대상군과 대조군에서 혈장 5-HT(serotonin) 함량비교에서는 혼합형(combined type)은 평균치 17.81 ng/ml, 주의력결핍형은 22.00ng/ml, 과잉운동 ${\cdot}$ 충동형은 24.75ng/ml 대조군의 평균치는 31.83 ng/ml로서 분산분석 결과 의미있는 차이가 관찰되었으며(F=4.33, df 3, 60, p<0.05), Scheffe씨법을 이용한 사후검증에서는 혼합형과 대조군간에는 의미있는 차이가 관찰되었으나 나머지 군간에는 의미있는 차이가 관찰되지 않았다. 2) 대상군과 대조군에서 혈장 5-HIAA함량 비교에서는 의미있는 차이가 관찰되지 않았다(F=2.08, df 3, 60, p>0.05). 3) ADHD전체군과 대조군간의 5-HT 및 5-HIAA함량비교에서는 5-HT함량의 평균치는 21.74ng/ml 이었으며, 대조군의 평균치는 31.83ng/ml로서 의미있는 차이가 관찰되었다(T=3.10, df 62, p<0.05). 그러나, 5-HIAA함량치는 의미있는 차이가 없었다(T=1.90, df 62, p>0.05). 4) ADHD군에서 TOVA소견과 5-HT 및 5-HIAA 함량간의 상관성에 있어서는 의미있는 상관성이 관찰되지 않았다. 5) ADHD군에서 DSM-IV소견과 5-HT 및 5-HIAA함량간의 상관성에 있어서도 의미있는 상관성은 관찰되지 않았다. 6) 대상군과 대조군에서 혈장 5-HT(serotonin)함량과 연령간의 상관성에 있어서는 모두 의미있는 상관성은 관찰되지 않았다(ADHD군:Pearson 상관계수 -0.20, 대조군:-0.05, 모두 p>0.05). 7) 대상군과 대조군에서 혈장 5-HIAA함량과 연령간의 상관성에 있어서도 모두 의미있는 상관성은 관찰되지 않았다(ADHD군:Pearson 상관계수 -0.17, 대조군:0.04, 모두 p>0.05). 이러한 소견으로 미루어 5-HT의 농도저하가 ADHD의 원인중의 일부로 작용할 가능성이 있으며, 5-HT 의 개체발생적인 측면에 있어서는 이상소견이 없을 가능성이 있다. 향후 생화학적인 소견에 의한 아형으로의 분류가 가능할 것으로 생각되면 이러한 분류에 의하여 약물에 대한 반응, 가족력의 차이 또는 예후에 대한 연구가 시행되어야 한다.

  • PDF

수면장애와 동반된 임상증상에 따른 심박변이도와 수면 관련 설문간의 상관성 연구 (The Relation of HRV, PSQI and IQ with Symptoms of Insomnia Patients)

  • 양동호;오경민;김보경
    • 동의신경정신과학회지
    • /
    • 제19권3호
    • /
    • pp.143-169
    • /
    • 2008
  • Objective: This study was to investigate the relation of what insomnia patients' symptoms with Heart Rate Variability(HRV), Pittsburgh Sleep Quality Index (PSQI) and Insomnia Questionaire(IQ). Methods: For this study, we carried out HRV, PSQI and IQ of 43 patients with insomnia who have come to Donguei oriental hospital of Donguei university from September 2007 to May 2008. We studied the relationship HRV with PSQI and IQ. For every symptom, this study classified the patients as none-symptom, symptom groups and analyzed each group's HRV, PSQI and IQ. Result : 1. The rate of xerostomia in female group is significantly higher than male group. 2. P1 (hypnagogic time), P3(awakening time), LF has negative correlation with age. The score of PI (hypnagogic time), P3(awakening time), P4(time of recoverable sleep) in the group whose age is more than 50 are higher than the group whose age is lesser than 50. 3. The frequency of using hypnotics(P6) in the group whose duration of onset are over 1 year was significantly higher than the score in the group whose duration of onset were below 1 year. 4. The accompanied symptoms in the group who started sleep disorder after stressor are higher than the number in the group whose duration of onset were poor-defined. 5. The drinking group have lower rate of anorexia, higher score of P1 (hypnagogic time) and TP than non-drinking group. 6. The score of P-total(general quality of sleep) have negative correlation with NNSO, HF, TP and positive correlation with LF/HF. 7. The score of P1(hypnagogic time) have negative correlation with LF. 8. The score of P5a(frequency of delayed hypnagogue) have positive correlation with LF/HF. 9. The score of P5e(frequency of cough or snoring) have negative correlation with RMSSD. 10. The score of I1a(degree of onset insomnia), I1b(degree of maintenance insomnia) have negative correlation with NN50, pNN50, HF. 11. The score of I1c(degree of termination insomnia) have negative correlation with RMSSD, NN50, pNN50, HF, TP and positive correlation with MHR, LF/HF. 12.. The score of P5i(frequency of night pain), LF/HF in the group with headache are higher than the group without headache. 13. The score of P5a(frequency of delayed hypnagogue), P5c(frequency of awakening for urine of stool) in the group with xerostomia are higher than the group without xerostomia. 14. The score of P9(degree of unrecoverable sleep) and I5a(degree of onset insomnia) in the group with chest discomfort are lower than the group without chest discomfort. 15. The score of P7(frequency of drowsiness for daytime) in the group with palpitation are lower than the group without palpitation. 16. The score of P5c(frequency of awakening for urine of stool), P5f(frequency of burning sensation), P5h(frequency of nightmare), I1c(degree of termination insomnia) in the group with burning sensation of upper trunk are higher than the group without burning sensation. 17. The score of NN50, pNN50, LF, TP in the group with anorexia are lower than the group without anorexia. 18. The score of P-total(general quality of sleep), P2(duration for sleep start) in the group with constipation are higher than the group without constipation. 19. The score of P4(recoverable sleep time) in the group with depression are higher than the group without depression. 20. In the group with anxiety, the score of P1(time of hypnagogue), P5a(frequency of onset insomnia) are lower and the score of I1b(degree of maintenance insomnia), I4(objective degree of impairment by insomnia) are higher than the group without anxiety. 21. The score of NN50, pNN50, HF, TP in the group with inertia are lower than the group without inertia. 22. The accompanied symptoms haver posirive correlation with P5f(frequency of burning sensation), P5h(frequency of nightmare), P5i(frequency of night pain) and LF/HF. Conclusion: This study provides insights into the complicated relations of the pattern of sleep disorder with many symptoms such as headache, chest discomfort, depression and any others. And this study showed that autonomic nervous system have important function in the regulation of sleep.

  • PDF

한국 서해 폐기물 투기해역의 하계 수질인자 분포특성과 장기 투기행위로 인한 영양염 농도 변화 (Variation of Nutrients due to Long-Term Effects of Ocean Dumping and Spatial Variability of Water Quality Parameters in Summer at the Ocean Waste Disposal Site Off the West Coast of Korea)

  • 고혁준;최영찬;박성은;차형기;장대수;윤한삼;이충일
    • 한국환경과학회지
    • /
    • 제22권11호
    • /
    • pp.1389-1402
    • /
    • 2013
  • This paper focuses on the impacts of waste dumping on inorganic nutrients in the dumping area of the Yellow Sea, and the effect of an governmental regulation of pollution in dumping areas. The environmental variables and parameters of the dumping and reference areas in the Yellow Sea were measured during July 2009 and analyzed. In addition, the analyzed data for inorganic nutrients over the last 10 years were obtained from the Korea Coast Guard (KCG) and the National Fisheries Research and Development Institute (NFRDI). The chemical environment of the study area revealed increases in concentrations of inorganic nutrients, Chemical Oxygen Demand (COD), and Volatile Suspended Solids (VSS) in the bottom layer. On the contrary, the pH level was decreased. Most notably, the time series data of inorganic nutrients showed gradual increase over time in the dumping area, and thus, the oligotrophic waters trend toward eutrophic waters. The increases appears to be due to the disposal of large amounts of organic waste. In recent times, the wastes disposed at the area were largely comprised of livestock wastewater, and food processing waste water. The liquefied waste, which contains an abundance of nutrients, causes a sharp increase in concentrations of inorganic nitrogen in the dumping area. On the one hand, the dumping sites have been deteriorated to such an extent that pollution has become a social problem. Consequentially, the government had a regulatory policy for improvement of marine environmental since 2007 in the dumping area. Hence, the quality of marine water in the dumping site has improved.

인공위성 고도계 자료로 추정한 동해 표층해류와 공간분포 변동성 (Satellite-altimeter-derived East Sea Surface Currents: Estimation, Description and Variability Pattern)

  • 최병주;변도성;이강호
    • 한국해양학회지:바다
    • /
    • 제17권4호
    • /
    • pp.225-242
    • /
    • 2012
  • 이 연구는 인공위성 고도계로 관측한 해수면 높이 자료를 이용하여 동해 표층해류를 생산하고, 동해 전체 영역에 대하여 동시성 있는 표층해류 분포를 동해 해류 정보 사용자에게 제공하기 위한 최초의 시도이다. 동해 전 영역에서 인공위성 고도계와 연안 조위관측소의 해수면 높이 자료를 동시에 얻을 수 있으므로 준실시간으로 넓은 해역에 대하여 동시성 있는 표층 지형류의 산출이 가능하다. 산출된 동해 표층 지형류로부터 주요 해류의 위치와 세기 그리고 중규모 이상의 소용돌이 발달 양상을 살펴볼 수 있다. 따라서 이들 해류의 이름과 평균적인 위치를 알 수 있도록 동해 해양지명에 대한 명칭과 위치를 기술하고, 개념적인 해류도를 제시하였다. 동해 해류정보가 실제 실용화될 수 있도록 인공위성 고도계 자료를 이용하여 산출한 지형류의 월 계절 연도별 해류 분포 예를 들고, 각 해류 분포를 설명하였다. 또한 시 공간적으로 변화하는 동해 표층해류 분포 형태를 객관적으로 분류하기 위하여 16년(1993~2008년)간의 표층해류 자료를 경험직교함수(Empirical Orthogonal Function, EOF)를 이용하여 분석하였다. EOF분석 제1모드(mode)는 주로 한국 동해안을 따라 북쪽으로 흐르는 동한난류와 야마도분지 남서부 시계방향 순환의 강화 또는 약화를 나타냈다. 제2모드는 동한난류가 동해 남부를 가로지르며 사행하는 정도를 나타냈으며, 해류 사행의 파장은 약 300 km이었다. 제1모드와 제2모드가 모두 해류의 연간 변동성을 나타냈으며, 제1모드와 제2모드의 시간계수에 따라 동해 표층해류 분포를 관성 경계류 패턴(pattern), 대마난류 패턴, 사행 패턴, 외해분지류 패턴으로 분류할 수 있다.

걸음걸이 분석 기반의 파킨슨병 분류를 위한 특징 추출 (Features Extraction for Classifying Parkinson's Disease Based on Gait Analysis)

  • 이상홍;임준식;신동근
    • 인터넷정보학회논문지
    • /
    • 제11권6호
    • /
    • pp.13-20
    • /
    • 2010
  • 본 논문은 걸음걸이 분석 기반의 특징 추출과 NEWFM(Neural Network with Weighted Fuzzy Membership Functions)을 이용하여 건강한 사람의 족압(foot pressure)과 파킨슨병 환자의 족압으로부터 건강한 사람과 파킨슨병 환자를 분류하는 방안을 제시하고 있다. NEWFM에서 사용할 입력을 추출하기 위해서 첫 번째 단계에서는 PhysioBank에서 제공하는 족압 데이터와 시간에 따른 족압의 변화를 이용하여 각각 4개의 특징을 추출하였다. 두 번째 단계에서는 웨이블릿 변환(wavelet transform, WT)을 이용하여 이전 단계에서 추출한 8개의 특징으로부터 웨이블릿 계수를 추출하였다. 마지막 단계에서는 추출된 웨이블릿 계수들을 이용하여 통계적 기법인 주파수 분포와 주파수 변동량을 이용하여 40개의 입력을 추출하였다. NEWFM은 족압 데이터로부터 8개의 특징을 추출하여 건강한 사람과 파킨슨병 환자를 분류하였을 때 왼쪽 족압과 오른쪽 족압의 차를 이용한 특징과 시간에 따른 족압의 변화에 대한 차를 이용한 특징의 경우에 높은 정확도(accuracy)가 나타났다. 이러한 결과를 통하여 걸음걸이에 있어서 질질 끄는 특징을 보이는 파킨슨병 환자의 양쪽 족압의 차가 건강한 사람의 양쪽 족압의 차보다는 상대적으로 적다는 특징을 본 실험을 통해 확인할 수 있었다.

히마와리 위성의 산불방사열에너지 자료를 이용한 산불배출가스 추정: 2017년 삼척 및 강릉 산불을 사례로 (Estimation of Fire Emissions Using Fire Radiative Power (FRP) Retrieved from Himawari-8 Satellite)

  • 김대선;원명수;이양원
    • 대한원격탐사학회지
    • /
    • 제33권6_1호
    • /
    • pp.1029-1040
    • /
    • 2017
  • 산불은 다량의 온실가스를 대기 중으로 방출하는 자연재해로서, 이를 효율적으로 감시하기 위해서는 정지궤도 위성의 산불방사열에너지(fire radiative power, FRP)를 활용하는 방법이 필요하다. 본 연구에서는 2017년 5월 6일에 발생한 우리나라 삼척과 강릉 산불을 사례로, 히마와리 위성의 중적외 채널자료를 이용하여 FRP를 산출하였으며, 이를 통해 MODIS(Moderate Resolution Imaging Spectroradiometer)의 제한적인 시간해상도로는 관측이 불가능한 10분 간격의 산불 피해강도의 실시간 모니터링이 가능함을 확인하였다. 또한 히마와리 FRP를 이용하여 강릉 산불의 배출가스를 계산하였으며, 에어코리아 실측치와 비교하였을 때 거리 차에 의한 1~3시간의 지연현상과 함께, 산불배출가스의 시계열 패턴이 매우 잘 일치함을 알 수 있었다. 또한 선행연구에서 고해상도 영상분석을 통해 제시한 산불배출가스 추정량과 비교하였을 때, 100 ha당 배출량이 삼척은 약 12%, 강릉은 약 2%의 차이로 매우 유사한 결과를 나타냈다. 이는 산불 피해면적과 피해강도에 대한 직접적인 분석 없이도, 정지궤도 위성의 FRP만을 이용하여 산불배출가스의 정밀한 추정이 가능함을 의미한다. 이 연구는 향후 발사될 우리나라 정지궤도 기상위성인 GK-2A(Geostationary Korea Multi-Purpose Satellite-2A)의 산불배출가스 추정 및 에어로솔 산출에 활용될 수 있을 것으로 사료된다.

USE OF NEAR-INFRARED SPECTROSCOPY TO PREDICT OIL CONTENT COMPONENTS AND FATTY ACID COMPOSITION IN OLIVE FRUIT

  • Lorenzo, Leon-Moreno;Ana, Garrido-Varo;Luis, Rallo-Romero
    • 한국근적외분광분석학회:학술대회논문집
    • /
    • 한국근적외분광분석학회 2001년도 NIR-2001
    • /
    • pp.1512-1512
    • /
    • 2001
  • The University of Cordoba conducts since 1991 a breeding program to obtain new olive cultivars from intraspecific crosses. The objective is to obtain new early bearing and high-quality cultivars. In plant breeding, many seedlings must be tested to increased the chance of getting desirable genotypes. Therefore, fast, cheap and accurate methods of analysis are necessary. The conventional laboratory techniques are costly and time-consuming. Near Infrared Spectroscopy (NIRS) can satisfy the characteristics requested by plant breeders and offers many advantages such as the simultaneous analysis of many traits and cheap cost. The objective of this work was to asses the performance of NIRS to estimate oil fruit components (fruit weight, flesh moisture, flesh/stone ratio and oil flesh content in dry weight basis) and fatty acid composition in olive fruit. Genotypes from reciprocal crosses between ‘Arbequina’, ‘Frantoio’ and ‘Picual’ cultivars have been used in this study. A total of 287 samples, each from a single plant, were scanned using a DA-7000 Diode Array VIS/NIR Analysis System (Perten Instruments), which covers the visible and NIR range from 400-1700 nm. All samples were analysed for fatty acid composition (gas chromatography) and 220 for oil fruit components (oil content by nuclear magnetic resonance), 70% and 30% of samples were randomly assign for the calibration and validation sets respectively. The preliminary results shows that calibration for palmitic, oleic and linoleic acids were highly accurate with calibration and validation values of $r^2$ from 0.85 to 0.95 and 0.76 to 0.91 respectively. Calibration for palmitoleic and estearic acids were less accurate, probably because of the narrow range of variability available for these fatty acids. For the oil fruit components, calibration were high accurate for flesh moisture and oil flesh content in dry weight basis ($r^2$ higher than 0.90 in both calibration and validation sets) and less accurate for the other characteristics evaluated. The first results obtained indicate that NIRS analysis could be an ideal technique to reduce the cost, time and chemical wasted necessary to evaluate a large number of genotypes and it is accurate enough to use for pre-selecting genotypes in a breeding program.

  • PDF

WZ Cephei: A Dynamically Active W UMa-Type Binary Star

  • Jeong, Jang-Hae;Kim, Chun-Hwey
    • Journal of Astronomy and Space Sciences
    • /
    • 제28권3호
    • /
    • pp.163-172
    • /
    • 2011
  • An intensive analysis of 185 timings of WZ Cep, including our new three timings, was made to understand the dynamical picture of this active W UMa-type binary. It was found that the orbital period of the system has complexly varied in two cyclical components superposed on a secularly downward parabola over about 80y. The downward parabola, corresponding to a secular period decrease of $-9.{^d}97{\times}10^{-8}y^{-1}$, is most probably produced by the action of both angular momentum loss (AML) due to magnetic braking and mass-transfer from the massive primary component to the secondary. The period decrease rate of $-6.^{d}72{\times}10^{-8}y^{-1}$ due to AML contributes about 67% to the observed period decrease. The mass flow of about $5.16{\times}10^{-8}M_{\odot}y^{-1}$ from the primary to the secondary results the remaining 33% period decrease. Two cyclical components have an $11.^{y}8$ period with amplitude of $0.^{d}0054$ and a $41.^{y}3$ period with amplitude of $0.^{d}0178$. It is very interesting that there seems to be exactly in a commensurable 7:2 relation between their mean motions. As the possible causes, two rival interpretations (i.e., light-time effects (LTE) by additional bodies and the Applegate model) were considered. In the LTE interpretation, the minimum masses of $0.30M_{\odot}$ for the shorter period and $0.49M_{\odot}$ for the longer one were calculated. Their contributions to the total light were at most within 2%, if they were assumed to be main-sequence stars. If the LTE explanation is true for the WZ Cep system, the 7:2 relation found between their mean motions would be interpreted as a stable 7:2 orbit resonance produced by a long-term gravitational interaction between two tertiary bodies. In the Applegate model interpretation, the deduced model parameters indicate that the mechanism could work only in the primary star for both of the two period modulations, but could not in the secondary. However, we couldn't find any meaningful relation between the light variation and the period variability from the historical light curve data. At present, we prefer the interpretation of the mechanical perturbation from the third and fourth stars as the possible cause of two cycling period changes.