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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab问题，菜鸟求指教
81/1返回列表
查看: 1777  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

北尘心

金虫 (著名写手)

[求助] matlab问题，菜鸟求指教 已有2人参与

L=100;
x=1:100;
for cas = 1:120
    scr = ['C:\Users\asus\Desktop\' num2str(cas) '.txt'];
f=fopen('scr','r');
g=fscanf(f,'%f',[L,L]);
fclose(f);
id=0;
K=0:0.001:pi;
   for k=0:0.001:pi;
        id=id+1;
        ans=0;
        for i=1:L
            for j = 1:L
                ans = ans + g(i,j)*cos(k*(i-j));            
            end
        end
        ans =ans / (2*100);
        pair(id)=ans;
    end
   
    pair=-pair;
    plot(K,pair,'LineWidth',2);
    xlabel('k_m');
    ylabel('Gpair');
        hold on;   
    end
end


??? Error using ==> fscanf
Invalid file identifier.  Use fopen to generate a valid file identifier.
求指导，想知道哪里错了，怎么改！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-14 19:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asmeng

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
北尘心: 金币+5, ★有帮助 2014-03-22 01:07:17
ben_ladeng: 金币+1, 谢谢指教 2014-03-22 06:57:11
f=fopen('scr','r');改为f=fopen(scr,'r');去掉引号试试。
赞一下 回复此楼
2楼2014-03-14 19:33:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

北尘心

金虫 (著名写手)
引用回帖:
2楼: Originally posted by asmeng at 2014-03-14 19:33:41
f=fopen('scr','r');改为f=fopen(scr,'r');去掉引号试试。

试了，但还是不行
回复此楼
3楼2014-03-14 19:47:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hppdyx

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
北尘心: 金币+5 2014-03-22 01:07:25
ben_ladeng: 金币+1, 谢谢指教 2014-03-22 06:57:21
1. 你把scr = ['C:\Users\asus\Desktop\' num2str(cas) '.txt'];
改为scr = stact('C:\Users\asus\Desktop\'num2str(cas)'.txt');

2. 如楼上所说，把fopen函数里面的scr头上的引号去掉
赞一下 回复此楼
不以风骚惊天下，但求淫荡动世人
4楼2014-03-14 20:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hppdyx

木虫 (知名作家)
引用回帖:
4楼: Originally posted by hppdyx at 2014-03-14 20:00:58
1. 你把scr = ;
改为scr = stact('C:\Users\asus\Desktop\'num2str(cas)'.txt');

2. 如楼上所说，把fopen函数里面的scr头上的引号去掉

第一个里面的函数是strcat，写错了。。。不好意思
赞一下 回复此楼
不以风骚惊天下，但求淫荡动世人
5楼2014-03-14 20:02:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asmeng

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
北尘心: 金币+3, ★有帮助 2014-03-22 01:07:40
ben_ladeng: 金币+1, 谢谢指教 2014-03-22 06:57:34
你的数据是两列吧，错误提示是fscanf,应该是fscanf(f,'%f%f',[L,L]);另外，[L,L]应该是一个区间吧？[-L,L]?
赞一下 回复此楼
6楼2014-03-14 20:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asmeng

金虫 (正式写手)
引用回帖:
6楼: Originally posted by asmeng at 2014-03-14 20:28:52
你的数据是两列吧，错误提示是fscanf,应该是fscanf(f,'%f%f',);另外，应该是一个区间吧？?

[L,L]应该是表示读的数据的个数，L×L个，按列的顺序。
赞一下 回复此楼
7楼2014-03-14 20:34:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

北尘心

金虫 (著名写手)
引用回帖:
7楼: Originally posted by asmeng at 2014-03-14 20:34:14
应该是表示读的数据的个数，L×L个，按列的顺序。...

还是不行啊
回复此楼
8楼2014-03-22 08:50:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 北尘心 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭