24小时热门版块排行榜    

查看: 2514  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

luqingq

银虫 (初入文坛)

[求助] c++求复数组成的n阶行列式的值的代码!急!急!

c++求复数组成的n阶行列式的值的代码!急!急!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)


jjdg(金币+1): 感谢参与 2011-10-15 01:57:20
引用回帖:
7楼: Originally posted by luqingq at 2011-10-14 15:40:24:
谢谢  
for( M = I;wave1[M][J].real()==0.0 && wave1[M][J].imag()==0.0;M++)cout<<"M="<<M<<endl;
这样应该可以限制M的值只能到两百吧
为什么这样比较浮点数是不对的 ...

1. 还是不能

2. 可以去网上搜索一下“浮点数比较”

最后还是建议用别人写好的库
8楼2011-10-14 16:09:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

luqingq

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by 锐利的碎片 at 2011-10-13 19:09:41:
非要c++吗?

恩 最好是c++
c应该也可以
3楼2011-10-14 13:19:01
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

luqingq

银虫 (初入文坛)

for( I=0,J=0;I         {
                if( wave1[I][J].real()==0.0 && wave1[I][J].imag()==0.0)
                {
                        for( M = I;wave1[I][J].real()==0.0 && wave1[I][J].imag()==0.0;M++)cout<<"M="<                         if(  M == m_system )
                        {
                                value_1_re=0.0;
                                value_1_im=0.0;
                                flag1=1;
                        }
                        else
                                for( L=J;L                                 {
                                        temp1=wave1[I][L];
                                        wave1[I][L]=wave1[M][L];
                                        wave1[M][L]=temp1;
                                }
                                if(flag1!=1)
                                        k*=(-1);
                                else
                                        break;
                }
                for( s=m_system-1;s>I;s--)
                {
                        xx=wave1[s][J];
                        for( T = J; T < m_system; T++)
                                wave1[s][T]=multiplication(wave1[I][T],multiplication(xx,conjg(wave1[I][J]))/multiplication(wave1[I][J],conjg(wave1[I][J])));
                }
        }
        if(flag1!=1)
        {
                for( I=0;I                         ff1=multiplication(ff1,wave1[I][I]);
                value_1_re=k * ff1.real();
                value_1_im=k * ff1.imag();
        }//value_1的值的计算
4楼2011-10-14 13:20:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

luqingq

银虫 (初入文坛)

引用回帖:
4楼: Originally posted by luqingq at 2011-10-14 13:20:26:
for( I=0,J=0;I<m_system && J<m_system;I++,J++)
        {
                if( wave1[I][J].real()==0.0 && wave1[I][J].imag()==0.0)
                {
                        for( M = I;wave1[I][J].real()==0.0 && wave1[I][J ...

补 充:::::::
各位大侠 这是我根据一个实数的行列式求值改编的,但是运行不出想要的结果。
其中:wave1[m_system][m_system]是我想要求的行列式,m_system的值是200,
value_1_re 、value_1_im分别是所求的行列式的值的实部和虚部,
mulitiplication(x,y)和plus(x,y)是我自己编的两个小函数求出实数和实数相乘和相加的值。
大家可以直接忽略掉 for( M = I;wave1[I][J].real()==0.0 && wave1[I][J].imag()==0.0;M++)cout<<"M="< 最后的两个分号我已经改掉了,但是每次运行到这个循环的时候M的值不是说到200就停止了,而是一直不停的增加下去,跟个死循环差不多,跳不出去了,不知道哪里出问题了,各位 帮帮忙吧!先谢谢了。。。。。。
5楼2011-10-14 14:08:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见