| 查看: 575 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】matlab模拟公式
|
|||
|
我想模拟附件中的公式，下面是我写的程序 clc close all; syms xi eta tau ; a=1; tau=1; c=1500; rho=1000; N=100; %t=[0:1:1]; x=0; xi=300; eta=300; p1=zeros(N,1); %采样点数 f=zeros(N,1); %xi eta粘滞 rho密度 tau 脉宽 omega=sqrt(4*a^2*c.^2-a^4*((xi+(4*eta/3))/rho)^2); f(t)=tau*exp(-((t/pi).^2))/sqrt(pi); s(t)=diff(f)-((xi+(4*eta/3))/rho)*diff(f,2); for t=1:N p1(t)=exp(-(a^2)*((xi+(4*eta/3))/rho)*t)*sin(omega*t-a*x)/(2*omega); p(t)=conv(p1(t),s(t)); end plot(p) 报错如下 ??? Error using ==> unknown Matrix dimensions must agree. Error in ==> viscosityl at 18 s(t)=diff(f)-((xi+(4*eta/3))/rho)*diff(f,2); 请高手帮忙改正，谢谢 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
woshilsh(金币+4,VIP+0):常来，吼吼 4-25 13:42
loulou85(金币+10,VIP+0): 4-25 15:43
woshilsh(金币+4,VIP+0):常来，吼吼 4-25 13:42
loulou85(金币+10,VIP+0): 4-25 15:43
|
>> a=1; tau=1; c=1500; rho=1000; N=100; t=linspace(0,1,100); x=0; xi=300; eta=300; p1=zeros(N,1); %采样点数 f=zeros(N,1); %xi eta粘滞 rho密度 tau 脉宽 omega=sqrt(4*a^2*c.^2-a^4*((xi+(4*eta/3))/rho)^2); f=tau*exp(-((t/pi).^2))/sqrt(pi);d1=diff(f);d2=diff(f,2); s=d1(1  end-1))-((xi+(4*eta/3))/rho)*d2;for t=1 N-2)p1(t)=exp(-(a^2)*((xi+(4*eta/3))/rho)*t)*sin(omega*t-a*x)/(2*omega); p(t)=p1(t)*s(t); end >> plot(p) >> |
2楼2009-04-25 12:36:47
3楼2009-04-25 15:49:08
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
4楼2009-04-25 19:06:52
5楼2009-04-26 11:07:36
forever_wc
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5124.4
- 散金: 120
- 红花: 2
- 帖子: 55
- 在线: 4210.7小时
- 虫号: 756344
- 注册: 2009-04-24
- 专业: 传热传质学
6楼2009-04-27 11:53:56
投票:
30