24小时热门版块排行榜    

Znn3bq.jpeg
查看: 768  |  回复: 3

youyouxiong

新虫 (小有名气)

[求助] matlab高人帮看看这个程序

matlab菜鸟,对这东西真不懂啊,请教各位高手
怎样用matlab编了个计算程序,能运行,但是和计算公式有差异,高手帮看看问题在哪,附件是应该得到图线形式
timeg=1.65

n=95;                           %selected  points
voltage=[0.65637
0.34749
0
4.51737
9.65251
13.16602
18.30116
23.74517
28.26255
34.32432
38.49421
42.66409
45.90734
50.3861
54.55598
58.72587
61.9305
65.44402
68.33977
71.89189
75.40541
77.64479
79.26641
76.67954
74.13127
71.23552
73.4749
76.67954
79.88417
84.71042
86.94981
88.22394
85.67568
89.84556
87.60618
86.29344
84.40154
82.12355
85.98456
84.71042
86.64093
88.57143
90.15444
92.08494
91.42857
89.84556
87.25869
90.81081
91.11969
89.53668
88.57143
86.94981
85.01931
86.29344
88.57143
91.42857
88.22394
85.98456
83.43629
80.23166
77.64479
74.44015
70.92664
67.72201
63.51351
57.76062
53.89961
51.00386
46.5251
42.04633
35.63707
31.46718
25.98456
21.50579
15.71429
9.30502
6.10039
2.58687
0.34749
-0.30888
0
2.89575
4.51737
2.58687
0.65637
0.34749
0
0.34749
3.55212
8.03089
11.23552
15.09653
18.30116
22.77992
25.67568];

%integration
pressure=voltage;
timemm=0:timegn-1)*timeg;
integfunction=pressure;
pressure(1)=0;
for i=1n-1)
    xx=[timemm(i),timemm(i+1)];
    s=[i-1,i]
    ff=[integfunction(i)*sin(7426.067244*(i-s)),integfunction(i+1)*sin(7426.067244*(i+1-s-1))];
    pp=spline(xx,ff);
    int_pp=fnint(pp);
    pressure(i+1)=ppval(int_pp,[timemm(i),timemm(i+1)])*[-1;1];
end
for i=1n-1)
    pressure(i+1)=pressure(i)+pressure(i+1);
end        
plot(timemm,pressure);
grid on

1


回复此楼

» 猜你喜欢

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

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

youyouxiong

新虫 (小有名气)

正确的图应该是图2,图1是voltage
2楼2012-07-08 14:48:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

白衣太史

新虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-07-09 21:23:31
你的程序编制没有逻辑问题,要解决问题你的公式需要贴出来。
3楼2012-07-09 19:22:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Moon8808

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
xiegangmai: 金币+2, 谢谢参与! 2012-08-09 22:42:35
youyouxiong: 金币+5, 有帮助 2012-08-09 22:49:13
timeg=1.65;w=7426.067244;
timemm=0:timegn-1)*timeg;
for ii=1:n
     tt=timeg*(ii-1);
     integfunction(ii)=(2*sin((tt*w)/2)^2)/w*p(ii);
     K(ii)=sum(integfunction(1:ii));
end
plot(timemm,K)
4楼2012-08-09 22:37:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 youyouxiong 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 8/400 2026-05-24 11:58 by 晓晓爱翠翠
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 7/350 2026-05-24 10:20 by 光伏爱好者
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 08:01 by 9ps9vgkqva
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 pmo95bazuy 2026-05-23 7/350 2026-05-24 06:35 by fpo5ljpv91
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见