clc;clear;
name1=['chan','.txt'];
[site,mjd,ah,aw,zhd,zwd,tm]=textread(name1,'%s%s%s%s%s%s%s%*[^\n]','delimiter',' ');%提取前5列資料
n=length(tm);
for i=1:n%在1-n行資料內判斷資料是否為空，如果為空則跳出本次迴圈
tm(i,:)=str2num(tm);
s=i/4+1;
j=floor(s);
if mod(i,4)==1
tm1(j,:)=tm(i,:);
elseif mod(i,4)==2
tm2(j,:)=tm(i,:);
elseif mod(i,4)==3
tm3(j,:)=tm(i,:);
elseif mod(i,4)==0
tm4(j-1,:)=tm(i,:);
endendgpt3_tm=xlsread('gpt3_tm.xlsx');
sou=xlsread('54161_2018.xlsx');
sou_tm=sou(:,1);
nn=length(sou_tm);
for i=1:nn%在1-n行資料內判斷資料是否為空，如果為空則跳出本次迴圈
tm(i,:)=str2num(tm);
s=i/2+1;
j=floor(s);
if mod(i,2)==1
sou_tm1(j,:)=sou_tm(i,:);
elseif mod(i,2)==0
sou_tm2(j-1,:)=sou_tm(i,:);
endendfor i=1:365
vmftm(i,:)=(tm1(i,:)+tm2(i,:)+tm3(i,:)+tm4(i,:))/4;
sou_tm(i,:)=(sou_tm1(i,:)+sou_tm2(i,:))/2;
errorvmf(i,:)=vmftm(i,:)-gpt3_tm(i,5);
errorsou(i,:)=sou_tm(i,:)-gpt3_tm(i,5);
endmean(errorvmf)
rms(errorvmf)
mean(errorsou)
rms(errorsou)
php的精度計算問題（bcadd和bcsub）
我們在進行php開發的時候經常會遇到浮點型的問題，特別是涉及金額的部分，常常需要進行加減運算。當小數點的位數比較多的時候，往往容易犯一些很低階的錯誤。這裡記錄一下php的精度計算和封裝的小demo。1 概念解釋 php高精度計算問題 2 高精度數值對比大小問題 臨時發乙個專案遇到的php浮點計算問題...
計算從整數1到n之間，所有1的個數的和
整數中1的個數 計算從整數1到n之間，所有1的個數的和 c include using namespace std int numberof1between1andn solution int n return ones int main 通過使用乙個 位置乘子m 遍歷數字的位置,m 分別為1,10...
Linux雲計算的重定向和管道（1）
fd簡介 file descriptors fd 檔案描述符，檔案控制代碼 程序使用檔案描述符來管理開啟的檔案 date 檢視當前系統的時間 date f5.txt 將系統的時間固定寫到乙個資料夾中。程序管道piping 管道命令可以將多條命令組合起來,一次性完成複雜的處理任務。語法 command...