24小时热门版块排行榜    

Znn3bq.jpeg
查看: 718  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 药学求调剂 +3 RussHu 2026-04-12 4/200 2026-04-12 17:49 by 陈皮皮
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 计算机22408 281分,求调剂 +7 17715607211 2026-04-06 7/350 2026-04-12 00:43 by xuxiang
[考研] 269电子信息求调剂,可转专业 +11 独酌wl 2026-04-06 11/550 2026-04-11 11:12 by 逆水乘风
[考研] 调剂 +4 电气300求调剂不 2026-04-08 7/350 2026-04-11 10:44 by 紫曦紫棋
[考研] 085500求调剂材料 +10 易11122 2026-04-09 10/500 2026-04-11 10:39 by maddjdld
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 一志愿西交机械专硕求调剂 +8 求上岸的小王 2026-04-10 8/400 2026-04-10 15:09 by hemengdong
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 生物学学硕,初试351分,求调剂 +4 …~、王…~ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 调剂求助(生物与医药) +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
信息提示
请填处理意见