 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2999 | 回復(fù): 23 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
maomao1210金蟲 (正式寫手)
|
[交流]
【分享】暑假閑來無事，與蟲蟲共同進(jìn)步 已有17人參與
|
||
|
暑假閑來無事, 有對寫學(xué)術(shù)型代碼感興趣的蟲子可以一起共同進(jìn)步 [ Last edited by maomao1210 on 2010-7-28 at 18:53 ] |
回復(fù)此樓» 猜你喜歡
一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑!
已經(jīng)有4人回復(fù)
-
298求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有11人回復(fù)
-
0856求調(diào)劑
已經(jīng)有7人回復(fù)
-
356求調(diào)劑
已經(jīng)有3人回復(fù)
-
332求調(diào)劑
已經(jīng)有7人回復(fù)
-
348求調(diào)劑
已經(jīng)有4人回復(fù)
-
305求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有13人回復(fù)
-
調(diào)劑310
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 閑來無事不從容
已經(jīng)有16人回復(fù)
- 【原創(chuàng)】請教一下夫妻兩人同專業(yè)共同申請海外博士后的問題
已經(jīng)有45人回復(fù)
sxlion811
金蟲 (正式寫手)
- 程序強(qiáng)帖: 3
- 應(yīng)助: 12 (小學(xué)生)
- 貴賓: 0.128
- 金幣: 1079.5
- 散金: 2188
- 紅花: 3
- 帖子: 870
- 在線: 445.3小時(shí)
- 蟲號: 479151
- 注冊: 2007-12-14
- 專業(yè): 食品科學(xué)
★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
|
既然大家開始有時(shí)間了，有個(gè)問題請教大家，實(shí)現(xiàn)語言不限。 有一組數(shù)據(jù)，每一行有四個(gè)變量id1-id4，如果任意有兩個(gè)值相等的話， 就認(rèn)為他們是同組數(shù)據(jù)。 比如說，取其中6行數(shù)據(jù)說明。 1 11 21 31 1 12 22 31 2 11 22 32 2 13 21 32 3 14 23 32 4 11 22 33 確定為同一組值后，需要對他們進(jìn)行編號。 例如第1，2行，有相同的1，31，所以為相同組，編號1； 又3和4，3和6行分別有相同的2，32和11，22，所以3，4，6為同組編號2。 第5行不和其他行有兩個(gè)或以上相同的值，編號3。 大家有什么好的解決方法嗎？由于數(shù)據(jù)可能上千條，所以需要考慮效率的問題。 |
4樓2010-07-22 20:05:44
feichengxyz
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 8530.5
- 散金: 60
- 紅花: 4
- 帖子: 1162
- 在線: 908.4小時(shí)
- 蟲號: 425858
- 注冊: 2007-07-28
- 性別: GG
- 專業(yè): 催化化學(xué)
3樓2010-07-22 18:28:59
sxlion811
金蟲 (正式寫手)
- 程序強(qiáng)帖: 3
- 應(yīng)助: 12 (小學(xué)生)
- 貴賓: 0.128
- 金幣: 1079.5
- 散金: 2188
- 紅花: 3
- 帖子: 870
- 在線: 445.3小時(shí)
- 蟲號: 479151
- 注冊: 2007-12-14
- 專業(yè): 食品科學(xué)
★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
|
還是簡化一下吧，那個(gè)有點(diǎn)復(fù)雜。 有一列數(shù)據(jù)，見附，任何一行數(shù)據(jù)有相等的sort1或sort2，分成一組。 如：1，2，14，18，3，28為一組。 4，15，6，11，10，15，20，17，24，12，22，7，30為第二組。 依此分組。 sort1 sort2 1 2 1 3 1 14 2 14 2 18 3 28 4 15 4 30 6 10 6 15 6 16 6 27 7 17 8 30 10 20 11 15 12 17 12 22 12 24 14 18 16 27 17 30 19 23 24 30 好像要用到圖論中的方法。 |
5樓2010-07-23 10:10:43
高人能給指點(diǎn)下不？
★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
|
第一次編矩陣?yán)镉袕?fù)數(shù)的程序！夠煩人啊！ 思路：在外部定義一系列變量，t,w,f,pi為實(shí)型，fi為復(fù)型；矩陣h0為復(fù)型，要求它的特征值，h0由內(nèi)部子程序調(diào)用得到，它又由幾個(gè)小矩陣：c,a,b組合得到。首先將h0定義為全0矩陣，然后再對它內(nèi)部片段矩陣賦值，c,a,b按照一定順序排列在h0里面，得到h0的值后將其傳回主程序求特征值。大家能否幫忙看下，不知道h0定義為復(fù)型好不好啊 program main use IMSL implicit none real,parameter::pi=3.1415927 complex,parameter::fi=(0.0,1.0) complex q real t,f,w complex, allocatable::h0(:,: ) real, allocatable:: eigenvalue(: ) integer :: i,n read*,n allocate(h0(4*n,4*n)) allocate(eigenvalue(n)) w=0.0 t=1.0 f=0.25 call hmatr(h0,n) eigenvalue=eig(h0) do i=1,4*n write(*,"('eigenvalue=',f10.7)" )eigenvalue(i) end do stop contains subroutine hmatr(h0,n) integer n,i,j,n0 complex h0 complex,dimension(:,: )::c(4,4),a(4,4),b(4,4) h0=0.0 a=reshape((/0.0,0.0,0.0,t,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/),(/4,4/)) b=transpose(a) do n0=1,2*n-1,2 c=((/w,t*cmplx(exp(-fi*(-pi+2*pi/3*f*(1.5*n0+1/4)))),0,0,t*exp(fi*(-pi+2*pi/3*f*(1.5*n0+1/4))),w,t,0,0,t,w,t*exp(fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),0,0,t*exp(-fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),w/),(/4,4/)) j=2*n0-1 do j=1,4*n-3,4 h0(j:j+3,j:j+3)=c if(j==1) then h0(j:j+3,j+4:j+7)=b; h0(j:j+3,4*n-3:4*n)=a; else if(j>=5.and.j<4*n-3) then h0(j:j+3,j+4:j+7)=b; h0(j:j+3,j-4:j-1)=a; else h0(4*n-3:4*n,1:4)=b; h0(j:j+3,j-4:j-1)=a; end if end do end do end subroutine hmatr end 出錯(cuò)信息： 1）Error: The shape matching rules of actual arguments and dummy arguments have been violated. [H0] call hmatr(h0,n) -----------^ 2）E:\matlabxuexi\yyuu.f90(36) : Error: A constant or named constant is required in this context. c=((/w,t*cmplx(exp(-fi*(-pi+2*pi/3*f*(1.5*n0+1/4)))),0,0,t*exp(fi*(-pi+2*pi/3*f*(1.5*n0+1/4))),w,t,0,0,t,w,t*exp(fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),0,0,t*exp(-fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),w/),(/4,4/)) ---------^ 3）E:\matlabxuexi\yyuu.f90(36) : Error: An INTEGER or REAL data type is required in this context. c=((/w,t*cmplx(exp(-fi*(-pi+2*pi/3*f*(1.5*n0+1/4)))),0,0,t*exp(fi*(-pi+2*pi/3*f*(1.5*n0+1/4))),w,t,0,0,t,w,t*exp(fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),0,0,t*exp(-fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),w/),(/4,4/)) ---------^ 4）E:\matlabxuexi\yyuu.f90(39) : Error: The syntax of this substring is invalid. [H0] h0(j:j+3,j:j+3)=c -------^ Error executing df.exe. yyuu.obj - 4 error(s), 0 warning(s) [ Last edited by 田山東 on 2010-7-23 at 11:32 ] |
6樓2010-07-23 11:30:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +7 | 蕉蕉123 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級已過 +7 | 幸運(yùn)哩哩 2026-03-22 | 11/550 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考博] 26申博 +3 | 加油沖��！ 2026-03-26 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|