应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【分享】暑假闲来无事，与虫虫共同进步
241/3返回列表123下一页
查看: 3288  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

maomao1210

金虫 (正式写手)

[交流] 【分享】暑假闲来无事，与虫虫共同进步 已有17人参与

暑假闲来无事,
有对写学术型代码感兴趣的虫子可以一起共同进步

[ Last edited by maomao1210 on 2010-7-28 at 18:53 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-07-22 16:41:55
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
高手啊，拜一个
回复此楼
2楼2010-07-22 17:04:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

feichengxyz

铁杆木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1):感谢无私帮助，可否新开一帖，跟maomao有所区分，方便大家求助！ 2010-07-22 21:09:16
插播一广告：本人提供Matlab代码编写服务
赞一下(2人) 回复此楼
3楼2010-07-22 18:28:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sxlion811

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
既然大家开始有时间了，有个问题请教大家，实现语言不限。

有一组数据，每一行有四个变量id1-id4，如果任意有两个值相等的话，
就认为他们是同组数据。
比如说，取其中6行数据说明。
1 11 21 31
1 12 22 31
2 11 22 32
2 13 21 32
3 14 23 32
4 11 22 33

确定为同一组值后，需要对他们进行编号。
例如第1，2行，有相同的1，31，所以为相同组，编号1；
又3和4，3和6行分别有相同的2，32和11，22，所以3，4，6为同组编号2。
第5行不和其他行有两个或以上相同的值，编号3。

大家有什么好的解决方法吗？由于数据可能上千条，所以需要考虑效率的问题。
赞一下(1人) 回复此楼
开心努力一辈子
4楼2010-07-22 20:05:44
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sxlion811

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
还是简化一下吧，那个有点复杂。
有一列数据，见附，任何一行数据有相等的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

好像要用到图论中的方法。
赞一下(1人) 回复此楼
开心努力一辈子
5楼2010-07-23 10:10:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

高人能给指点下不？

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
第一次编矩阵里有复数的程序！够烦人啊！
思路：在外部定义一系列变量，t,w,f,pi为实型，fi为复型；矩阵h0为复型，要求它的特征值，h0由内部子程序调用得到，它又由几个小矩阵：c,a,b组合得到。首先将h0定义为全0矩阵，然后再对它内部片段矩阵赋值，c,a,b按照一定顺序排列在h0里面，得到h0的值后将其传回主程序求特征值。大家能否帮忙看下，不知道h0定义为复型好不好啊

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

出错信息：
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 ]
赞一下(1人) 回复此楼
everythinghasitsseason.enjoyyourlife.
6楼2010-07-23 11:30:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhoudb

至尊木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
十分钦佩。

是否愿为用等效电路处理电化学数据做一拟合程序？
赞一下(1人) 回复此楼
7楼2010-07-25 11:50:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sxlion811

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+1):谢谢说明。 2010-07-26 07:58:05
随便说一下，上面的两个问题，经人指点后得知，需要用到深度优先搜索算法，谢谢。
赞一下(2人) 回复此楼
开心努力一辈子
8楼2010-07-25 21:06:46
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
maomao你怎么不在群里冒头，在这里出风头啊。
赞一下(1人) 回复此楼
9楼2010-07-25 22:24:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jaey

怎样联系？

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我的邮箱：zdj_zhangdaijuan@sohu.com
赞一下(1人) 回复此楼
10楼2010-07-26 08:01:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maomao1210 的主题更新
241/3返回列表123下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 人气不行了 +9 fansofjerry 2026-08-21 9/450 2026-08-24 17:48 by zhaozhengxin
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭