24小时热门版块排行榜    

Znn3bq.jpeg
查看: 723  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 307中医考研调剂 +4 于以采蘩 2026-04-14 4/200 2026-04-15 19:14 by AN流800
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 化工学硕294分,求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 290求调剂 +21 luoziheng 2026-04-10 23/1150 2026-04-14 15:49 by zs92450
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 电气专硕320求调剂 +6 小麻子111 2026-04-10 6/300 2026-04-12 10:54 by lemon6009
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 0859,337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
信息提示
请填处理意见