24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

当前主题已经存档。

loulou85

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 81
散金: 50
红花: 4
帖子: 737
在线: 140.9小时
虫号: 623844
注册: 2008-10-12
性别: GG
专业: 声学

[交流] 【求助】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);

请高手帮忙改正，谢谢

回复此楼

» 猜你喜欢

三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复

1楼2009-04-25 11:39:06

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

hitzhang

木虫 (正式写手)

仿真EPI: 1
应助: 0 (幼儿园)
贵宾: 2.15
金币: 1376.7
散金: 969
红花: 8
帖子: 862
在线: 226.4小时
虫号: 390575
注册: 2007-06-02
性别: GG
专业: 无机非金属类电介质与电解

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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)
>>

赞一下(13人)

回复此楼

2楼2009-04-25 12:36:47

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

loulou85

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 81
散金: 50
红花: 4
帖子: 737
在线: 140.9小时
虫号: 623844
注册: 2008-10-12
性别: GG
专业: 声学

t=linspace(0,1,100); 这里的100是秒吧，我改成t=linspace(0,0.001,10);咋报错??? Attempted to access s(9); index out of bounds because numel(s)=8.

Error in ==> viscosity2 at 22
p(t)=p1(t)*s(t);
呢