-
求標準偏差
> function c=myfunction(x)
> [m,n]=size(x)
> t=0
> for i=1:numel(x)
> t=t+x(i)*x(i)
> end
> c=sqrt(t/(m*n-1))
function c=myfunction(x)
[m,n]=size(x)
t=0
for i=1:m
for j=1:n
t=t+x(i,j)*x(i,j)
end
end
c=sqrt(t/(m*n-1
標簽：
gt
myfunction
function
numel
上傳時間：
2014-09-03
上傳用戶：jjj0202
-
動態規劃的方程大家都知道，就是
f[i,j]=min{f[i-1,j-1],f[i-1,j],f[i,j-1],f[i,j+1]}+a[i,j]
但是很多人會懷疑這道題的后效性而放棄動規做法。
本來我還想做Dijkstra，后來變了沒二十行pascal就告訴我數組越界了……（dist:array[1..1000*1001
div 2]...）
無奈之余看了xj_kidb1的題解，剛開始還覺得有問題，后來豁然開朗……
反復動規。上山容易下山難，我們可以從上往下走，最后輸出f[n][1]。
xj_kidb1的一個技巧很重要，每次令f[i][0]=f[i][i],f[i][i+1]=f[i][1]（xj_kidb1的題解還寫錯了）
標簽：
動態規劃
方程
家
上傳時間：
2014-07-16
上傳用戶：libinxny
-
一本很特別的matlab教材，慢慢探索吧
本書第一章系統講述MATLAB6.5的集成工作平臺，引導讀者初會MATLAB基本用法。第二、三章系統敘述MATLAB6.1的四大基本數據類型（數值數組、字符串數組、元胞數組、構架數組），至于其它擴展數據類型（函數句柄、符號數據、內聯函數、unit數組、稀疏類）則另辟章節專述。此后，本書用九個獨立章分述MATLAB6.5的數值計算、符號計算、函數和數據可視、面向對象編程、GUI交互操作界面設計、EXE獨立應用程序生成、實現不同軟件平臺交互的API、M-book數據圖形文字環境集成等八大通用功能。本書用專門的一章深入淺出地闡明SIMULINK的分層建模、仿真功能、與MATLAB交互的功能。
全書包含390多個計算范例。所有算例的程序是可靠、完整的。讀者可以在自己的MATLAB環境中完全準確地重現本書所提供的算例結果。算例的舉一反三，將使讀者很快掌握要領，從模仿走向靈活應用。
書中正文（包括算例）所涉及的指令全部羅列在附錄的“A索引”中，用戶很容易根據關鍵字符找到相關的正文說明和算例應用。
標簽：
matlab
教材
上傳時間：
2013-12-27
上傳用戶：lmeeworm
-
Euler函數:
m = p1^r1 * p2^r2 * …… * pn^rn ai >= 1 , 1 <= i <= n
Euler函數：
定義：phi(m) 表示小于等于m并且與m互質的正整數的個數。
phi(m) = p1^(r1-1)*(p1-1) * p2^(r2-1)*(p2-1) * …… * pn^(rn-1)*(pn-1)
= m*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pn)
= p1^(r1-1)*p2^(r2-1)* …… * pn^(rn-1)*phi(p1*p2*……*pn)
定理：若(a , m) = 1 則有 a^phi(m) = 1 (mod m) 即a^phi(m) - 1 整出m
在實際代碼中可以用類似素數篩法求出
for (i = 1 i < MAXN i++)
phi[i] = i
for (i = 2 i < MAXN i++)
if (phi[i] == i)
{
for (j = i j < MAXN j += i)
{
phi[j] /= i
phi[j] *= i - 1
}
}
容斥原理：定義phi(p) 為比p小的與p互素的數的個數
設n的素因子有p1, p2, p3, … pk
包含p1, p2…的個數為n/p1, n/p2…
包含p1*p2, p2*p3…的個數為n/(p1*p2)…
phi(n) = n - sigm_[i = 1](n/pi) + sigm_[i!=j](n/(pi*pj)) - …… +- n/(p1*p2……pk)
= n*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pk)
標簽：
Euler
lt
phi
函數
上傳時間：
2014-01-10
上傳用戶：wkchong
-
//Euler 函數前n項和
/*
phi(n) 為n的Euler原函數
if( (n/p) ％ i == 0 ) phi(n)=phi(n/p)*i
else phi(n)=phi(n/p)*(i-1)
對于約數：divnum
如果i|pr[j] 那么 divnum[i*pr[j]]=divsum[i]/(e[i]+1)*(e[i]+2) //最小素因子次數加1
否則 divnum[i*pr[j]]=divnum[i]*divnum[pr[j]] //滿足積性函數條件
對于素因子的冪次 e[i]
如果i|pr[j] e[i*pr[j]]=e[i]+1 //最小素因子次數加1
否則 e[i*pr[j]]=1 //pr[j]為1次
對于本題：
1. 篩素數的時候首先會判斷i是否是素數。
根據定義，當 x 是素數時 phi[x] = x-1
因此這里我們可以直接寫上 phi[i] = i-1
2. 接著我們會看prime[j]是否是i的約數
如果是，那么根據上述推導，我們有：phi[ i * prime[j] ] = phi[i] * prime[j]
否則
phi[ i * prime[j] ] = phi[i] * (prime[j]-1)
(其實這里prime[j]-1就是phi[prime[j]]，利用了歐拉函數的積性)
經過以上改良，在篩完素數后，我們就計算出了phi[]的所有值。
我們求出phi[]的前綴和
*/
標簽：
phi
Euler
else
函數
上傳時間：
2016-12-31
上傳用戶：gyq
-
遙控解碼通過電腦串口顯示
/* 晶振：11.0569MHz */
#include <REGX52.h>
#define uchar unsigned char
uchar data IRcode[4] //定義一個4字節的數組用來存儲代碼
uchar CodeTemp //編碼字節緩存變量
uchar i,j,k //延時用的循環變量
sbit IRsignal=P3^2 //HS0038接收頭OUT端直接連P3.2(INT0)
/**************************延時0.9ms子程序**********************/
void Delay0_9ms(void)
{uchar j,k
for(j=18 j>0 j--)
for(k=20 k>0 k--)
}
/***************************延時1ms子程序**********************/
void Delay1ms(void)
{uchar i,j
for(i=2 i>0 i--)
for(j=230 j>0 j--)
}
標簽：
uchar
unsigned
11.0569
include
上傳時間：
2013-12-12
上傳用戶：Breathe0125
-
Instead of finding the longest common
subsequence, let us try to determine the
length of the LCS.
Then tracking back to find the LCS.
Consider a1a2…am and b1b2…bn.
Case 1: am=bn. The LCS must contain am,
we have to find the LCS of a1a2…am-1 and
b1b2…bn-1.
Case 2: am≠bn. Wehave to find the LCS of
a1a2…am-1 and b1b2…bn, and a1a2…am and
b b b
b1b2…bn-1
Let A = a1 a2 … am and B = b1 b2 … bn
Let Li j denote the length of the longest i,g g
common subsequence of a1 a2 … ai and b1 b2
… bj.
Li,j = Li-1,j-1 + 1 if ai=bj
max{ L L } a≠b i-1,j, i,j-1 if ai≠j
L0,0 = L0,j = Li,0 = 0 for 1≤i≤m, 1≤j≤n.
標簽：
the
subsequence
determine
Instead
上傳時間：
2013-12-17
上傳用戶：evil
-
//初始化
initscr()
//獲得屏幕尺寸
getmaxyx(stdscr, h, w)
//畫背景
for(i=0 i<h i++)
for(j=0 j<w j++){
mvaddch(i, j, ACS_CKBOARD)
}
refresh()
//建立窗口
pad = newpad(80, 128)
for(i=0 i<80 i++){
char line[128]
sprintf(line, "This line in pad is numbered d\n", i)
mvwprintw(pad, i, 0, line)
}
//刷新屏幕
refresh()
prefresh(pad, 0, 1, 5, 10, 20, 45)
for(i=0 i<50 i++){
prefresh(pad, i+1, 1, 5, 10, 20, 45)
usleep(30000)
}
//等待按鍵
getch()
標簽：
getmaxyx
initscr
stdscr
for
上傳時間：
2014-08-30
上傳用戶：龍飛艇
-
實驗源代碼
//Warshall.cpp #include<stdio.h> void warshall(int k,int n) { int i , j, t; int temp[20][20]; for(int a=0;a<k;a++) { printf("請輸入矩陣第%d 行元素:",a); for(int b=0;b<n;b++) { scanf ("%d",&temp[a][b]); } } for(i=0;i<k;i++){ for( j=0;j<k;j++){ if(temp[ j][i]==1) { for(t=0;t<n;t++) { temp[ j][t]=temp[i][t]||temp[ j][t]; } } } } printf("可傳遞閉包關系矩陣是:\n"); for(i=0;i<k;i++) { for( j=0;j<n;j++) { printf("%d", temp[i][ j]); } printf("\n"); } } void main() { printf("利用 Warshall 算法求二元關系的可傳遞閉包\n"); void warshall(int,int); int k , n; printf("請輸入矩陣的行數 i: "); scanf("%d",&k);
四川大學實驗報告 printf("請輸入矩陣的列數 j: "); scanf("%d",&n); warshall(k,n); }
標簽：
warshall
離散
實驗
上傳時間：
2016-06-27
上傳用戶：梁雪文以
-
Computation of loudness (Zwicker model) according to ISO 532B / DIN 45631 norms.
This model is valid for steady sounds.
Code based on BASIC program published in the following article:
Program for calculating loudness according to DIN 45 631 (ISO 532B)",
E.Zwicker and H.Fastl, J.A.S.J (E) 12, 1 (1991).
標簽：
計算
聲學
上傳時間：
2016-11-14
上傳用戶：zztony16
